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Journal of the 
HYDRAULICS DIVISION 


Proceedings of the American Society of Civil Engineers 


THE APPLICATION OF SEDIMENT-TRANSPORT MECHANICS 
TO STABLE-CHANNEL DESIGN 


Emmett M. Laursen,! A.M. ASCE 
(Proc. Paper 1034) 


SYNOPSIS 


The three requisites for a stable alluvial channel are explicitly stated and 
the role of sediment transport in each assessed. The similarity of the many 
sediment transport formulas is demonstrated and the general method of ap- 
plication to design illustrated, The use of these formulas as scaling rela- 
tionships between different channels is advocated, 


Not enough is known about the mechanics of sediment transportation for 
theory alone to be used with confidence in the design of an alluvial channel. 
However, enough is known for transport theory to be an aid tothe designer 
in the exercise of engineering judgment. In part this aid consists of dispel- 
ling confusion by demonstrating that design methods and criteria which are 
seemingly different are instead basically equivalent. 

By definition, a stable sediment-bearing canal in erodible material must 
not only convey the required water discharge, but must also transport the 
sediment load supplied, and must not scour or silt either the bed or the banks. 
These three factors neatly delineate the designer’s problem, the solution to 
which is almost single valued—“almost” only because there may be a narrow 
but finite range of flow conditions above which the banks will not silt and be- 
low which they will not erode. The designer’s problem is the selection of the 
“correct” combination of slope, size, and shape of channel which will result 
in a stable canal for a particular combination of water discharge and sedi- 
ment load (or range thereof) and the particular materials through which the 
canal must pass. The selection must be based on personal experience guided 
by such relationships as have been deduced from general experience with 
existing channels and with laboratory flumes. The difficulty in selecting the 


Note: Discussion open until January 1, 1957, Paper 1034 is part of the copyrighted 
Journal of the Hydraulics Division of the American Society of Civil Engineers, Vol. 
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correct combination stems from an always limited personal experience and 
from a lack of confidence in the guiding relationships. 

In order to design a stable channel, three relationships are needed, one 
for each of the factors in the previous definition—a flow equation, a sediment- 
transport equation, and a bank-erosion criterion. Although the flow condi- 
tions in the stable canal must satisfy all three of these relationships simul- 
taneously, for clarity in exposition they will be considered separately in the 
following discussion. 


The Flow Equation 


Of the three relationships needed for the design of a stable alluvial chan- 
nel, the best understood is the flow equation. Basically this relationship is 
the Chezy equation 


Q=CAVRS (1) 


In the use of this equation the designer must assign, arbitrarily, a numerical 
value to the discharge coefficient C . In order to aid his judgment, he will 
usually prefer to use some empirical formula such as that of Manning: 


This form of the flow equation was widely accepted because it was generally 
believed that the coefficient n was primarily if not entirely a function of the 
boundary material. Although it is now realized that n can be affected by 
several factors, a wealth of information has been accumulated over the years 
on the n values for existing channels covering a wide range of size, shape, 
alignment, bank and bed material, and condition. Because of this store of in- 
formation, n values for channels similar to that being designed can usually 
be found, The importance of this fact is self evident, but it is noteworthy 
here because the same sort of aid to judgment is needed for the other two 
relationships. 

Sediment transportation will affect the flow, and it will hence be reflected 
in the flow equation insofar as it affects the discharge coefficient. Dunes are 
almost always associated with the transport of sediment, and the resistance 
coefficient, which is the inverse of discharge coefficient, is affected by the 
resulting increase in boundary roughness. That sediment transport may in- 
fluence the resistance coefficient in other ways has been claimed(1) but has 
not been proved unequivocally. No matter what the causal relation, however, 
the rate of sediment transport is one of the several criteria which should be 
considered in appraising the similarity of channels. 

Two attempts to relate the resistance coefficient to sediment transport 
should be mentioned. Einstein and Barbarossa(2) divided the bed resistance 
into two parts, that due to the sand roughness and that due to bar, or dune, 
roughness. To evaluate the sand roughness they used Strickler’s formula 


n=VK, /29.3 


4 
= 
Fi 
ime 


ASCE LAURSEN 1034-3 be 
where K, is the diameter (measured in feet) of the particle size coarser : 
than 65% of the sediment. After eliminating bank effects in narrow channels, 

the effect of the sand roughness was subtracted from the total resistance and 

this remainder was related to a modification of Einstein’s y function for the 
transportation of bed load. | 

Barton and Lin(3) have proposed a relationship between the Chezy coef- : 
ficient, the Reynolds number of the flow, and the ratio between the sediment 
size and the hydraulic radius. Although sediment transport does not enter the 
relationship explicitly, the data were obtained from flows transporting sedi- 
ment. 

Both of these proposals are so new that there has not been sufficient time 
to evaluate their usefulness or their weaknesses. Both of them, however, by- 
pass the problem of predicting dune size, which is probably the major effect 
of sediment transport on resistance to flow. Anderson(4) recently obtained 
a solution for the length of dune as a function of the Froude number and pre- 
sented a limited amount of data corroborating his mathematical solution. His 
solution had a limit of a plane bed where the dune length became infinite. In 
addition he presented one example of a plane bed at a Froude number of 1.6. 
In measurements at the Iowa Institute, however, a plane bed (i.e., without 
dunes) has been obtained at a Froude number as low as 0.85. The sand in this 
experiment was 0.1 mm in diameter, finer than the sand of Anderson’s exam- 
ple. Of course, it is possible for dunes to disappear because the dune height 
reduces to zero as well as because the dune length increases in infinity. The 
height of the dune, which was not and cannot be a part of Anderson’s analysis, 
is fully as important as the length. Very likely it is more a function of the 
particle size and the Reynolds number of the flow than of the Froude number. 

The empirical equations describing stable channels which have been de- * 
veloped in India(5) include relationships which are in essence evaluations of . ; 
the Chezy C , just as the Manning equation is. These flow relationships 4 
naturally are subject to the same criticism as is often made of the other q : 
Indian criteria—that since they are based on data from Indian canals they are ns 
tailored to fit Indian conditions. This criticism, although valid, does not is 
mean that due consideration should not be given to modifying the relation- 
ships to fit the conditions in other areas of the world. 

Whatever flow equations and whatever numerical value of a resistance co- 
efficient the designer chooses, he will arbitrarily, though to the best of his 
judgment, have placed a restriction on the physical character of the channel 
eventually to be selected. Before the final decision two more restrictions in- 
volving the same type of evaluation must be explicitly stated. ag 


The Sediment- Transport Equation 


If the flow entering the canal is not clear, but carries with it a sediment 2 
load, the capacity of the channel must exactly equal the rate at which sediment . 
is supplied to the channel for regime conditions to prevail. This requirement, oS 
then, places a second restriction on the physical characteristics of the chan- a 
nel to be selected. Some of the many sediment-transport equations for bed 2 
load that have been proposed are listed in Table I, This table is not com- 
prehensive nor is it intended to rate or judge the many expressions that can 
be found in the literature, In fact, the problem of picking a sediment-trans- 
port equation is an unhappy one that must be faced by each individual designer. 
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All of the equations are slightly different in form; all are based on limited 
data, either entirely or largely from the laboratory; each fits some data bet- 
ter than the others; none fits all the data (which is why there are so many); 
none has gained the wide acceptance shared by the relatively few flow equa- 
tions. If unable to decide with confidence on one “best” equation, the designer 
may use more than one, or even all, in the same manner as the hydrologist is 
wont to do with the many flood formulas at hand. 

In Table I, it should be noted that the equations suggested by Einstein(12) 
and Kalinske(13) in their original papers have not been used, but instead the 
modifications suggested in the sediment chapter of Engineering Hydraulics. 
In these modifications the original parameters are used, but the forms of the 
equations are changed to fit the data for the high rates of sediment transport. 
That these modified equations are not strikingly different from the other em- 
pirical sediment-transport equations should not be surprising, because in 
their development the assumption is made that the velocity at the edge of the 
laminar sublayer vz (whether or not actually existent) can be substituted 
for the velocity at the level of the moving grain of sediment vg . Since vj 
is, as usual, equated to 11.6 /7/p and bears little relation to Vg except 
that both are small and both are at levels near the bed, this is tantamount to 
assuming that qs is a function of T. 

In order to compare the different transport equations more easily, they 
have been put into a form indicating the dependence of qs on the mean veloc- 
ity V and the depth of flow y by neglecting the critical shear T,. and using 
the Manning equation to eliminate the slope. The similarity between not only 
these but all the commonly used sediment-transport equations can thus be 
made apparent. One might also suspect, along with Gilbert, (14) that the ex- 
ponent as well as the coefficients in the transport equation should be variable. 

If a sizable bed load is to be expected, the simplified forms can be used to 
estimate the transport characteristics of the canal. However, since it is dif- 
ficult to measure the bed load or to predict the bed-material load entering 
through the headworks of proposed canal, a further modification of the forms 
may be helpful. For example, if the first equation in Table I is written for 
the sediment discharge ratio, and then the ratios of the various parameters 
are taken for the proposed and an existing “similar” channel, relative rather 
than absolute values can be estimated: 


qs/q (existing) B,(e) n*(e) y 5/3 (e) 


It is interesting to note that Chezy proposed his flow equation in this ratio 
form,(15) and that the accumulation of n values together with descriptions of 
channels is essentially only a scmewhat more sophisticated manner of han- 
dling the problem, 

A further rearrangement of the simplified transport equations will put 
them in the form of the Kennedy equation and the variations thereof as shown 
in Table IL The similarity between the modified transport equations and the 
Kennedy type of equations indicates that the physical basis for the latter equa- 
tions is the requirement for sediment-transport capacity. That only Inglis 
has considered this factor explicitly is probably due to the fact that the equa- 
tions were derived from data of canal systems in which the sediment- 
discharge ratio did not vary greatly. 
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On the basis of personal judgment, aided to some extent by the sediment- 
transport equations of one form or another, the requirement for sediment 
capacity must be satisfied, just as the requirement for flow capacity was 
satisfied, and a second restriction on the combination of size, shape, and 
slope of channel is imposed. 


The Bank-Erosion Criterion 


If the banks of the canal were vertical and non-erodible, the final selection 
of channel characteristics satisfying only the first two requirements could be 
made either arbitrarily or on the basis of economics. An alluvial channel 
will have sloping banks of erodible material. Flow conditions along the banks, 
therefore, must be such that sediment will not deposit on the banks and that 
bank material will not be scoured. The first of these requirements will 
govern for canals designed for a minimum slope; the second for canals de- 
signed for minimum area. Satisfactory solution of the non-silting problem 
requires an understanding of sediment-transport mechanics in much greater 
detail than is presently available, especially of the suspended-sediment load. 
Knowledge concerning the non-scouring problem is somewhat more advanced, 
but the bank-erosion criterion is still the least satisfactory of the three re- 
lationships needed for design. 

In the case of banks composed of granular material (e.g., sands and 
gravels) erosion will occur if the forces tending to move the particles are 
sufficiently large to overcome the resisting force of gravity. In the case of 
banks composed of cohesive material (e.g., clays) the force of cohesion will 
also resist movement and in fact may be much greater than the resisting 
force of gravity. 

In associating the flow conditions in the canal to the bank erodibility the 
canal shape enters the problem on both sides of the relationship. The slope 
of the banks will obviously affect the susceptibility of the banks to erosion 
and the cross-sectional shape will affect the distribution of the boundary 
shear. 

Lane and Carlson(20) have presented the relationship for the ratio of the 
critical tractive force on a bank slope Tcs to the critical tractive force Tep 
as a function of the angle of the side slope and the angle of repose of the ma- 
terial. That is, 


tan@d 
K= tap = cos (5) 


where ¢ is the angle of the side slope, and @ is the angle of repose. Other 
studies under the direction of Lane(21) have attempted to obtain the distribu- 
’ tion of shear in a trapezoidal section. This knowledge of the effect of cross- 
sectional shape allows data or equations on the critical tractive force for bed 
material to be transferred to the banks. 
The two formulas for critical tractive force most widely quoted are those 
of White(22) 


Teb = 12 D (6) 
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and of Shields(10) 


Teb = (Sg!) D D/8) (7) 


Any of the relationships for tT, of the equations in Table I can also be used— 
at the discretion of the designer. It should be noted, however, that these re- 
lationships may result in a low rate of movement rather than zero movement. 
A criterion based on zero movement is, perhaps, unduly restrictive, because 
a small amount of scour can usually be tolerated and because suspended sedi- 
ment will tend to add cohesive material to the banks. Application of the rela- 
tionship from the first equation in Table I will result in values of permissible 
canal velocities approximately the same as given by Fortier and Scobey; (23) 
application of the relationships of White and Shields will result in smaller 
permissible velocities. 

Very little is known concerning the resistance to movement of clays—ex- 
cept that some clays can stand up under very high boundary shear. One might 
speculate that the resistance to boundary shear should be related to the in- 
ternal shearing characteristics of the clay mass in an undisturbed saturated 
condition. For naturally deposited material, which is not likely to be com- 
pletely homogeneous, the weakest part of the mass would govern. Erosion of 
seams or pockets of weak material would, of course, expose the adjacent 
stronger material to greater forces. The only published work concerning the 
resistance of cohesive materials is contained in Fortier and Scobey, and even 
there insufficient information is given about the character of the material to 
serve as more than a rough guide for the designer. 

This last relationship, together with the first and second, fixes the width, 
depth, and slope of the channel. No method is available to determine the best 
geometrical shape of the cross section—except examination of other channels. 
In fact, canals are usually dug to a shape easily handled by construction ma- 
chinery, and then allowed to modify that shape by a moderate amount of scour 
and fill. 


CONCLUSIONS 


Three relationships are needed to specify a stable alluvial channel, The 
three relationships may be those directly stating the three requirements for 
the stable canal as presented herein, or they may be any three independent 
forms derivable therefrom. The three independent regime-theory criteria 
implicitly contain these same factors. 

The difficulties of designing a stable canal lie not in using the general re- 
lationships, but in assigning numerical values to coefficients and even in 
choosing a formal equation. The essence of design procedure is shown in 
Fig. 1. Choice of the equation form and numerical values to describe the 
sediment-transport capacity will fix a relationship between the mean velocity 
of flow and the depth of flow as indicated by the heavy curve in the schema- 
tized diagram. Choice of the equation form and numerical values to describe 
the flow capacity will result in a family of curves relating velocity and depth 
for various canal slopes as indicated by the dashed curves. The dot-dash 
curve of constant tractive force describing the bank competence where it in- 
tersects the transport-capacity curve will specify one of the family of flow- 
capacity curves, Velocity, depth, and slope are thereby fixed, and, of course, 
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the width can be obtained from the continuity principle. 

The interrelationship shown above constitutes the general application of 
sediment-transport mechanics to stable-canal design. Only a complete, com- 
prehensive, rigorous theory, however, could substitute for a full scale model 
canal. It is by observation of similar canals that the designer can assign the 
needed numerical values, and sediment-transport mechanics can aid in as- 
sessing the degree of similarity of different canals. The parameters of the 
various equations can be used as comparative criteria and the equations 
themselves as approximate scaling relationships. If the proposed and exist- 
ing conditions are treated as ratios, the discrepancies between absolute 
values as computed by the equations should be greatly reduced. 

Practically speaking, few if any canals are completely stable. In planning 
remedial and maintenance work the principles of sediment-transport 
mechanics can be used in the same manner as in design, the principal dif- 
ference being that a canal almost completely similar is then readily available 
for observation. 
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Fig.! Requisites for stability of an alluvial channel 


| 
| | 
| 
\ 
b 
i 
| 
a 
q 


© 
= 
Po 
a 
3 
< 


HY 4 


1034-10 


a 89 ay =Sb 
Lg = = Sb fii] 
b = Sb fo] spiatus 
= = Sb [6] S3m 
gu &g ,¢(0 YY -S = Sb 
'g= (2-2)2 'y =Sb [9] (qnowys) skogng 


: 
| 
| 
4 
. 
: 
4 
= 
| 
| 
| 
| 
| 
a 
4 
i 
: 


LAURSEN 


Duboys (Straub) 
Schoklitsch 
Meyer-Peter 
WES 

Shields 
Brown-Einstein 
Brown-Kalinske 
Kennedy [6] 
Lacey 
Blench [Ig] 
Inglis [i9] 


m’ =(m+3)/(6m- 3) 


iA2m-1) 
Ss 


wil4 


q 


1034-11 


| 
Table I. Similarity of Bed-Load and Kennedy Formulas oa 

| 

‘| 

1/3 

2n q 

4 

U4 

| y 

| 3/10 

| 

V 

| = 

vb 

4 

V 
* C9 

| 

4 


a! 
| 
| 
a 
4 
4 
4 
4 
4 


Paper 1035 HY 4 


Journal of the 
HYDRAULICS DIVISION 
Proceedings of the American Society of Civil Engineers 


MEASURING EVAPOTRANSPIRATION FROM ATMOSPHERIC DATA 


George S. Benton,! A.M. ASCE and Jack Dominitz2 
(Proc. Paper 1035) 


SYNOPSIS 


A method of evaluating evapotranspiration is presented which is based on 
the mass balance of water vapor in the atmosphere. The method is applied 
to various regions of the North American continent for the calendar year 
1949 and the results are compared with hydrologic data and with an empirical 
method of estimating evapotranspiration. 

It is shown that although accurate values of evapotranspiration can be ob- 
tained for large areas such as the entire continent, the accuracy of the meth- 
od is reduced as the size of the area under consideration is decreased. The 
potentialities and limitations of the proposed method are evaluated, 


INTRODUCTION 


Evapotranspiration—the loss of water from natural surfaces by evapora- 
tion and transpiration—has resisted measurement more strongly and is con- 
sequently less well known than most hydrologic parameters. The addition of 
moisture to the land by precipitation and its depletion through surface runoff 
may be measured, but the rate at which the soil moisture supply is lost to the 
atmosphere under natural conditions cannot be measured directly by existing 
instruments. 

Yet, a quantitative understanding of the process of evapotranspiration is of 
considerable importance to the meteorologist and the hydrologist. The mete- 
orologist would like more information concerning the water balance of the at- 
mosphere, as well as the turbulent exchange processes that take place in the 
microclimatic layer. On the other hand, the hydrologist would like to be able 
to maintain a balance sheet of soil moisture content to assist him in evaluat- 
ing drought hazard or irrigation requirements, or to improve methods of 
river forecasting or reservoir management. Evapotranspiration is also of 


Note: Discussion open until January 1, 1957. Paper 1035 is part of the copyrighted 
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importance to the military engineer, for off-the-road trafficability is strong- 
ly affected by the water content of the soil. 

Since instruments do not yet exist for measuring actual evapotranspiration 
under field conditions, scientists have turned to expressing evapotranspira- 
tion from natural regions by various crude formulae. Many such formulae 
have been developed. Some are based on simplified physical theory (e.g., 
Richardson and Montgomery(9))3 but most are purely empirical (e.g., 
Meyer;(7) Lioyd(6)), Naturally, one of the greatest difficulties in devising 
and using such formulae is the lack of objective evapotranspiration measure- 
ments to use as control data. Tank estimates of evapotranspiration and soil 
moisture determinations in experimental areas provide the bulk of such data, 
although in many regions the mean annual evapotranspiration can be computed 
with reasonable accuracy from precipitation and runoff measurements. 

A recent empirical formula for evapotranspiration which has received 
widespread attention was devised by Thornthwaite.(11) Thornthwaite devel- 
oped a formula for potential evapotranspiration, which he defined as the maxi- 
mum amount of water that could be lost from a region under a given set of 
climatic conditions. He also developed a simple water budget for estimating 
actual evapotranspiration from potential evapotranspiration. Thornthwaite’s 
method generally yields estimates of evapotranspiration which are reason- 
able in their seasonal distribution and annual total (Carnahan and Benton). (3) 
However, such empirical estimates cannot be properly evaluated until accur- 
ate measurements of evapotranspiration become available for natural re- 
gions. 

One possible method of determining evapotranspiration, which has not 
heretofore been adequately explored, involves an evaluation of the atmos- 
pheric water balance over the continents. The latter subject has been dis- 
cussed by many authors, but little quantitative information was available until 
Holzman(5) presented his analysis of the relationship between the air mass 
and the hydrologic cycles. The concepts advanced in this paper were later 
extended and developed by Benton, Blackburn, and Snead.(1) In recent years 
our understanding of the atmospheric water balance has been greatly in- 
creased by the extensive use of electronic computers in meteorological re- 
search. As an outgrowth of this program, the flow of water vapor over the 
North American continent was analyzed and discussed by Benton and 
Estoque. (2) 

The availability of extensive meteorological data and the ability to process 
this information numerically in an efficient manner, led to a growing realiza- 
tion that our knowledge of atmospheric vapor flow might be sufficient to allow 
the direct computation of evapotranspiration from a statement of the mass 
balance of atmospheric water, developed in a manner analogous to the fa- 
miliar hydrologic equation. Possibly the first attempt to evaluate evapotrans- 
piration in this manner was presented by Moller,(8) who used a sparse net- 
work of meteorological stations over northwest Europe. A more extensive 
study of evapotranspiration over the North American continent was presented 
by Benton and Estoque. (2) 

The investigations of Benton and Estoque showed that for land areas of 
continental extent, estimates of evapotranspiration obtained from atmospheric 
data are at least as accurate as those which can be obtained from hydrologic 
estimates or from empirical relationships. The major objective of the 


3. Numbers in parentheses refer to references at the end of the paper. 
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present study is to evaluate the extent to which atmospheric data can be suc- 
cessfully used to determine evapotranspiration from smaller areas. The at- 
mospheric balance equation is derived in the following section, and applica- 
tions of this equation to specific areas of various sizes are then presented. 
Computed evapotranspiration is compared with estimates obtained using 
Thornthwaite’s empirical procedure, as well as with annual evapotranspira- 
tion totals obtained from precipitation and runoff records. 

It is shown that although the procedure for evaluating evapotranspiration 
from atmospheric data is physically sound, practical difficulties arise. 
These difficulties, related to problems of statistical sampling, are compara- 
tively unimportant when'large areas and long time intervals are considered 
but become more critical as the size of the area and time interval is reduced. 
Nevertheless, the atmospheric method offers a fundamentally simple ap- 
proach to the evaluation of evapotranspiration which may become of increas- 
ing practical value as meteorological observations and techniques of elec- 
tronic data processing are improved. 


Development of the Atmospheric Balance Equation 


The conventional equation, expressing the basic relationship between 
evapotranspiration and other hydrologic parameters for any time interval, is 
the hydrologic equation: 


P=E+R+AS 


(inches) (1) 


Here P, E, R, and AS are precipitation, evapotranspiration, net runoff, and 
change in storage of surface and sub-surface water all expressed in inches 
depth, respectively. Ground water flow into or out of the region can usually 
be neglected if a large area is considered. Unfortunately, equation (1) cannot 
be solved for E for periods less than a few years in length, since AS is not 
negligible for short time periods. 

A similar mass balance equation may be developed for the atmosphere, 
expressing the relationship between evapotranspiration and meteorological 
parameters. The following expression is obtained for an arbitrary time in- 
terval and region: 


E=F+P+AN (gm/em*) (2) 


Here F is the net outflow of atmospheric water vapor (F = outflow minus in- 
flow); E is evapotranspiration; P is precipitation; and AN is the change in the 
precipitable water content of the atmosphere, Two terms, AN and F, require 
further explanation. 

The change in precipitable water content, 4 N, is defined as the difference 
between the atmospheric water content at the end of the time period and that 
at the beginning. It is analogous to the change in storage term in the hydro- 
logic equation, although fortunately the moisture holding capacity of the at- 
mosphere is so low that it has none of the numerical importance of the latter 
term, Precipitable water, N, is defined as 


(gm/cm”) (3) 
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Here °F is the specific humidity, in gm/kg; p is the pressure in millibars 


which varies from age to its surface value, pp; and g is the acceleration of 
gravity, in cm/ sec2, Over an arbitrary time interval T : 


AN = N,-N, (gm/em”) (4) 


where Ng is the precipitable water at the end of the interval, and Ny is the 
precipitable water at the beginning. 

The net outflow of water, F, from the arbitrary area is more difficult to 
evaluate. Let the vector w be the mean transfer of water vapor at an arbi- 
trary point in the atmosphere. Then: 


W= (4) dt (gm/cm-mb-sec) (5) : 


Here c is the horizontal wind vector, in cm/sec, expressed as a function of 
time t. Note that w is the flux of vapor in gm/sec through a cross-sectional 
area one centimeter in width and one millibar in thickness. The mean trans- 
fer of water vapor integrated vertically through the atmosphere is: 


w=/ wdp (gm/cm-sec) (6) 


where p, is now the mean surface pressure, Finally the net outflow of water 
vapor, F, from an arbitrary area A is given by: 


(em/em®) 


where n is a vector of unit length normal to an infinitely small segment of 
the perimenter, dA , and pointing out of the area; and the line integral is = 
evaluated around the perimeter of the area. Note that the transfer of liquid Fe 
water is neglected, It can be shown that the resulting error is negligible in 
all save extremely unusual local situations. 

Substituting equation (7) into (2): 


E=P+AN+ fa-war (gm/em’) 


All of the variables on the right hand side of equation (8) can be evaluated 
provided a rawinsonde or radiosonde network is available to provide esti- it. 
mates of wind velocity and humidity through the atmosphere. Such a network am 
is operated by U. S. Weather Bureau and Armed Services personnel, and par- 
ticipating weather stations are scattered throughout the United States and 
North America. Equation (8) thus provides a physically sound method of 
computing evapotranspiration from natural areas. Furthermore, if it is 
possible to evaluate the various terms in the atmospheric balance equation, 
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the change in storage of moisture in the soil could also be determined by 


equation (1). 


termined with considerable accuracy if satisfactory results are to be ob- 


plication of equation (8) to actual data. 


Description of Data 


Twenty-two radiosonde stations were utilized in the study, as listed in 


and 500 mb. The resulting data were entered on punchcards. 
Geostrophic winds were used in evaluating the mean transfer of water 


affected by local topography and by small scale turbulence since they are 


tered into the punch-cards along with pressure and humidity data. 


of water vapor transfer were summarized and averaged for each station- 


ponent relationship: 


1035-5 


substituting the appropriate value for evapotranspiration in the hydrologic 


It must be pointed out that although equation (8) is almost exact, and in- 
volves no hazardous assumptions, other difficulties do arise. The net inflow 
or outflow of water vapor through the atmosphere is for most areas and time 
intervals a small difference between two large numbers, which must be de- 


tained. This limits the size of the area and the length of the time interval 
which can efficiently be considered. The following sections describe the ap- 


The calendar year 1949 was chosen for analysis because the necessary 
meteorological data were readily obtainable from the research project with 
which this study was associated. Monthly evapotranspiration was evaluated 
using equation (8) for each of four sub-continental areas, represented by four 
polygons with meteorological stations at each vertex. These polygons were 
chosen to correspond approximately to the following regions: Canada, the 
central United States, north-central United States, and the Ohio River Basin. 


Table 1. At each station measurements of the elevation of constant pressure 
surfaces and of humidity were recorded twice a day, 0300Z and 1500Z (i.e., 
Greenwich Civil Time), at each of three levels in the atmosphere: 850, 700, 


vapor, W, at each station. The geostrophic wind is a theoretical wind which 
can be computed from the slope of constant pressure surfaces even in the ab- 
sence of actual wind observations. Geostrophic winds are comparatively un- 


computed from the average slope of the pressure surface over a large region, 
often several hundred miles in width. For vhese and other reasons, they are 
generally used in computing atmospheric transfers of energy, momentum, 
and other physical quantities. In the present study, geostrophic winds were 
computed graphically from maps routinely prepared by the Weather Bureau- 
Air Force-Navy Analysis Center in Washington, D. C. These winds were en- 


Once the basic atmospheric data were entered into the punch-cards, the 
following extensions were performed. For each observation, the eastward 
and northward components of the geostrophic wind velocity were evaluated 
and each was multiplied by the water vapor content. These two components 


month at each level to give w (equation (5)), and the integrated water vapor 
transfer at each station, W, was then evaluated graphically from the com- 


Ww w dp tf dp+}f dp (gm/cm-sec) (9) 
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TABLE 1 


RADIOSONDE STATIONS USED IN THIS STUDY 


Sta. Lat. Long. 2 


Station No. North West 
Miami, Fla. 01 25° 40 FF 
New Orleans, La. 02 30 00 9 16 
Brownsville, Texas 03 25 55 7 26 
El Paso, Texas 04 31 48 106 245 
Little Rock, Ark. 06 34 44 92 «(214 
Atlanta, Georgia 07 33 39 84 
Greensboro, N. C. 08 36 05 79 S77 
Portland, Me. 10 43 39 70 «(19% 
Buffalo, N. Y. 1] 42 56 78 44m 
St. Cloud, Minn, 13 45 35 94 ~«Ol 
Lander, Wyo. 14 42 48 108 43 
Medford, Ore. 15 42 23 122 525 
Santa Maria, Calif. 16 34 56 120 £255 
Port Hardy, B. C. 19 50 41 127 2 
Edmonton, Alta. 20 53 34 113 
Churchill, Man. 21 56 21 94 468 
Stephenville, Newfoundland 23 48 32 58 33m 
Upper Frobisher, Baffin Is. 24 63 44 68 335 
Yakutat, Alaska 25 59 31 139 40) 
McGrath, Alaska 27 62 56 155 37§ 
Barrow, Alaska 28 71 20 156 245 
Coppermine, Dist. of Mack. 29 67 47 115 158 

| 
— 
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4. Calculations of the change in storage for the year 1949 for ten sample 
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where x and y are directed towards the east and north respectively. Vapor 
transfer above 400 mb was neglected, One possible error in this evaluation 
involves the layer between 850 mb (the lowest level of observation) and the 
surface, Previous work has shown that the maximum flux of vapor may oc- 
cur in this layer. (2) 

As a last step in computing the net vapor flux out of each region, the ap- 
propriate line integral of the component of W normal to the boundary of each 
polygon was evaluated as indicated in equation (7). A linear variation of the 
integrated water vapor transfer was assumed between stations on the perime- 
ter, although such an assumption can lead locally to non-negligible errors. 
Finally, the remaining quantities required for the determination of evapo- 
transpiration were computed, The precipitable water content of the atmos- 
phere at the beginning and end of each station month was estimated graphical- 
ly in a manner similar to that of vapor transfer, and the mean change in 
precipitable water for each region was then found by averaging weighted sta- 
tion values using the method developed by Thiessen.(10) Thiessen’s method 
was also used in evaluating precipitation from station measurements. With 
net vapor flux, change in precipitable water, and mean precipitation deter- 
mined, evapotranspiration was computed for each month for each of the four 
regions under consideration using the atmospheric balance equation. 

As indicated previously, evapotranspiration was also estimated by two 
other methods. The first of these was Thornthwaite’s empirical equation. 
Although Gilbert and Van Bavel (4) have shown that Thornthwaite’s equation 
yields somewhat better results when daily temperatures and precipitation 
are used, the present study was limited for simplicity to the use of monthly 
data in accordance with the procedure outlined in the original reference, (11) 
Thirty-seven stations scattered through North America, were employed in 
order to get adequate coverage. Thiessen’s method was again used to find 
the mean value of evapotranspiration from station values. 

A further check on the values of the total evapotranspiration for any year 
may be obtained by evaluating the runoff, neglecting the change in storage, 
after equation (1). It is anticipated that such estimates would usually be with- 
in two inches of the real values.4 River discharge data were compiled from 
U. S. Geological Survey records.(12) Since only United States records were 
readily available, three polygons (central United States, north-central United 
States, and the Ohio River Basin) were checked. This involved the examina- 
tion of sixty natural watersheds located within these polygons. 


Presentation of Results 
North American Continent 


As noted previously, the field of water vapor flux over the North American 
continent has been described in detail by Benton and Estoque.(2) As part of 
their study, the water balance of the continent was computed and measure- 
ments of evapotranspiration using atmospheric data were compared with hy- 
drologic estimates and with precipitation. For the convenience of the reader, 
the monthly distribution of evapotranspiration for the continent is reproduced 


stations, using a method developed by Thornthwaite,(11) indicates that 
two inches is usually a conservative estimate for large areas. 
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in Figure 1, For a discussion, the reader is referred to the original refer- 
ence. 


Canada 


The distribution of the monthly evapotranspiration for the area as illus- 
trated in Figure 2 is shown in Figure 3. A sample tabulation of the various 
quantities in equation (8) is shown in Table 2. The relationship between the 
evapotranspiration evaluated from atmospheric data and that determined from 
Thornthwaite’s empirical expression is very similar to that found for the 
North American continent, Values computed by use of the atmospheric bal- 
ance method are higher in the winter; reach a maximum later in the year; 
and result in a larger annual total. These differences may in fact be caused 
by certain assumptions made in the empirical formula. Thus, Thornthwaite 
assumed that evapotranspiration ceases whenever the mean monthly tempera- 
ture falls below freezing, and that potential evapotranspiration is fully real- 
ized as long as any water remains in storage in the soil. Neither of these 
assumptions is realistic. They would result in errors similar to the differ- 
ence found in Figure 3. 

Due to the unavailability of river discharge data for most of northern 
Canada, it is difficult to estimate the annual evapotranspiration using runoff 
measurements. It may be noted, however, that the annual evapotranspiration 
determined from atmospheric transfers of water vapor was 17.0 inches, 
whereas precipitation totalled 21.7 inches. This would require runoff for the 
year to total only 4.7 inches, if change in storage were neglected. Also, an 
increase in storage might be expected since precipitation in the last few 
months of 1949 exceeded that in corresponding months of 1948. Assuming the 
value of the precipitation to be correct, it becomes evident that evapotrans- 
piration may be slightly overestimated. This apparent discrepancy and its 
cause will be investigated in a later section. 


Central United States 


The monthly distribution of evapotranspiration for this area (Figure 2) is 
shown in Figure 4. In general, the relationship between the curves obtained 
from the present method and from Thornthwaite’s empirical formula remains 
the same, although the atmospheric balance gives consistently higher monthly 
values. A check on the annual total may be made since the runoff is known. 
For the year, precipitation amounted to 38.0 inches and runoff to 10.9 inches. 
Annual evapotranspiration should therefore be 27.1 inches plus a correction 
due to changes in storage between the beginning and end of 1949. This change 
was probably negative for this area primarily because November, 1949, was 
a very dry month (U. S. Weather Bureau).(13) A reasonable estimate of 
evapotranspiration might therefore be between 27 and 29 inches for the year. 
The computed value of 34.6 is high, and cannot be explained by inadequacies 
in the precipitation and runoff data. Nevertheless, the seasonal distribution 
of evapotranspiration is good. 

It may be noted in passing that the empirical method yields an estimate of 
the annual evapotranspiration equal to 24.7 inches. This figure is apparently 
too low. The smooth curves in Figure 4 should not, therefore, be regarded 
as accurate. The correct evapotranspiration probably lies somewhere be- 
tween the two sets of estimates. 
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North-central United States 


The area chosen is presented in Figure 5 and the computed monthly dis- 
tribution of evapotranspiration is compared with that of the empirical method 
in Figure 6. Although both curves display the same general characteristics, 
the difference found in the examples presented above are again apparent. 
Thus, the empirical method yields lower values and an earlier maximum. 

The annual evapotranspiration may be roughly estimated from the hydro- 
logic equation. Precipitation for this area was found to be 44.8 inches while 
runoff was 16.2 inches. This results in an annual total of approximately 29 
inches plus a correction due to change in storage. It is estimated that this 
change in storage was negative since a large portion of this polygon (the Ohio 
River Basin) experienced a fairly large negative change.(3) The annual 
evapotranspiration may, therefore, be between 29 and 31 inches. The com- 
puted value of 32.8 inches may be slightly too high; however, it is probably 
more accurate than the low estimate obtained by use of Thornthwaite’s 
method (25.4 inches). It is evident that although the individual monthly values 
of evapotranspiration obtained from the atmospheric moisture balance may 
not be too precise, the seasonal distribution and the annual total are reason- 
ably adequate. 


Ohio River Basin 


A volygon representing the Ohio River Basin (Figure 7) is next considered, 
This area was chosen primarily because much pertinent work had been done 
by Carnahan and Benton(3) in evaluating hydrologic parameters, including 
evapotranspiration, for this area and year. The comparison between the two 
methods (atmospheric transfer and Thornthwaite) is presented in. Figure 8. 
The resemblance between the two sets of data is poor and the trend of the 
values obtained from the proposed method during the first six months seems 
to be opposite to what is expected. The maximum, although of the proper 
magnitude, occurs late in the summer season. The annual evapotranspiration, 
however, seems to be of the proper magnitude, 28.5 inches, compared to the 
“best estimate” of 28.2 inches, found by Carnahan and Benton employing many 
different methods. Therefore, although the monthly data and seasonal distri- 
bution are unsatisfactory the annual value seems to be about correct. 

Whether the discrepancies in the monthly data are caused by the fact that this 
is the smallest of all areas considered (826,000 km2), will be discussed in a 
later section. It is interesting to note, however, that this region is not very 
different in location or size from the north-central United States, for which 
much better results were obtained. This might suggest that some systematic 
error, rather than a random sampling error, is involved. 


Statistical Analysis of the Atmospheric Balance Equation 


The atmospheric balance method as developed earlier, is based on sound 
physical principles and theoretically should give correct values for the rate 
of evapotranspiration from any selected area. From the limited number of 
examples furnished in the previous section, however, it is evident tha cer- 
tain difficulties are encountered in practice which must be removed before 
this method can be applied to sub-continental regions. To clarify these dif- 
ficulties, statistical analysis of an idealized system has been used to deter- 
mine the factors affecting the accuracy of the atmospheric balance equation. 
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Such analysis was limited to the divergence term since the possibility of er- 
ror is greatest in evaluating the net flux of water vapor and all errors in 
computing the other terms are negligible in comparison. 

Consider an idealized area, A, on whose perimeter, P', n stations are 
located with equal spacing, Each station reports N observations at 7 levels, 
each representing equal pressure increments through the atmosphere. Such 
a model is illustrated in Figure 9. 

Let the standard error of the mean value of any component of w (the mean 
atmospheric transfer of water vapor at any level and station) be given by s. 
Integrating the mean transfer through the height of the atmosphere: 


W= w. h, (gm/cm-sec) (10) 
i 


where hj is the height, in millibars, of the atmospheric layer represented by 
the value of the mean transfer, wj. Since all levels are spaced equally, hj is 
a constant and will hereafter be denoted as h. 

The standard error of estimate of any component of W (the integrated va- 
por transfer) is made up of two parts: error due to sampling, since only a 
finite number of levels are used through the atmosphere; and error due to 
inaccuracies in the computed values of w. For averages over a moderate 
period of time (e.g., a week or more), w varies smoothly through the atmos- 
phere; and hence if several levels are used the representative sampling error 
will be very small. Neglecting the sampling error, and assuming for sim- 
plicity that s is the same at each level, the standard error of any component 
of W becomes: 


r= n (gm/cm-sec) (11) 


After converting the divergence of water vapor over an area into a line in- 
tegral, and after assuming linear variation of W between stations on the 
perimeter: 


JL,  (gm/em*-sec) (12) 


where W., = W, and Ly is a unit vector normal to the line L, and directed 


outward from the polygon. Since the stations were equally spaced around the 
perimeter, Lj is a constant L. It will also be assumed that the standard er- 
ror of estimate of each component of W does not vary from station to station. 
The standard error of estimate of V-W is then given by: 


‘he (gm/em?-sec) (13) 


Sampling error has again been neglected: i.e., the assumption of.linear vari- 
ation of W between stations has been assumed valid. 
Introducing the perimeter P', and the number of levels n from equation 
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Fig. 9. An idealized equilateral polygon used in evaluating errors 
in water vapor divergence. Vertical observations are for 


equal increments of pressure. 
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(11), the standard error in V-W becomes: 


S.6.= SHE ls (gm/em?-sec) (14) 
n 


where H is the depth of the moisture bearing atmosphere, The 0.95 confi- 
dence interval for V-W is given by: 


A nn 


Although this formula was computed for an ideal case, and certainly not 
for any of the polygons treated in the text, valuable information may, never- 
theless, be obtained from it. Thus it can be seen that the accuracy of the 
divergence term may be increased by employing a greater number of stations, 
levels, and observations. It is also evident that the size and shape of the area 
also play an important role in obtaining a true estimate of the divergence. An 
optimum region would be a large area whose shape approximates a circle. 
Let us now consider the problem of obtaining a numerical value of the 
confidence interval from expression (15). The only unknown term is s, the 
standard error of estimate of the components of w (the mean transfer of 
water vapor at any station-level). It is obvious that s includes three types of 
error: representative sampling error, random instrumentation error, and 
consistent error. The last type of error is comparatively unimportant, since 
its percentage effect will not increase upon differentiation. Thus we may 
write: 


$ (gm/cm-mb-sec) (16) 


where sj is the contribution due to representative sampling error, and so is 
the contribution due to random instrumentation error. 

By estimating the average standard deviation of wind speed and mixing 
ratio within a twelve hour period (the interval between successive observa- 
tions) it was found that a reasonable figure for sj was 0.7/ /N gm/cm-sec. 
Also, if the observational error in measuring wind speed is about 1.5 m/sec 
and in measuring mixing ratio is about 0.3 gm/kg, sg will be approximately 
0.8//'N gm/cm-mb-sec. Thus the value of s is about 1.1///N. 

Subject to the stringent assumptions inherent in the derivation, the confi- 
dence interval for V-W thus becomes:5 


where N is the number of observations at each station-level. 


5. For example, it must be noted that expression (17) would give a confidence 
interval equal to zero if an infinite number of stations, levels, or observa- 
tions were employed. It is clear, therefore, that expression (17) is only 

approximate as a result of the idealization of the problem and must be in- 
terpreted carefully. 
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The Accuracy of Evapotranspiration Estimates 


In the present section, the confidence interval for divergence of water 
vapor has been computed for four sub-continental regions and for North 
America using equation (17). The data are presented in Table 3. E is evi- 
dent that the P/A ratio is the dominant factor. As an approximate rule, the 
accuracy of the method varies inversely as the square root of the area. 

Table 4 compares the values of the annual evapotranspiration obtained by 
the proposed method with those obtained from other methods. In four out of 
eight comparisons the confidence interval does not include the estimate de- 
rived by use of other methods. Although this is not conclusive, since the 
other methods are themselves subject to considerable question, the possibili- 
ty must be considered that the confidence intervals are too low. This is 
likely, since the assumptions made in deriving equation (17) were not always 
realistic. 

In particular, two assumptions are open to question. First, it was as- 
sumed that there were no sampling errors introduced due to interpolations 
between levels or in extrapolation to the surface. However, water vapor 
transfers were computed only at 500, 700, and 850 mb, and as has been noted 
earlier the extrapolation to lower levels is very critical, since maximum 
transfers sometimes occur below the 850 mb level. Although every effort 
has been made to achieve accurate extrapolations, there is no assurance that 
the results are completely satisfactory. 

A second assumption which may often be invalid is that the variation of 
water vapor transfer from one station to another can always be adequately 
described by a linear function. Assuming the average spacing of stations to 
be more or less constant, errors due to non-linearity would tend to increase 
as the size of area (and hence the number of stations) is decreased. It is en- 
tirely possible that errors of this type also appear in Table 4. 

It is of interest to note that inaccuracies in the evaluation of the atmos- 
pheric balance equation which arise due to lack of adequate representative 
sampling in the vertical and around the periphery of the area can largely be 
eliminated by use of more levels and a larger number of radiosonde stations. 
The improvement in evapotranspiration estimates to be obtained by the use 
of more extensive meteorological data has yet to be determined by actual 
computation. However, it is entirely possible that the success achieved for 
the North American continent can be extended to areas as small as the Ohio 
River Basin and to time intervals as short as ten days. 


CONCLUSIONS 


This study has presented a method of evaluating evapotranspiration based 
on sound physical principles. The method has been applied to various sub- 
regions of the North American continent; and the estimates obtained have 
been compared with hydrologic data where available, and with Thornthwaite’s 
empirical method. The proposed method yields reasonable monthly and annual 
totals of evapotranspiration when very large areas of almost continental size 
are considered. However, the technique may become unsatisfactory when 
substantially smaller areas (e.g., of the size of the Ohio River Basin) are 
treated. 

One obvious possibility is that this decreased accuracy for smaller areas 
is related to statistical sampling, and the fact that the divergence of water 
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vapor represents a small difference between two large numbers (the total in- 
flow and outflow). Statistical analysis of an idealized system leads toa 
formula fcr the confidence interval for evapotranspiration which verifies that 
the error should increase as the area is decreased. However, numerical 
evaluation indicates that for smaller areas this estimate of the confidence 
interval is probably too small. This may, of course, be due to inadequacies 
in the expression for the confidence interval. For example, it was assumed 
that errors in extropolating moisture transfer from 850 mb to the surface 
and in assuming linear variation of transfer between stations were negligible. 
Despite the difficulties encountered, the potentialities of this method can- 
not be discounted. Improved accuracy can surely be obtained by use of more 
meteorological stations and more levels of observations. The use of ob- 
served winds for smaller areas is also a possibility, although the adequacy 
of moisture divergence fields obtained from such transfers is open to ques- 
tion. In any case, the atmospheric balance method has yielded monthly esti- 
mates of evapotranspiration for large areas which are at least as accurate 
as those obtained from any other known method. With further refinements 
of meteorological observations and techniques of data processing, this 
method of measuring evapotranspiration may prove to be of considerable 
practical interest to the engineer. 
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SYNOPSIS 


Chemical coagulation is a basic water treatment process, which prepares 
the water for subsequent treatment operations. Coagulation includes addition 
of chemicals, their distribution, floc formation and floc build-up. The con- 
ditioning process that leads to floc formation and build-up is “flocculation.” 
Flocculation is usually accomplished by prolonged agitation at low velocity. 

Commonly used coagulants are alum and ferric sulfate. Alum’s ionization, 
hydrolysis, and reaction with alkalinity are of interest but the process of 
coagulation cannot be explained on the basis of simple chemical relationships. 
Colloidal particles play an important role in coagulation. Some of the signifi- 
cant characteristics of colloidal dispersions are large specific surface area, 
high adsorptive capacity, and particle charge. Colloidal particles commonly 
encountered in natural waters carry a negative charge, whereas those formed 
during the initial phase of the coagulation process often have a plus charge. 
Ions possessing a charge opposite to those borne by the colloidal particles 
tend to precipitate the latter. An ion’s precipitating power increases tre- 
mendously with its valence. 

Important factors in the coagulation process are pH, dissolved solids, 
temperature, time-concentration, coagulant aids, exchange capacity, and 
mixing. The last named is essential for efficient coagulation and consists of 
two stages: rapid mixing for a short period followed by slow mixing for a 
much longer time. Flocculation is accomplished in the latter mixing stage. 

Laboratory studies of the coagulation process are essential. Flexibility 
is an important aspect of plant design. 


INTRODUCTION 
Clarity is one of several highly-regarded characteristics of water. This 


Note: Discussion open until January 1, 1957. Paper 1036 is part of the copyrighted 
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is especially true if the water is to be used for municipal purposes. Some 
industries also require process water virtually free of turbidity. The treat- 
ment of waste water—sewage and industrial wastes—involves concentration 
and separation of suspended matter, that is clarification of the waste water. 

One of the oldest methods of water treatment is the removal of suspended 
solids by subsidence,(1) However, it is a well known fact that the nature of 
the suspended matter in most natural waters is such that subsidence, un- 
assisted, will not, even in combination with rapid sand filtration, render the 
water suitable for general use. Chemical coagulation, a water treatment key- 
stone, is an almost universal requirement. Coagulation prepares the water 
for sedimentation and greatly increases the efficiency of the latter. Coagula- 
tion plus sedimentation prepare the water for filtration. Coagulation ranks 
high on the list of things that must be done to render water safe, clear and 
palatable. 

Much has been learned about the whole process of coagulation especially 
in the last 25 to 30 years; yet, today the practice of coagulation has certain 
of the empirical characteristics of an art. The process is complex; it is 
sensitive; and natural waters are endowed with rugged individuality. We are 
probably just beginning to obtain a fairly clear picture of a process that has 
been practiced for centuries. There is still room for exploration. 


Coagulation - Flocculation 


The coagulation process includes “addition of chemicals to the water, 
their uniform distribution through the water, and the building up of a readily 
settleable floc, usually by prolonged agitation”.(2) The process of condition- 
ing that leads to floc formation and build-up is termed “flocculation,” the 
subject of this discussion. Obviously it is difficult to separate “flocculation” 
from the overall process of “coagulation.” Any discussion of the former can- 
not help but include the latter. Therefore much of this discussion relates to 
the process of coagulation which includes flocculation. 

The aim of flocculation is the promotion of coalescence of the precipitates 
derived from the coagulating chemicals and the natural suspended matter into 
fast settling precipitates termed “flocs” in water treatment parlance. As will 
be mentioned subsequently in this paper “conditioning” or “promotion-of- 
coalescence” is commonly accomplished by slow mixing by hydraulic or, 
more frequently, mechanical means. An important step in the coagulation 
process, flocculation appears to merit the attention it has received especially 
in recent years. 


HISTORICAL NOTES 


Mankind because of his need for water and his curiosity regarding it has, 
through the ages, exhibited much interest in its qualities and the means by 
which they might be altered the better to serve his current needs. Conse- 
quently it is not surprising to note that coagulation of water was practiced in 
ancient times. Apparently there was appreciation of the need for mixing as 
early as the sixteenth century for an Italian physician who visited Egypt then 
described the use of almonds for coagulation as follows: 

“,.. they smeared the edges of the vessel with five sweet almonds, prop- 
erly crushed, and grasping the almonds in the hand suddenly plunged hand 
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and arm into the water up to the elbow and moved elbow and hand vigorously 
and violently this way and that through the water stirring it up until they had 
made it far more turbid than before. Then they withdrew the arm from the 
vessel leaving the almonds in the water. They let it clarify and it was prop- 
erly clear in three hours time.”(1) 

In the United States coagulation was used in connection with a New Jersey 
public water supply in 1885 and studies of aluminum sulfate as a coagulant 
were undertaken at Rutgers University during the same year.(1) Since 1885 
coagulation “has been, with few exceptions, part and parcel of rapid filtration 
in America.”(1) 


Some Chemical Aspects of Coagulation 


Chemicals(3)(17) 


Numerous chemicals have been considered for use as coagulants. How- 
ever, relatively few are in common use today. These are listed in the follow- 
ing table. 


TABLE I 
CHEMICAL COAGULANTS 


Name Approximate Formula 


Aluminum Sulfate Alg(SO4)3 15H20 
(Alum, Filter Alum) 


Ferric Sulfate Feo(S04)39H20 
(Ferrifloc, Ferrisul) 


Ferric Chloride or FeCl36H20 
Chlorinated Copperas FeCl and Fe2(SO4)3 


Ferrous Sulfate FeS0O47H20 
(Copperas) 


Lime 
Hydrated lime Ca(OH)2 
Quicklime CaO 


Activated Silica SiO, 


An important characteristic of aluminum and iron salts is their ability, 
under appropriate conditions, to form gelatinous-appearing precipitates or 
“flocs.” A more precise term would be “chemical flocs” to distinguish them 
from “biological flocs” produced in sewage treatment processes. Of the 
chemicals listed in Table I aluminum sulfate (alum) is probably the most 
widely used, It is cheap and generally available. Lime is frequently applied 
merely as a coagulant aid but in some waters it alone can produce effective 
coagulation. Lime is a widely used softening and coagulating agent for the 
hard waters encountered throughout much of the mid-western section of the 
United States. Consequently its inclusion in the list of commonly used coagu- 
lants is believed justified. Activated silica is a newcomer and is generally 
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considered a coagulant aid. However, at least one authority (3) has classed it 
as a coagulant. 


Chemical Reactions 


Although alum has certain inherent chemical peculiarities it is widely 
used and therefore alum has been selected to illustrate certain coagulation 
reactions. The commercial product is approximately Alg(SO4)3 15H20. 
(Detailed specifications have been prepared by the American Water Works 
Association).(4) When alum is dissolved in water it ionizes yielding alumi- 
num ions (Al+++) bearing a high positive charge and negatively charged 
(SO4~~) ions. This may be illustrated by the equation 


Aly (SO4)3 2Al+++ + 38047~ 
Alum also hydrolyzes readily as follows: 
Alo(SO4)3 + 2A1(OH)3 + 


In most natural waters the sulfuric acid formed by this reaction is prompt- . 
ly neutralized by the natural alkalinity of the water. (This is assumed to be 3 
calcium bicarbonate.) 


3Ca(HCO3)o + 3CaSO4 + + 6H2O 


(The acid may also be neutralized by the addition of an alkaline agent, com- 
monly lime.) 

The overall equation, assuming 15 molecules of water of crystallization, is 
often written as follows: 2 


Alp (SO4)3 15 HgO + 3Ca(HCO3)2 2A1(OH)3 + 3CaSOq + 6CO2 + 15H20 


Aluminum hydroxide (more accurately the hydrous oxide of aluminum) is 
an amphoteric substance. This means that under acid conditions (low pH) it 
behaves as a base and under alkaline conditions it behaves as an acid, Con- 
sequently it will be dissolved by too great a surplus of either H+ or OH™ ions. 


AL(OH)3 + 3H* — Al*++ + 
(Solution at low pH) Zz 
Al(OH)3 +OH” AlO,” + 


(Solution at high pH) 


Precipitates formed by iron coagulants do not exhibit amphoteric behavior, 
at least not the extent shown by those derived from alum, 

These reactions are intended to serve as simplified illustrations and do 
not purport to indicate the exact chemistry of the coagulation processes nor 
the composition of the alum floc.(5) The chemistry involved is much more 
complex than these simple equations suggest. However, certain observa- 
tions of a qualitative nature may be made based on a consideration of these 
equations. There are: 


a) Alum yields both plus and minus ions of fairly high charge. 

b) Alum produces a precipitate. (The composition varies). (5) 

c) Alum is an acid salt. Its hydrolysis can result in a reduction in the 
natural alkalinity of the water and an increase in the carbon dioxide content. 
The result is an increase in the hydrogen ion (H+) concentration or a reduc- 
tion in the pH value of the water. 


4 
~ 
| 
| 
. 
vat 
| 
of 


ASCE 1036-5 


d) Because of its amphoteric nature the solubility of alum floc is sensitive 
to the pH value of the water. pH adjustments may be required to achieve op- 
timum conditions for precipitation. 


Similar equations leading to generally similar conclusions could be written 
for most of the iron coagulants (ferric sulfate, ferric chloride, chlorinated 
copperas) Ferrous sulfate (copperas) behaves similarly; however, its reac- 
tion is further complicated by the necessity of oxidizing the iron from the 
“ferrous” to the less soluble “ferric” form. As previously noted an impor- 
tant difference between alum and the iron coagulants is the fact that the pre- 
cipitates or flocs formed by the latter are not amphoteric. 


Colloidal Particles 


The coagulation process cannot be described in terms of simple, stoichio- 
metric relationships. Coagulation is intimately associated with the science 
of colloidal materials. Toa great extent understanding of the mechanism of 
coagulation is dependent upon knowledge of the characteristics and behavior 
of matter in the colloidal state. Colloid chemistry has been characterized as 
“the twilight between chemistry and physics.” (6) 

Colloidal dispersions may be thought of as consisting of a homogeneous 
medium such as water and particles dispersed therein. The term “colloid” 
does not refer to a particular kind of matter but to a somewhat indefinite 
state of subdivision. In general colloidal particles are not microscopically 
resolvable nor can they be readily separated from water by filtration methods. 
The author is not aware of any satisfactory, short definition of colloidal sys- 
tems. No true dividing line exists between what may be termed suspensions, 
colloids, and solutions. Colloidal dispersions have been referred to as 
“colloidal suspensions” and as “colloidal solutions.” Actually there is a 
“continuous transition” from suspensions through colloids and on to so-called 
true solutions. (6) 


Particle Size 

Colloidal particles are often arbitrarily considered to be those falling in 
the approximate diameter range 1.0 to 100 millimicrons (10-7 to 10-5cm) al- 
though some authorities suggest a range of “over 1.0 millimicron to some- 
thing less than 500 millimicrons.”(6) Such distinctions are purely arbitrary 
and are based on expediency. 

For comparative purposes the following table of approximate particle 
sizes is presented, 


TABLE II 


PARTICLE DIAMETERS (7) 


Material Diameter 
(Order of Magnitude) 


Gravel 1 cm upward 

Coarse Sand cm (10 

Fine Sand Oyl.cm (10° “cm. ) 

Silt 107 j,cm. 

Bacteria 10° ‘cm. (1 micron) 
Clay 107 

Virus 10° gem. (1 millimicron) 
Oxygen molecule 10 “cm. 
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Surface Area 

Some of the distinctive properties of colloidal particles can be explained 
by their large specific surface area. Assume that water contains 200 parts 
per million of a material having a specific gravity of 2.5 and that this mate- 
rial is colloidally dispersed in the water as spherical particles having a di- 
ameter of 1.0 millimicrons (10-6cm). Each 1000 gallons of water would con- 
tain about 6 x 1020 particles and their surface area would be approximately 
45 acres. Obviously surface activity plays a prominent role in colloidal be- 
havior. 


Charge 

Colloidal particles normally bear a charge. This charge probably can be 
obtained in various ways, including ionization of a part of the colloidal par- 
ticle. The charge may also be obtained as a result of the colloid’s surface 
activity; i.e., by preferential adsorption of various ions from the water in- 
cluding hydrogen (H+) or hydroxyl (OH™) ions.(8) The stability of many col- 
loids can be explained by the fact that the particles carry like charges. — 
Similarly charged particles repel each other; thus, continued dispersion 
rather than precipitation is favored. The degree of stability of a particle is 
determined by a value termed the “zeta potential” which is related to par- 
ticle charge and the distance through which the charge is effective. (8) 

Particle charge can be demonstrated experimentally using colloidal ma- 
terial found in natural waters. For example Saville,(9) working with a 
colored water, showed that the color was caused by negatively charged col- 
loidal particles. Saville actually demonstrated electrophoresis, a character- 
istic property of colloidal particles having to do with their migration under 
the influence of an electrical potential. Colloidal particles encountered in 
natural waters such as clay and organic color are commonly negatively 
charged. On the other hand colloids formed in the initial stages of the coagu- 
lation by aluminum and iron salts generally carry a positive charge. 


Precipitation of Colloidal Particles 


Colloidal particles can be precipitated by ions of opposite charge. The 
relative precipitating power of ions increases tremendously with increase in 
the ion’s charge or valence, For example sulfate (SO4--) is more than 100 
times as effective as chloride (Cl~) in precipitating positively charged col- 
loidal particles. The aluminum ion (Al***) is roughly 500 times as effective 
as the sodium ion (Na*) for precipitating a negatively charged colloidal par- 
ticle. The greater the concentration of dissolved mineral solids (ionic 
strength) the more readily precipitation occurs. (An increase in ion concen-. 
tration lowers the zeta potential), Therefore it is generally easier to accom- 
plish coagulation in waters of high mineral content than in those containing 
low concentrations of mineral matter. 

Inasmuch as ions responsible for charge on colloidal particles are often 
hydroxyl (OH-) or hydrogen ions (H+) colloidal-particle precipitation should 
be influenced by the pH value of the water, Water treatment experience sup- 
ports this statement. The process of coagulation is often quite sensitive to 
pH changes. Recognition of the importance of the pH factor in coagulation 
may mean the difference between good and poor results. 
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Coagulation Process 


Mechanism of Coagulation ay 


Black(3) whose studies of coagulation phenomenon have contributed greatly 
to our knowledge of the subject has divided the process into three phases. 

a) The coagulant furnishes trivalent aluminum or ferric ions (high plus 
charge) capable of neutralizing the negative charge on colloidal particles of 
color and turbidity. 

b) If the coagulation reaction occurs below about pH 7.0 the result of the 
initial action of the coagulant is the formation of positively charged particles a 
of colloidal dimensions. Black has termed these “micro-flocs.” Their plus ; 
charge at low pH is probably due to adsorption of hydrogen ions and surplus 
aluminum or ferric ions. These positively charged particles have adsorptive 
capacity and neutralizing power and can continue to remove negatively 
charged particles of turbidity and color. 

c) The third or flocculation phase involves further adsorption, charge 
neutralization and finally agglomeration of the micro-flocs into larger and 
larger particles which will settle readily. 

It is apparent that two types of colloids are involved in coagulation, name- 
ly, those present in the water undergoing treatment and those formed by the 
coagulating chemicals added to the water. In general the former are nega- 
tively charged whereas the latter commonly carry a plus charge. However i 
reversal of the charge on particles formed by the initial coagulation reactions oI 
may occur at approximately pH 7.0 and above. This is probably due to ad- a 
sorption of hydroxyl ions. 


Factors Affecting the Process 


pH Value 

The pH value of the water affects the solubility of precipitates formed by 
iron and aluminum, time required for floc formation, and charge on colloidal Z 
particles. Effect of pH value on the time required for floc formation has been iy 
studied by Theriault and Clark;(10) Bartow, Black and Sansbury;(11) Black, 
Rice and Bartow;(12) Baylis(13) and numerous other investigators. Based on 
theoretical considerations and certain experimental evidence Theriault and 
Clark concluded that when other influences are not considered optimum con- 
ditions for alum-floc formation will be found at about pH 5.5. From a prac- 
tical standpoint other influences are important but it is often the case that 
alum coagulation is best accomplished at fairly low pH values (5.0 - 7.0). 
This is particularly true of soft, colored waters, 

The pH range of relative insolubility for aluminum is approximately 5.0 to 
7.0, above pH 4.0 for ferric iron and above pH 9.5 for ferrous iron. ( 


Dissolved Mineral Matter 

Dissolved mineral matter in the water is a source of ions which can play 
an important part in the coagulation process. Coagulation of the positive 
colloidal particles (micro-flocs) formed by iron or aluminum coagulants is 
influenced considerably by other ions. The importance of any particular ion 
is to a great extent related to its valence or charge. As previously men- 
tioned, the sulfate ion which bears a fairly-high negative charge strongly 
favors precipitation of positively charged micro-flocs—much more so than 
the chloride ion which carries a lower charge. Various investigators have 
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demonstrated that sulfate definitely extends the zone of optimum floc forma- 
tion in the lower pH range. The magnitude of this effect is somewhat depend- 
ent on sulfate concentration. Sulfate is a constituent of two commonly used 
coagulants, alum and ferric sulfate, and is found in a wide range of concen- 
trations in natural waizis. 

At pH values of approximately 7.0 and above the charge on the micro-flocs 
may be negative rather than positive. If such is the case then positively 
charged ions particularly calcium and magnesium (high plus charge) would 
play a part in floc formation. 

In the complex process of coagulation probably all ions have some influ- 
ence, It is therefore apparent that natural waters, which vary widely in 
chemical make-up, will respond in various ways to the addition of coagulating 
chemicals, 

Effects of synthetic detergents on coagulation have been investigated by 
Vaughn(25) and associates at the South District Filtration Plant, Chicago, and 
by Smith(26) and others at the Robert A, Taft Sanitary Engineering Center, 
Cincinnati, 


Temperature 

The principal effect of temperature is its influence on the rate of floc 
formation. Generally the lower the temperature the longer the time required 
to produce good flocculation. Higher coagulant dosages may be required for 
cold water than for warm. Powell(14) has presented curves illustrating in- 
creased chemical requirements at temperatures below 50° F. Camp’s 
studies(15) showed that decreasing the temperature shifted the optimum pH 
to slightly higher values. Thuma(16) has reported recently on the adverse 
effects of low water temperature on water softening. 


Time-Concentration 

Various authorities(8)(17) have called attention to the fact that the coagula- 
tion process is a time-concentration phenomenon. High concentrations of 
coagulant produce results in shorter time than low concentrations. 


Coagulant Aids 

Discovery that partially neutralized sodium silicate solutions would sub- 
stantially improve coagulation was a notable advance in water treatment.(13) 
The silica activation process yields a negatively charged colloidal particle 
which not only hastens the rate of coagulation but also results in the forma- 
tion of larger floc particles. Activated silica is also beneficial in connection 
with water softening by the lime-soda ash method and has been successfully 
used as a substitute for the conventional coagulants, (18 

Lime, soda ash and acids (sulfuric and phosphoric) are coagulant aids in 
the sense that they can be used when necessary for pH adjustment, 

Suspended particles can assist flocculation by acting as nuclei for floc 
formation, Coagulation in some instances may be improved by the addition 
of small amounts of finely divided, insoluble material. Clay, bentonite, and 
activated carbon have been used for this purpose.(17) Nuclei may.also be 
provided by returning sludge from previously treated water. 

Application of chlorine ahead of the coagulation process has been reported 
to be advantageous.(14)(19) Chlorine may alter or destroy materials which 
inhibit floc formation and render other materials more susceptible to coagu- 


lation. 
Polyelectrolytes are a recent development of considerable interest. These 
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are synthetic organic compounds and are said to be useful in flocculating a 
wide variety of suspensions. One application has been in the treatment of 
suspensions encountered in ore processing. Apparently they can be used 
alone or in combination with conventional coagulants. 


Exchange Capacity of Turbidity Particles 

Investigations by Langelier and associates(20)(21) have indicated that an 
important factor in the coagulation process is the cation exchange capacity 
of the turbidity particles. When exchange capacity is fairly low, flocculation 
requires a binder material such as alum or iron floc. An amount of coagu- 
lant proportional to the exchange capacity must be added to provide colloid- 
destabilizing cations (Al+++ or Fe+++) and an additional amount is required 
to furnish binder material. Formation of the latter by coagulant hydrolysis 
occurs only after the turbidity particles have been destabilized by the coagu- 
lant cations. 


Mixing 

The importance of appropriate mixing cannot be over-emphasized. With- 
out it efficient coagulation is virtually impossible. It is important to recog- 
nize and distinguish between the two types of mixing commonly employed in 
coagulation. These are “rapid” or “flash” mixing and “slow” mixing or floc- 
culation. Rapid mixing requires only a few seconds using a high-speed mix- 
ing device, Flocculation on the other hand requires a much longer period of 
mixing at fairly low velocity. 

The first and probably one of the most important phases of coagulation in- 
volves the action of aluminum or ferric ions provided by the coagulant. It is 
therefore important that rapid dispersion of the coagulant through the water 
be achieved. The function of a rapid mixing device is that of quick and 
thorough dispersion of the coagulant. 

The initial phase of coagulation produces finely divided particles and op- 
portunity for coalescence must be provided if a readily settleable floc is to 
be obtained. Coalescence requires contact between particles and there is 
greatly increased opportunity for such contact if the water is slowly mixed, 
The result of slow-mixing is floc “build-up” or “flocculation.” (The process 
is often referred to as “conditioning.”) 

Flocculating velocities employed vary but often range from about 0.9 to 
1.5 feet per second.(22) Too high a velocity will prevent the formation of 
large floc particles or break up those that do form. Flexibility in mixing 
speed is desirable in most instances. Time allowed for slow mixing varies 
but an allowance of 30 to 45 minutes is common practice. Camp(23) has re- 
lated the rate of flocculation to velocity gradient and detention time. 

Following the flocculation phase of the coagulation process turbulence 
must be avoided, The floc particles are fragile and can be broken up in pas- 
sage from the flocculation chamber to the settling basin. The transition from 
flocculation to sedimentation should be as direct as possible and be accom- 
plished at low velocity. If conduits are required velocity in them should in 
general be less than 1.0 foot per second. 


Coagulation Control 


In view of the complexity of the coagulation process and the well-known 
variability of natural waters, it is not surprising that coagulation poses a 
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difficult control problem for the water-plant operator. For practical operat- 
ing control there is no substitute for trial determinations of coagulant dosage 
conducted with laboratory stirring equipment.(24) Optimum dosage, pH ef- 
fects, and use of coagulation aids can be readily studied in the laboratory and 
it is essential that this be done. The results of laboratory studies, to be use- 
ful, have to be interpreted in terms of actual plant performance. For this, 
there is no substitute for experience, In this connection one further point 
needs emphasis. Laboratory data indicating the need for adjustments are of 
little value to an operator if the means of making the required changes are 
not available to him. Obviously there is often need for the provision of flexi- 
bility in the way of chemicals and facilities for their application, points of 
chemical feed, and mixing. The coagulation process cannot be forced into a 
rigid pattern, Although coagulating chemicals have been used in large-scale 
water treatment for more than a half century, we still have much to learn 
about their behavior and the nature of the substances dispersed in natural 
waters. 
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A FLOW CONTROLLER FOR OPEN OR CLOSED CONDUITS 
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(Proc. Paper 1037) 


SYNOPSIS 


A general principle of flow control, utilizing a non-linear resistance, is 
developed that is applicable to both open and closed conduits. The control is 
infinitely adjustable for discharge over its design range, and holds the flow 


close to the predetermined discharge over its head range. The control is 3 

readily converted into a flow meter with a linear head-discharge curve and a 

adjustable range and sensitivity. Experimental work is presented to confirm a 


the theory. 


INTRODUCTION 


The combination of a disc moving within a profiled throat, and a nonlinear 
resistance to support the disc against the pressure drop is shown to yield a 
principle of flow control that is infinitely variable and theoretically exact (i.e., 
with no change in discharge or drag coefficients). By reversing the throat 
section and measuring the head across the disc, the device is converted into a 
flow meter. The underlying principles of single orifice flow control are out- 
lined briefly before developing the specific relationships and design data. A 
model has been constructed and tested in the Civil Engineering hydraulic 
laboratory of the University of Michigan. The results are presented for both 
the flow control and the flow meter. 


Single Orifice Flow Control Concepts 


An orifice flow control must act to hold the discharge, Q, constant in 
Q = Cy av2g h (1) 
by changing the orifice area, A, as the head, h, varies over its design range. 


Note: Discussion open until January 1, 1957. Paper 1037 is part of the copyrighted 
Journal of the Hydraulics Division of the American Society of Civil Engineers, Vol. 
82, No. HY 4, August, 1956. 


* Prof. of Hydraulics, Civ. Eng. Dept., Univ. of Michigan, Ann Arbor, Mich. 
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Cq is the discharge coefficient, and g the acceleration due to gravity. Assum- 
ing the discharge coefficient constant, the area must vary inversely as the 
square root of the head. For a control meant to operate at a single discharge 
only, this may be accomplished by means of a linear spring assembly, Fig. 1. 
Head across the device displaces the disc, thereby altering the net area. The 
throat may be profiled so that the discharge is held constant. 

If one attempts to make a variable flow control from this device, by varying 
the position, z, of the spring base BB, Fig. 1, the dimensionless head-dis- 
charge curves appear as in Fig. 2, with Qq the design discharge, hp the mimi- 
mum design head, and 7) the ratio of maximum design head to minimum design 
head. 

By selecting another approach, although still retaining a linear spring, the 
throat profile may be constructed so that the discharge is the same at mini- 
mum and maximum heads, for any setting, z, Fig. 1, within its size range. 
The discharge at intermediate heads varies from that at the end points, how- 
ever, as indicated in Fig. 3. To obtain the area required, as a function of x' 
(Fig. 1), let Y be the displacement of the disc for head change from ho to hoM, 
i.e., from minimum to maximum design head. The area must change from 


A, to Ag as given by 
Bo = AgV (2) 


which must hold for any setting, z, of the spring base BB. As the area of ori- 
fice opening is a function of x only, A = f(x"), 


f(x') = f(x - Y) (3) 


Since this relation must hold for any initial x', a law of variation of A with x' 
of the following form is indicated 


A = Age ~* (4) 
where Ao, &, are undetermined constants. Substituting into Eq. (3) 


Y) 
= Age a(x ron 


or, simplifying 


a = 
2Y 


where ln is the natural logarithm. 
The variation of A with x' then becomes 


x'1n 
A = Ace (6) 


with Ap the opening between disc and throat when x = 0, i.e., A, = Amin- 

To improve the accuracy of this control for the same head and discharge 
range, two linear springs may be introduced into the assembly, as in Fig. 4, 
altering the first spring so that the disc is displaced $ by the first spring only 
as the head varies from hy to hg¥# and is displaced the remaining distance af 
as both springs act with the head changing from ho v7 to hon. This results in 
a dimensionless head-discharge curve, Fig. 5, in which the control is accurate 
for three heads and the maximum error has been reduced. 

By taking n linear springs, the first acting over the distance Y/n, the first 
three over 2Y¥/n to 3Y/n, etc. with the corresponding head ranges 
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hy to hon!/2, honi/n ho?2/” to hon 3/2, ete. 


the discharge is maintained the same at n + 1 heads with the maximum varia- 
tion in discharge decreasing as n increases. By increasing n indefinitely, a 


theoretically exact flow controller is obtained. The spring system then be- 
comes a non-linear spring. The relationships for this case are developed in 


the following section. 


Basic Equations for Nonlinear Resistance Flow Control 


The basic equations for the flow controller have been derived by the author 
in a different manner in a previous paper.! One form of the nonlinear resist- 
ance flow control is shown in Fig. 6. The spring element is of the restrained- 
tip cantilever type, in which the two leaf springs are forced together by a pro- 
filed backing that shortens their effective length. The spring element is a 
distance y = Y apart at minimum design head and contacts at y = 0 for maxi- 
mum design head. By rotating the milled knob, the discharge setting, z, can 
be varied. From the geometry of Fig. 6 it is evident that 

x'=y+z (7) 
where x' gives the disc position (the area of opening between disc and throat 
is a function of x' only); and y is a function of the pressure drop across the 
disc only, i.e., it is a function of h only, for a given fluid. 
For any setting z 


Avh = (A+ AA) vh +4h (8) 
or, since A = f(x") 
f(x’) = f(x" + ay) (9) 


because a change in Ah causes a change in Ay and Ax' = Ay for constant z. 
Expressing the area of opening, A, as an exponential function of x’, Eq. (4) 


f(x') = A= Amin e** (10) 
where & is as yet undetermined. Substituting Eq. (10) into Eq. (9), and solving 


for a, 


Expanding the logarithmic term into a series and taking the limit as Ay ap- 
proaches zero, 


-A _ 1 ah* 
2 2htay * 


1 @dh 
2h dy 


Since @ is a function of x' only, and y is a function of h only, Eq. (11) shows a 
function of x' equal to a function of y'. This is impossible because of Eq. (7), 


(11) 


1. “A Quick Response Variable Flow Control Device,” V. L. Streeter, Instru- 
ment Society of America, Vol. 2, No. 2, Feb. 1955, pp. 48-51. 
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hence both sides of Eq. (11) must be constant. Integrating Eq. (11) 
-24y = 


Using the end conditions: h = ho, y = Y; h = hon, y = 0, botha and c are de- 
termined, 


(12) 


and 


(13) 
Equation (10) becomes 


x in”? 
A= e (14) 


Now, substituting into Eq. (1), using h from Eq. (13), A from Eq. (14) and x' 
from Eq. (7) 


z in 
Q= Cg 26 hoy (15) 


which yields the discharge for the given arbitrary setting z, and shows that Q 
is independent of head h. 

The three equations, Eqs. (13), (14), and (15) determine the characteristics 
of the flow control. The nonlinear resistance must follow the head-deflection 
relation of Eq. (13), and the throat profile must be made so that the net area 
of opening is given by Eq. (14). By providing for the base of the resistance to 
be adjusted axially, the quantity z may be changed and variable flow control 
results, as given by Eq. (15). 


Equation for Throat Profile 


If an axially symmetric throat section is desired, the throat profile, as 
given by the cylindrical coordinates (r, x) in Fig. 7, must yield the net area 
for any x' required by Eq. (14). Using Pappus theorem this area is 


x'1n 
Amin e (16) 


——"¥ 
(r ro )e (ro+ r) V(r - ro)” + (x' - x)? (17) 


For a given x' this equation yields a curve of r as a function of x. Drawing 
segments of these curves for various values of x' , Fig. 8, shows that their 
envelope yields the desired throat profile. The partial differential of Eq. (17) 
with respect to x' is taken and solved simultaneously with Eq. (17) to yield the 
cylindrical coordinates of the throat profile. Dividing the partial of Eq. (17) 
with respect to x' into Eq. (17) yields, 


X'-xX = | (18) 
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Substituting Eq. (18) into Eq. (17) gives the desired equation, 


This equation permits the direct computation of x for any r in terms of the 
design quantities Y, 7, r', and ro. 

There is an inherent limitation to the value of r in Eq. (19), obtained by 
setting the radical term equal to zero, 


Y 
Tmax =To + Inq (20) 
For larger values of r the envelope solution breaks down. For r = ryyax, from 
Eq. (19) 


Xmax = Ra v2 (21) 
The corresponding value of x', from Eqs. (18) and (21), is 


Discharge Limitations 


For the flow control to operate the disc must stay within the range 
© <X'<Xtyax- For control over the whole head range this restricts z to the 
range 0 < Z = Xmax - Y. To find the maximum discharge over the whole head 
range, Z = X'max - Y is substituted into Eq. (15) using Eq. (22), 
Y Y 

Qmax = 27 VE ho my (2r,g + ) (23) 
making use of the fact that Amin = n(r'2 - r,2). The minimum discharge is 
obtained from Eq. (15) by setting z = 0, 


Qmin = Sq Amin ¥28 (24) 


Equations (23) and (24) permit the selection of the three design quantities ro, 
r', and Y for given head and discharge range. It is of interest to note that 
the maximum and minimum discharges are not related. The minimum throat 
radius r' does not come into Eq. (23). 

For flow control over partial head ranges, the limitation is that z varies 
from -Y to o in the lower extended zone and from xjmax - Y tO x max in the 
upper extended zone, with the disc (x') remaining between 0 and x max at all 
times. The spring travel is limited to portions of the range Y to 0, which re- 
stricts the head range to part of ho to ho). 

For z less than zero, the “threshold” value of y, to maintain x' = 0, is 


y=-2 
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Substituting this value of z into Eq. (15), using Eq. (13) for y, 


This is the lower limit line of the lower extended zone. 
For z greater than x{nax - Y, and with the disc at its limiting position 


x'= X'max, 
Z=Xinax 
Again, using Eqs. (15) and (13), the upper extended zone is limited by 
X max 127 
Q = Cq AminV2gh e” 2Y (26) 


The total flow control range is shown in Fig. 9 in dimensionless form, with 
the z range indicated on the right side. 


Design of Nonlinear Resistance 


The nonlinear resistance to displacement of the disc, as required by Eq. 
(13), may be expressed in terms of force exerted on the disc. This force, F, 
is 


F = wh (27) 
where Cg is the dimensionless drag coefficient, and w is the unit weight of 
liquid. The minimum design force, Fo, is then 

Fo = Caro" who (28) 
Cg is assumed to be constant over the range of travel of the disc. Substituting 
for h and ho in Eq. (13) 


(29) 


The resistance may be obtained in many ways, three of which are dis- 
cussed here, namely 1) fluid resistance, 2) dead weight resistance, and 
3) spring resistance. 


Fluid Resistance. Referring to Fig. 10, a bellows with area, Ap forces 
manometer liquid into a manometer leg of varying cross-section. The area, 
a, of manometer is determined as a function of s. The necessary equations 
are Eq. (29), 


-Ap dy = ayds; =a ds (30) 


and 


F = Ab (Wo - w) (s + 8) (31) 


where a, is the constant cross-sectional area and Wg is the unit weight of 
manometer fluid. F and y may be eliminated from the last three equations to 
yield a and s in terms of s; 


(32) 
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and 


where Spo is the value of s for sy = 0 and F = Fy. The discharge setting is 
altered by changing z. For vertical installations the weight of moving parts 
would change the setting z for a given discharge, but would not otherwise dis- 
turb proper functioning of the controller. 


Dead Weight Resistance. A weight suspended on a cable or metal tape such 

that it is lifted out by a cam can easily be arranged to yield the nonlinear 
force-displacement relation required by Eq. (29). Figure 11 illustrates one 
method, with the cam rotating about a transverse axis. To determine the 

shape of the cam, reference is made to Fig. 12, where the cam rotates through | 
the angle, «, from o when F = Fo, toa, when F = F,7. The moment of the 
weight, W, about the shaft for any angle o& is 


In 
Wp cos(@-a) = Fro = (34) 
With constant, 
Bae 


since ppa= Y- y. Letting W = Fo, Eq. (34) reduces to 


pcos (O@-X) = p,& % (35) 


For any &, Eq. (35) is a straight line in polar coordinates (p,@). The envelope 
of the family of straight lines is the cam profile equation. Taking the partial 
differential of Eq. (35) with respect to & and solving simultaneously with 

Eq. (35) to eliminate a yields 


tan’ 
r= e (36) 


which is a logarithmic spiral. By dividing the partial differential of Eq. (35) 
with respect to a by Eq. (35) 


tan a) = (37) 


showing that @ -c is a constant. @ varies from tan~! to + tan=! 
(Inn in Eq. (36). 

If a cam is to be designed to support the submerged weight, G, of the mov- 
ing parts of the controller, 


a inn 
Fo po Wp cos (@ &) (38) 
Taking the partial differential with respect to « and dividing by Eq. (38) 
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For a series of values of & between o and &,, corresponding values of 6 are 
found. Using these values the corresponding values of p are found in Eq. (38), 
yielding polar coordinates of the cam profile. 


Spring Resistance. A restrained-tip cantilever spring may be designed to 
provide the required resistance law by following a method due to Mr. S. P 


Clurman.2 Two identical, flat springs of constant width and thickness are 
placed as shown in Fig. 13. As the force, P, is increased the springs are 
forced against the profiled backings, thereby reducing their effective lengths 
and increasing their stiffness. The springs are fastened together at each end 
by a spacer block that prevents rotation. As there are actually four springs, 
each one takes half the load and is deflected through half the total displace- 
ment. The equation of the backing profile is determined as follows: 

Letting G = KF, be the submerged weight of moving parts of the controller 
for vertical stem applications, (G = 0 for horizontal applications) 


-y 1 
P= F-G=F, (qe - 


The sign of G has been taken for the case where the weight of moving parts 
acts in the opposite direction to the pressure force on the disc. With 6 the 
deflection of one spring from the position for minimum design head 


where P' is the load on a single spring (P/2). For the minimum design load 
the spring is cantilevered from its full length, 2, with the backing touching at 
Xo = 0, Fig. 14. The variables xp, yp are coordinates of the backing profile. 
With no load on the spring it will be displaced from the 6 = 0 position a 
distance 3 
L 
where E is the modulus of elasticity and I the section moment of inertia about 
the neutral axis. 
With (xo, yo) the last point of contact between spring and backing for a 
given active load P"', the deflection from no-load position is 


5 (t-x0)? | 2 


The right-hand side of the equation is obtained by consideration of deflection 
of cantilever beams, as follows: In Fig. 15, by superposition of the two simple 
loadings 


aYo 
(1-X9) ax, + Yo (41) 


dyo _ P'(L-xo)” M(2-x0) 
(42) 


Defl. = P (4-x0)? M(L-xo)? 
3 2 EI 


Replacing M by EI d?y9/dxo? and eliminating P’ in the two equations 


2. “The Design of Nonlinear Leaf Springs”, by S. P. Clurman, Transactions, 
American Society of Mechanical Engineering, Feb. 1951, pp. 155-161. 
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which, with the addition of yo, is the right-hand side of Eq. (41). 
6 and F, may be expressed in terms of the stiffness of the spring. From 


Eq. (40) 
261n7 
P = = (e - k) 


2¢é1n 
dP Folnq 


which is the spring stiffness, For § = 0 


dp’ Folny (44) 
ae 


é=0 
The stiffness of a spring varies inversely as the cube of its length, hence 


ap’ 
m4 1-%,)° 


Solving Eq. (43) for 6, using the last two equations, 


2 Foing dé 2 in7 ° 


From Eq. (44) 
Fo(1-k) Y 
24 EI 2 Inq 


Substituting the last two equations in Eq. (41) 


46 
21ny 2 in7 + 46) 


Xo 
Inserting the dimensionless quantities 


x - 


1n(1-x) = (1-x) + 4(1-x) + 6y 
In} dx dx 


This equation is converted into a linear equation with constant coefficients by 
the substitution of 


t = In(l-x) 
Then 


2- 
dy ay ze 
ag 
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and Eq. (47) becomes 


ay y= -k- 
¢ 6y ing (1-k-3¢) (48) 


The solution of this equation for the conditions 


¥=0,f=0, 0 


ay ing= -k - - + k) 


Substituting back for x, 
2y Ing = -k - 3 - 3 in(1-%) + (4g) 


This is the dimensionless equation for spring backing profile. 
From Eq. (45) the maximum value of X¥ = x,/ is, for § = ¥/2, 


(50) 


The cross-section of the spring must satisfy Eq. (44). 
Consideration of stresses in the spring indicates the maximum moment at 
the spring extremity. From Fig. 16 the moment Mj is given by 


M) = P L(1-%) - Mg (51) 
Mg may be found from Eq. (42), by obtaining dy,/dx, from Eq. (49) 


1-x 


Y | 3(142k) (1-%) + 


22 


Y 3 
2 | 1-x 


) 


(52) 


and from Eq. (44) 


Substituting the last two equations into Eq. (52) and simplifying 
1 
Replacing Mpg in Eq. (51) by this value, 
Fol 3 “a 
M) = 1 2k(1 »| 


M, is a maximum for x = 1 -7 -1/3, 


4 
4 
= 
aft 
= 
Xmax = 1 
a 
iq 
dyo _ Y dy 2 
dXxo dx 
and 
From Eqs. (40) and (45) 
Fof? 
= 
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(53) 


The maximum fibre stress, S,,,, is 


Sa abe 

where t is the spring thickness and b the spring width. Eqs. ie and (54) per- 

mit the thickness and width to be expressed in terms of S,,, 4, and the design 
quantities n, Y, k, and F,, 


By (37275 ~ 1 - 79) (55) 


7 


(39/3 - 1 - (54) 


F, (EY)? 


The restrained-tip cantilever spring is relatively simple to construct and 
has very small hysteresis loss. Its main objection, from the standpoint of a 
closed conduit controller, is its size. 


Flow Meter with Adjustable Sensitivity and Range : 
‘ If the throat casting in the flow control unit is reversed, as in Fig. 17, and 
a manometer is connected to the two ends of the throat, the pressure drop 
across the disc, as given by R the gage difference, is a linear function of the 
1 discharge. The relation among the variables x’, y, z' is 
x'+y=2z' (57) 
where z' + z is a constant. Substituting into Eq. (1), using Eqs. (13), (14) and iz 
; (57) to eliminate x', y, and A, . 
Q = Cq Amin [26 e 2Y h (58) 7 
Ro? 
showing Q to be a linear function of h. 
The sensitivity and range is changed by adjusting z, or z'. Equation (58) 
may be arranged into the dimensionless form 
‘In 
CqAminv2ghoy 7 he 
which is plotted in Fig. 18. 
The upper threshold of the metering range is given by a at 


may In 


(60) 


The lower threshold of the metering range is 


ae 
| 
and 
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Q h 
CqAmin V2 Roy Roy (61) 


Large settings, z, result in large gage differences for a relatively small dis- 
charge range, while small settings permit metering over a much wider flow 
range. The sensitivity of the meter, as given by Q/h, is an exponential func- 
tion of the setting z. 


Experimental Program 


An experimental flow control unit was designed to operate in a concrete 
flume in the floor of the Civil Engineering Hydraulic Laboratory. The unit 
was placed into a circular opening in the horizontal shelf of a barrier, Figs. 
19 and 20. A 45° V-notch weir was placed a short distance downstream from 
the controller. Hook gages were mounted upstream and downstream from the 
barrier, so that the head across the controller, as well as the discharge, 
could be determined. The disc has a displacement of 0.75 inches as the head 
varies from a minimum of 2 inches to a maximum of 8 inches. The disc di- 
ameter is 7.50 inches and the minimum throat diameter was designed to be 
7.60 inches. The minimum throat section was later measured to be 7.605 
inches. 

With the design quantities Y = 0.75 inches, h, = 2 inches, n = 4, rg = 3.75 
inches, and r' = 3.80 inches, cylindrical coordinates of the throat profile were 
computed from Eq. (19). The restrained tip cantilever spring, Fig. 19, was 
used for the nonlinear resistance, and was computed from Eqs. (49), (55), and 
(56) for 9 = 3 inches and k = 0, i.e. not allowing for weight of the disc assem- 
bly. The submerged weight of disc assembly was balanced by use of styro- 
foam floats placed on a downstream extension of the shaft. The amount of 
buoyancy was adjusted so that y = 0.75 inches when h = 2 inches. Originally 
the drag coefficient for the disc was assumed to be unity. Based on this co- 
efficient the spring bottomed at h = 6.5 inches. By adding extra width to the 
spring the resistance was increased so that the disc displacement was 0.75 
inches for head change from 2 inches to 8 inches. 

The styrofoam floats also served another purpose by stabilizing the disc 
for the larger discharges. Without stabilization the disc tended to oscillate 
through its whole stroke for discharges greater than about 50 g.p.m. A steel 
can was anchored to the floor of the flume concentric with the disc and floats, 
with about 0.25 inches radial clearance between can and floats. The resist- 
ance to forcing water into and out of the can provided adequate stabilization 
for all the experiments. Limitations in flume height restricted the range of 
heads for the larger discharges. 

Fins were placed in the throat section to guide the disc. These were found 
to cause excessive friction, and were cut back radially to avoid contact with 
the disc. A teflon bearing in a spider at the small end of the throat provided 
excellent guidance with a minimum of resistance. 

Test results for the flow controller are shown in dimensionless form in 
Fig. 21, with a summary of results in Table I. 

The throat casting was inverted to convert the unit into a flow meter. Head 
across the disc was measured with hook gages and discharge by means of the 
weir. Test results are shown in dimensionless form in Fig. 22, with a sum- 
mary of results in Table II. 
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TABLE I 


SUMMARY OF RESULTS---FLOW CONTROLLER 


Number of 

turns on Maximum . 

adjusting Head Range Discharge Divergence 
screw inches g-Pp-.m. percent 


0 2.32-8.27 20.0 1.4 


3 2.17-8.03 24.5 1.6 


2.0 


2.22-8.12 30.0 


2.12-8.72 35.8 2.2 


2.34-8.06 44.8 1.3 


2.8 


2.20-8.88 54.5 


2.29-8.46 66.6 1.9 


2.14-7.91 80.8 2.4 


1.92-7.59 


99.4 2.8 


1.88-7.17 123.8 


TABLE II 


SUMMARY OF RESULTS---FLOW METER 


Number of 


turns on Qgpm Maximum 
adjusting n= Divergence 
screw Rinches percent 


0 4.229 2.8 


3 


4.2 


5.225 


3.3 


6.559 


8.096 3.2 


9.876 


3.2 
11.972 


14.773 


18.334 


22.728 


28.372 
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Discussion of Experimental Data 


As only one model was constructed, there has been no opportunity to take 
advantage of experimental findings to increase accuracy. If a change in dis- 
charge coefficient occurs due to change of disc position within the throat, it 
may be compensated by redesigning the throat so that CA varies with x'as 
prescribed by the theory. In the present model no iaabrumentehien was pro- 
vided to determine the exact position of the disc for a particular setting. The 
setting z could be altered either by rotating the disc on its shaft (4 threads to 
the inch) or by rotating the adjusting knob over the spring (13 threads to the 
inch). The turns referred to in Figs. 21 and 22 are for the adjusting knob 
above the spring. 


Flow Controller. For the flow controller the “0 turn” disc position at h = 8 
inches was about at the minimum throat section. The average value of dis- 
charge for this setting, 20.05 g.p.m. was taken as Qmin. As the exact position 
of the disc in the throat is not known for this setting the area is unknown and 
likewise the discharge coefficient. 

Separate measurements of discharge coefficient were made with the disc 
in known, fixed positions, and without the upper supporting framework for the 
spring assembly. These results indicated a coefficient of about 0.62 within 
the throat and about 0.63 with the disc slightly above the throat casting. 

In each series of observations for flow controller or flow meter measure- 
ments were taken about equally for increasing heads and for decreasing heads. 
Hysteresis effects appear to be very small, much less than the s-shaped trend 
of the data, which apparently is caused by the spring. The width of spring was 
varied in the model, as well as the amount of styrofoam, until the displace- 
ment Y was 0.75 inches at h = 2 inches, and was zero at h = 8 inches. The s- 
shaped trend in the data may be explained by the spring backing being incor- 
rect to the extent that the spring was too stiff for the head range 2 to 5 inches, 
thereby holding the disc higher in the throat and increasing the flow area 
slightly. Also the spring was not quite as stiff as required by the theory for 
the range 5 to 7.5 inches and permitted the disc to be lower than required by 
theory, thereby decreasing the discharge. 

Figure 21 is plotted with log Q/ Qmin 25 Ordinate so that the same relative 
percentage error scale occurs for all settings, and so that the average dis- 
charge lines would be equally spaced. A change of three turns represents a 
displacement Az of 3/13 inch. 

Log Q/Qmin is plotted against number of turns in Fig. 23, where Q is the 
average of the discharges for a given setting. The solid line is determined by 
least squares and yields the equation 


0.994 e9-06704t (62) 
Qmin 


Comparing the exponent with z Inn/2Y, z = 13.78t. As z = 13t by direct meas- 
urement, the discrepancy can best be explained as due to a change in discharge 
coefficient. Tests with a fixed disc yielded a change in discharge coefficient 
when the spring support stands were in place that conforms generally to the 
relationship z = 13.78t. 

The disc was stabilized by placing the steel can around the styrofoam 
floats for the settings 15 to 27 turns. 
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Flow Meter. The “0 turn” position for the flow meter was near the minimum 
throat section for h = 2 inches. The discharge Q,,jn was taken as 20.05 g.p.m. 
as in the case of the flow controller. The data as plotted in Fig. 22 generally 
confirms the relationship that discharge varies linearly with head for any 
setting. A slight s-shaped trend in the data, similar to that of the flow con- 
troller, but reversed, can be explained by the spring being too stiff for small 
heads and not stiff enough for the larger heads. 

The slopes m = Q/h (Q in g.p.m., h in inches) for the best radial lines 
through each set of observations have been computed by the method of least 
squares and are shown in Table I, with the maximum percentage divergence 
of any single observation. The taking of data for the flow meter was rather 
difficult due to the large volume of storage upstream from the meter. Water 
was admitted to the flume from a constant head tank and the level in the flume 
came to equilibrium very slowly. Measurements were usually taken before 
steady flow was established, which could introduce minor errors in the ob- 
servations. The flow meter is inherently stable and did not require any form 
of stabilization. 

It was necessary to shift the disc on its shaft between the runs 0 - 12 and 
15 - 2%. Unfortunately, it was not adjusted so that the change from 12 turns 
to 15 turns was actually 3/13 inch. This is evident in the plot of log m/ Mo 
against number of turns in Fig. 24, where mg is the value of Q/h for 0 turns 
(mo = 4.229). By constructing the best line through the points for 0 - 12 turns, 
and for 15 - 27 turns, the relations between t and z were found to be z = 12.97t 
and z = 12.87t respectively, which confirm the known relation z = 13t. This 
indicates a substantially constant coefficient of discharge over the whole 
range of tests on the flow meter. 


CONCLUSIONS 


Theory and design data for both the flow controller and the flow meter have 
been developed, with several forms of nonlinear resistance. An open channel 
model was constructed using a restrained-tip cantilever spring for nonlinear 
resistance, and data was taken for ten settings each as a flow controller and 
as a flow meter. The results generally confirm the theory and point the way 
to the construction of accurate and relatively simple control and metering 
equipment. 
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NOTATION 


A},A, = area of opening between disc and throat 
a,aj = manometer cross-sectional areas 
ap = area of bellows 
b = width of leaf spring 
Cq = drag coefficient 
Cq = discharge coefficient 
E = modulus of elasticity 
e = base of natural logarithms 
F = pressure force on disc 
F, = minimum design force on disc 
f = functional relation 
G = submerged weight of disc assembly 
g = acceleration due to gravity 
h = head drop across disc 
hp = minimum design head loss across disc 


I = moment of inertia of spring section about neutral axis 


K, = stiffness of spring of length i 
k = ratio G/F, 
\ = length of spring for load F,/2 
Q' = length of disc support shaft 


In = natural logarithm 
M = moment 
P = force on spring assembly 
P' = force on single spring element 
Q = discharge 
R = manometer gage difference 
r = radial coordinate of throat 
r' = minimum throat radius 
= disc radius 
S,, = maximum fibre stress 
S,8,,8; = dimensions of manometer 
t = thickness of spring 
W = weight 
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w = unit weight of fluid - 
xX = axial coordinate of throat profile 
x' = position of disc in throat 
Xp = coordinate of spring backing 


x = dimensionless length x/Q 
| ¥ = value of y for h =h, 


y = position of spring 
y = dimensionless quantity yo/Y 
Yo = coordinate of spring backing 


z = discharge setting 
a = constant, angle 
B = angular deflection 


5 = displacement of single spring element from minimum design head 
position 


¢ = dimensionless variable 
7 = ratio of maximum design head to minimum design head 


@ = polar coordinate 


P, Po = polar coordinate 
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Fig.! Linear Spring 
Flow Controller 
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Fig.2 Head-Discharge 
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i Profiled Spring Backing 


Spring Block 


a Fig.6 Nonlinear Spring Flow Controller 


Fig.7 Notation for Disc Fig.8 Envelope 
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Fig.l! Profiled Cam Fig.12 Notation for Cam 
Equation 


Fig.13 Restrained Tip Cantilever Spring 
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/ 
Position for Minimum Design Head 


Fig.14 Notation for Cantilever Spring 
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Fig. 15 Deflections of Cantilever Spring 


Fig.16 Free body Diagram of Spring 
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Flow Controller in Place 


Fig. 20. 
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Fig.2! Flow Controller Test Data 
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Fig.23 Discharge Related to Number of Turns 
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HYDRAULICS DIVISION 
Proceedings of the American Society of Civil Engineers 


SEVEN EXPLORATORY STUDIES IN HYDRAULICS 


Hunter Rouse,! M. ASCE 
(Proc. Paper 1038) 


FOREWORD 


Immediately following the Sixth Hdraulics Conference, held at Iowa City 
in June of 1955, the Iowa Institute of Hydraulic Research conducted an inten- 
sive three-week course in experimental techniques, in which twenty-one 
graduate engineers from various organizations participated. The purpose of 
the course was not only to acquaint the class with new types of instrumenta- 
tion, but also to instill in them an appreciation for the methods of research. 

During the first part of the course, eight of the elementary experiments in 
mechanics of fluids taught at Iowa to engineering students were carried out 
and reported upon by the class in teams of three. These experiments includ- 
ed curvilinear flow with cavitation in a closed conduit; curvilinear flow with a 
free surface; viscous flow in a pipe; turbulent flow in smooth and rough pipes; 
form effects at conduit transitions; diffusion of a submerged jet; fall velocity 
of spheres; and flow past a circular cylinder. Messrs. D. W. Appel, P. G. 
Hubbard, E. M. Laursen, and D. E. Metzler of the Institute staff acted as 
instructors during this introductory period. 

The second part of the course involved eight experiments utilizing ad- 
vanced methods of measurement: the Pitot cylinder for constructing a two- 
dimensional flow pattern; resistance wires, pressure cell, and Sanborn re- 
corder for wave studies; the electrolytic tank for pressure distribution at a 
three-dimensional inlet; an electrical analog computer for the cross-thrust 
on neighboring cylinders; the hot-wire anemometer for submerged-jet turbu- 
lence; stagnation-tube surveys for boundary-iayer shear on a plate; multiple- 
resistance gage for scour at a bridge pier; and a strain-gage ¢:namometer 
for lift and drag on a foil, The foregoing instructors were assisted in these 
experiments by Messrs. M. Hug, S. C. Ling, T. T. Siao, and A. Toch. 

In the third part of the course, which is of primary concern herein, the 
members of the class were encouraged to undertake short investigations of 
an original nature, to the end of determining new and publishable information 
through application of the foregoing techniques. The following seven topics 


Note: Discussion open until January 1, 1957. Paper 1038 is part of the copyrighted 
Journal of the Hydraulics Division of the American Society of Civil Engineers, Vol. 
82, No, HY 4, August, 1956. 

1. Director, lowa Inst. of Hydr. Research, State Univ. of Iowa, lowa City, 


Iowa. 


1038-1 


| 


Paper 1038 HY 4 
: 
2 
4 
| 


1038-2 HY 4 August, 1956 
were suggested: 


Development of the non-circulatory waterspout 

Penetration of a jet into a counterflow 

Surface profiles at a submerged overfall 

Coefficient of contraction for a submerged jet 

Maximum amplitude of oscillation of a cylinder due to its vortex wake 
Deflection of a jet by a normal wall 

Influence of slope and roughness on the free overfall 


Experimental facilities for the investigations were improvised from equip- 
ment at hand, and a member of the Institute staff was assigned for consulta- 
tion to each of the seven teams, In the short time available, each team per- 
formed the required experimental measurements and computations, and six 
of them submitted preliminary reports to serve as a record for reference 
purposes. The seventh team, under the leadership of Mr. J. C. L Dooge, was 
able to give somewhat more time to its project, and the report that he finally 
submitted was well beyond the preliminary stage. Such amplification of the 
other reports as proved necessary to prepare them for publication in the fol- 
lowing form was carried out subsequently by the advisers under the writer’s 
editorship, and the illustrations were brought to their present state of uni- 
formity by Mr. Metzler. Although none of the articles even approaches the 
final word upon the subject, it is believed that each will provide both useful 
information and an incentive to further research, 
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DEVELOPMENT OF THE NON-CIRCULATORY WATERSPOUT 


J. Davidian,? J. M, ASCE, and J. E. Glover3 
(D. W. Appel,* J. M. ASCE, Adviser) 


Statement of the Problem 


The r:ost familiar occurrences of waterspouts in nature are those accom- 
panying cyclonic storms over the oceans. A vortex is formed by the inter- 
action of air masses, into the low-pressure core of which water is drawn, 
sometimes to great heights. The circulatory motion of the vortex is evident- 
ly an eseential part of the phenomenon. Less well known is the non- 
circulatory waterspout, which may form in the low-pressure z@ne produced 
by radial flow toward an intake in either of two fluid layers of different den- 
sity. In plumbing systems the phenomenon is known as back-siphonage, 
water as well as air being drawn into the unsubmerged end of an open supply 
pipe that is momentarily under vacuum. The purpose of this study was to 
determine experimentally the rate of flow necessary to develop a waterspout 
just to the point at which the second as well as the first fluid is carried into 
the pipe. This represents the critical limit for safeguarding water-supply 
lines from contamination, and for the general separation of liquids of nearly 
the same density by decantation. 

The geometric variables in the problem are the diameter D of the suction 
pipe and the vertical distance y from its open end to the interface between 
the fluids, If one assumes for the purpose of a preliminary analysis that the 
effects of viscosity and surface tension are negligible, the only significant 
fluid properties are the density of the upper (or moving) fluid P and the dif- 
ference between the specific weights of the two fluids Ay. As the flow vari- 
able, the mean inlet velocity V at which a waterspout forms has been cho- 
sen. When these variables are combined by using the techniques of 
dimensional analysis, the two dimensionless parameters obtained are y/D 
and pv2/a y D. The first of these is a geometric variable, whereas the latter 
represents the ratio of inertial and gravitational forces and hence is a form 
of the Froude number. This analysis indicates that a single series of mea- 
surements of the critical velocity for the formation of a waterspout over the 
practical range of relative heights of the suction pipe should suffice to deter- 
mine the function 


¢ 


However, to check the validity of the assumptions made in the analysis and to 


2. Hydraulic Engineer, U. S. Geological Survey, Atlanta, Georgia. 

3. Hydraulic Engineer, Waterways Experiment Station, Vicksburg, Missis- 
sippi. 

4, Research Engineer, lowa Institute of Hydraulic Research, State University 

of Iowa, Iowa City. 
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demonstrate the generality of the results, two series of tests were conducted. 
In the first series a waterspout was formed in water under air, and in the 
second the waterspout was formed in saline water under fresh. 


The Experiments 


To simplify observation and extend the range of the first series of tests, 
three suction pipes 0.5, 1.5, and 3.0 inches in nominal inside diameter were 
used interchangeably. Each of these in turn was mounted vertically above an 
open tank of water 3 feet in diameter and 3 feet deep in which the elevation of 
the water surface could be regulated by means of supply and drain lines and 
measured by a hook gage (see Fig. 1). The inflow of air necessary to form 
the waterspout was produced by connecting the pipe to a large tank from 
which the air had been exhausted to a pressure of 10 inches of mercury abso- 
lute. 

As shown in the figure, a differential manometer was connected between a 
stagnation tube and a wall piezometer in the horizontal portion of the suction 
pipe for measurement of the centerline velocity at the gaging section, This 
was related to the mean velocity V = Q/}7D2 at the inlet section by means 
of the Bernoulli equation, through the assumption of adiabatic expansion be- 
tween the inlet and the gaging section. The pressure at the inlet was read on 
a barometer and the pressure at the gaging section was obtained from a 
Bourdon gage. The coefficient of the flow meter was determined by calibra- 
tion with water. 

With the water surface set at the desired elevation for a test, the control 
valve was slowly and continuously opened until the waterspout that developed 
just began to carry water into the pipe. At this instant the two gages were 
read at the flow-measuring section. The value in the suction line was then 
closed, the vacuum tank pumped down, and a different water level established 
in preparation for another test. 

In the second series of experiments a glass-walled tank 1 square foot in 
cross section and 2 feet high was used. This permitted visual observation 
of the interface between clear fresh water above and dyed salt water (specific 
gravity = 1.1) below. The arrangement of the suction pipe and vacuum tank 
was the same as before. The rate of flow, however, was determined by timing 
the fall of the free surface in the tank, which was quite feasible with the low 
rates required. Because of the limited size of the tank, only the 0.5-inch and 
1.5-inch suction pipes were used. 


The Results 


All experimental data are presented graphically in Fig. 2 in the dimension- 
less form suggested by the preliminary analysis. It can be seen that the tests 
using different sizes of pipe as well as entirely different fluids produced re- 
markably consistent results. Some scatter was to be expected, particularly 
for the 1/2-inch pipe in the two-liquid test, because of the difficulty in regu- 
lating the rate of flow to obtain critical conditions. The tendency of the points 
to follow a single curve (approximated by the power relationship 


y/D = 0.42[V/YDA¥/p ]*) indicates that the assumption of negligible viscous 
and capillary effects was justified for such an exploratory study. A 
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photograph of the waterspout in air is shown in Fig. 3, and one of the saline 
waterspout in fresh water is presented for comparison in Fig. 4. 

In conclusion, the limitations on the applicability of these results should 
be emphasized. The tests were run with fluids of low viscosity in relatively 
large tanks. Ripples, pendulations in the level of the interface, and the prox- 
imity of walls or other solid obstacles can be expected to increase the gap 
that is necessary to prevent formation of a waterspout. Thus, the results 
presented represent the minimum spacing required. It is evident that the 
value of y/D of 2 sometimes used in design is not sufficiently great to in- 
sure against the formation of a waterspout at high rates of siphonage. 
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Fig.2. Composite plot of measured data 
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Fig. 3. Surface profile of waterspout in air. 


Fig. 4. Saline waterspout (dyed) in fresh water. 


a 
=] 
1038-7 te 
ASCE ROUSE 
q 
i 
: 
q 
| 
| 
: 
= 


August, 1956 


PENETRATION OF A JET INTO A COUNTERFLOW 


J. Arendt,5 H. A. Babcock, ® and J. C. Schuster? 
(H. Rouse,l ASCE, Adviser) 


If a fluid jet is directed upstream against a uniform flow of the same fluid, 
it is to be presumed that the jet will penetrate some distance into the oncom- 
ing flow and then be turned back upon itself and carried downstream, as 
shown schematically in Fig. 5. The degree of penetration—in fact, the geome- 
try and kinematics of the entire flow pattern—could be expected to depend as 
a first approximation upon the diameter and velocity of the jet and the veloci- 
ty of the counterflow; that is, 


r, d, Yo) =0 


which may be reduced to the non-dimensional form 


03 +B) 


According to this relationship, two jets will yield similar patterns of flow 
(i.e., similar distributions of v/vg) only if the ratio vj/vo is the same for 
both, under which circumstances the distance upstream to the point of stag- 
nation (or the diameter of the enclosed jet flow, or any other characteristic 
length) will be a specific number of jet diameters. However, available infor- 
mation on submerged jets® would lead one to conclude that approximate simi- 
larity may be attained for jets of sufficiently small diameter and high velocity 
through the control of a single jet variable rather than two together. This 
variable is the momentum flux: 


= ind’ pw, 


Under these circumstances, 


»P,M)=0 


5. Assistant Professor, Civil Engineering Department, Fenn College, 
Cleveland, Ohio. 
6. Assistant Professor, Civil Engineering Department, Colorado School of 
Mines, Golden, Colorado. 
. Hydraulic Engineer, Bureau of Reclamation, Denver, Colorado. 
Albertson, M. L., Dai, Y. B., Jensen, R. A., and Rouse, H., “Diffusion of 
Submerged Jets,” Transactions A.S.C.E., vol. 115, 1950. 
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whence, in non-dimensional terms, 


x r Vv 


in which m=/M/p /y, . 


Experiments to permit an evaluation of this functional relationship were 
performed in an air tunnel having a 3x3-foot cross section and a velocity 
range from 15 to 45 feet per second. Measuring equipment consisted of a 
differential water manometer read by vernier to 0.001 foot, a hypodermic- 
needle stagnation tube, and calibrated Pitot cylinders with pairs of holes at 
30° and 180° from each other for determining the velocity magnitude and 
direction, The jet was produced by discharging compressed air from a cop- 
per tube of 0.015-foot inside diameter mounted at the axis of the tunnel test 
section approximately midway along its 12-foot length. 

Preliminary observations of the flow pattern were made by means of 
threads giving a visual indication of the distance that the jet penetrated into 
the counterflow and the radial extent of the zone of mixing. Once a sufficient 
length of penetration relative to the initial diameter of the jet had been estab- 
lished, measurements of the velocity distribution were made at a series of 
cross sections. Whereas a proper investigation would require the precise 
determination of both the magnitude and the direction of the velocity vector, 
in the limited time available it was possible to obtain only the magnitude. 
Even this, however, permitted an approximate evaluation of the flow pattern 
by the following technique. 

Assuming as a first approximation that v, = v, diagrams are prepared of 


27r Vx as a function of r , as in Fig. 6. Since Q = an/r Vx dr , use of a 


planimeter then permits the preparation of diagrams for Q as a function of 
r. Nowa stream tube (which is formed of streamlines) represents by defi- 
nition a surface through which no fluid passes. It follows that Q must be the 
same at successive cross sections of the same stream tube, and as well that 
the incremental rate of flow AQ between coaxial tubes must also be constant. 
It is hence possible to approximate the location of streamlines in an axial 
plane by noting the values of r for the same flow rates 0 , 4Q, 24Q, etc., 
at successive longitudinal distances x. From the direction of flow indicated 
by these streamlines, a second approximation may be obtained, which is usu- 
ally sufficient. 
The experimental results were evaluated in this general fashion, the sole 

difference being that non-dimensional rather than dimensional notation was 

used throughout. In other words, 7 quantity 27rv/mvg was integrated 
with respect to r/m to 5 Q/m? vo , and for the construction of stream- 
lines the increment A Q/m2vg = 0.2 was chosen, The resulting streamline 
configuration is shown in Fig. 7. This may be presumed to indicate, at least 
as a first approximation, the pattern of flow that would be produced by any jet 
of high enough relative velocity and small enough diameter for the latter to 
be inappreciable in comparison with the length of penetration (i.e., xs/d must 
be greater than perhaps //25). Under such circumstances, it follows that, 


since x,/m = -2.7, m=) M/? /vo , and M= tnd2pv;2 


x n Vv 
> - 2.7 | 
tT q 


4 
= 
: 
We 
pe 
Ra 
Jay 


1038-10 HY 4 August, 1956 


Although the measurements leading to the pattern shown in Fig. 7, lacking 
as they were in both precision and completeness, must be regarded as ex- 
ploratory only, they point at once to fruitful directions of more extended 
study. Three are readily apparent: first, verification of the assumption that 
m is a widely useful characteristic of jet performance by tests under a con- 
siderable range of conditions; second, determination of the pattern of second- 
ary flow—i.e., turbulence—in the same generalized manner; third, evaluation 
of the effect of viscosity in terms of a Reynolds number such as vgm/y , and 
perhaps eventually a completely analytic study of the pattern for the limiting 
case of laminar flow. 
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point 
countertlow and jet 


fig.5. Definition sketch 
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fig.6. Determination of streamline locations 


fig. 7 Generelized representation of 
the flow pattern 
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SURFACE PROFILES AT A SUBMERGED OVERFALL 


L. F. Ingram,? R. E. Oltman,!9 A, M, ASCE, and H. J. Tracy! 
(E. M, Laursen,12 A, M. ASCE, Adviser) 


Depending upon the interrelationship of a number of variables, the free 
surface at a submerged overfall will assume one of several basically differ- 
ent forms. In addition to the coordinate position, the pertinent variables in- 
clude the velocity and depth of the oncoming flow; the depth of the tailwater in 
the forebay; the channel slope, roughness, and cross-sectional dimensions; 
and the dimensions of the forebay. For the purpose of this exploratory study, 
the channel was made smooth and horizontal and the drop as great as possi- 
ble, changes in boundary alignment were limited to the vertical, and the sur- 
face configuration was observed in the crest vicinity only. The function under 
consideration was thus reduced to one between the following major 
parameters: 


9(5, 


Herein x and y represent horizontal and vertical distances from the over- 
fall crest, yg and yz are the initial and tailwater depths, and F is the 
Froude number gyo. 

The experimental investigation of this function was conducted in a glass- 
walled flume 1 foot wide, 2 feet deep, and 10 feet long. Falsework in the 
flume created a 1-foot drop, the horizontal portion of the channel being 2.5 
feet long and having a quarter-circle transition at its upper end. Discharge 
through the flume was measured by a calibrated elbow meter and an air-water 
manometer. The depth of the approaching flow could be controlled by a sluice 
gate at the beginning of the horizontal channel and measured by a point gage 
at a section 0.6 foot (i.e., 4yg) upstream from the drxp. The depth of the tail- 
water was controlled by vertical slats at the end of the flume, 7.5 feet beyond 
the drop, and was measured by floor piezometers connected to open manome- 
ter tubes. The water-surface profile was measured with a point gage gradu- 
ated to 0.001 foot and mounted on a carriage sliding on level stainless-steel 
rails. 

The first step in the standard operating procedure was to adjust the control 
valve and sluice gate until the desired rate of discharge and depth of approach 
were obtained. In all runs the approach depth was maintained at 0.15 foot, 


9. Hydraulic Engineer, Waterways Experiment Station, Vicksburg, Missis- 
sippi. 
10. Chairman, Branch Training Committee, U. S. Geological Survey, Wash- 
ington, D. C. 
11. Hydraulic Engineer, U. S. Geological Survey, Atlanta, Georgia. 
12. Research Engineer, lowa Institute of Hydraulic Research, State Universi- 
ty of Iowa, Iowa City. 
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and Froude numbers of 1, 2, 3, and 4 characterized the different series of 
runs, In each series successively greater tailwaters were established and 
the water-surface profiles measured. Except for the maximum tailwater, 
values of y¢/yoQ which were either positive or negative integers were chosen. 
The maximum tailwater used was that which formed a hydraulic jump just 
upstream from the drop; any slight increase in the tailwater above this value 
resulted in submergence of the sluice gate. Observations were also made of 
the tailwater elevation corresponding to the transitions between a plunging 
and a riding nappe and between a riding nappe and a hydraulic jump; because 
of something akin to a hysterisis effect, these observations were made with 
both a rising and falling tailwater. 

The water-surface profiles as measured are presented in dimensionless 
coordinates in Fig. 8. A more graphic representation of the typical flow pat- 
terns can be obtained from Fig. $, wherein photographs of the four runs with 
the highest tailwaters at F = 2 in Fig. 8 are shown, Three distinct regimes 
of flow can be seen, At the low tailwaters the nappe plunges below the sur- 
face, The tailwater elevation in this regime of flow affects the nappe prin- 
cipally insofar as it determines the pressure under the nappe: the lower the 
tailwater, the more the nappe is depressed, Above a certain tailwater eleva- 
tion, depending on the Froude number, a distinctly different flow pattern ob- 
tains in which the nappe does not plunge, but rides the surface downstream 
from the drop. As the tailwater elevation is further increased, the surface 
waves which characterize this regime of flow become progressively larger 
until the leading wave begins to break. Above this second critical tailwater 
elevation, a hydraulic jump forms with its toe slightly upstream from the 
drop. A small increase in the tailwater elevation is then sufficient to cause 
the jump to move entirely into the channel. For Froude numbers appreciably 
less than 2 the jump is of the undular rather than the breaking type, and the 
upper regime becomes indistinct from the intermediate. 4‘ a Froude number 
approaching unity, at which limit the sluice gate could no longer be used as a 
control, it is obvious that even the undular form of jump could not occur; tail- 
water elevations appreciably above the line of critical depth then resulted in 
complete submergence of the overfall and consequently in flow with a Froude 
number less than unity. 

The relative tailwater elevations at which the transitions between flow 
regimes occurred are shown in Fig. 10 as functions of the Froude number. 
Although no instability of the nappe was noted, there is evidently a significant 
difference between the values for rising and for falling tailwater. Because 
the flow tended to be stable, the critical elevations for falling tailwater are 
invariably below those for rising tai\water. 
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Fig 8. Measured variation in surface profile with 
Froude number and depth of tailwater 
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Fig. 9. Photographs of surface profiles for F = 2. 
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fig.!0. Lines of demarcation between 
regimes of flow 
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COEFFICIENT OF CONTRACTION FOR A SUBMERGED JET 


D. D. Curtis,13 M. ASCE, J. E. Martinez,!4 J, M. ASCE, and V. Vazquez15 
(H. Rouse,! M. ASCE, Adviser) 


In computing the rate of discharge for a submerged orifice, it is custom- 
ary to assume that the coefficient of jet contraction is the same as it would 
be if the jet were free. This assumption apparently yields satisfactory re- 
sults, else it would not have been made so consistently for the past century 
or more, It is not conceivable, however, that the flow patterns of the free 
and submerged jets would be identical or even closely similar. The present 
experiments were hence undertaken to compare both the geometric and the 
kinematic characteristics of a free and a submerged jet issuing from the 
same orifice under the same differential head. 

The tests were conducted in a glass-walled flume having a length of 25 
feet, a width of 2.5 feet, and a depth of 3 feet. A 2-inch sharp-edged orifice 
was placed 1.25 feet above the floor of the flume in the middle of a vertical 
steel plate separating the flume from a 6-foot stilling tank at its upstream 
end. Closure of the flume at its downstream end by a plate 2 feet high pro- 
vided not only a skimming weir for maintaining a constant depth of submer- 
gence but also a volumetric tank for preliminary determination of the head- 
discharge relationships. A gage carriage on rails, point and hook gages 
reading by vernier to 0.001 foot, a differential manometer with the same pre- 
cision of indication, a hypodermic-needle stagnation tube, and a stop watch 
were the measuring instruments used. 

Figure 11 shows the results of the calibration runs for both the free and 
the submerged jet, from which it will be seen that essentially the same coef- 
ficient of discharge—and hence essentially the same coefficient of contrac- 
tion—applies to both cases under the same differential head. The actual 
contraction of the free jet under a head of 2 feet was determined by measur- 
ing the surface profile in the vertical plane, computing the elevation of the jet 
axis therefrom, and eliminating the effect of gravitational deflection by sub- 
tracting from the measured profile elevations the corresponding deviation of 
the axis from the horizontal. The superposed results are shown in non- 
dimensional form in Fig. 12, together with computed values for the case of 
irrotational flow.16 Although the profiles differ somewhat in form, the limit- 
ing ordinates are very nearly the same, from which it is found that 
Ce = (0.78)2 = 0.61. 


13. Professor and Head, Department of Mechanics, Clemson Agricultural 
College, Clemson, South Carolina. 

14, Assistant Professor, University of New Mexico, Albuquerque, New Mexico. 

15, Assistant Professor, Monterrey Institute of Technology, Monterrey, N. L., 
Mexico. 

16. Rouse, H., and Abul-Fetouh, A. H., “Characteristics of Irrotational Flow 
through Axially Symmetric Orifices,” Journal of Applied Mechanics, Vol. 

17, No. 4, 1950. 
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From the stagnation-tube traverses of the submerged jet a family of 
velocity-distribution curves, selected ones of which are shown in Fig. 13, was 
obtained. From this sequence it is seen that the shear around the periphery 
of the jet both decelerates the flow within the jet and brings into motion the 
fluid immediately surrounding it; the region in which this effect occurs 
necessarily expands steadily with distance from the orifice. The extent of 
the jet contraction must, as a result, be defined rather arbitrarily. Now the 
most logical designation of the jet profile is that within which the rate of flow 
is the same at successive sections. The quantity 27rv at each section was 
hence plotted against the radius r , the integral Q= 22/rvdr of the area 
enclosed by each curve thereafter being plotted against r. From each of 
these curves the magnitude of r for which Q equals Qo , the rate of orifice 
discharge, could readily be determined, A series of the resulting values will 
be seen in Fig. 14, together with the profile of the free jet from Fig. 12. 

Comparison of the trend of the points with that of the curve will show that 
the profile of the submerged jet attains minimum proportions approximately 
a half diameter from the orifice but then gradually expands, whereas the free 
jet—although it theoretically approaches its limiting diameter asymptotically— 
actually undergoes little further change beyond a half diameter from the ori- 
fice. The minimum dimension of the submerged-jet profile, moreover, is 
somewhat larger than that of the free jet, the coefficient of contraction being 
correspondingly higher—i.e., Cc. = (0.79)2 = 0.63. Since the discharge coef- 
ficients of the two jets have already been seen to be essentially the same, one 
is forced to conclude that the rate of flow through the orifice is essentially 
independent of the tangential forces exerted upon the jet by the fluid body into 
which it emerges. In other words, the action of the turbulent shear in causing 
the submerged jet to expand has no appreciable effect upon the pressure field 
that controls the rate of efflux. This is also evident from the experimental 
indication in Fig. 13 that the velocity (and hence the pressure) in the irrota- 
tional core of the jet remains unaffected by the shear in the mixing zone. 
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free jet 
° Submerged jet 
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Fig./2. Comparison of measured and 
computed profiles of free jet 
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G04 Q06 Q08 Of 
r (ft) 
f71g./3 Variation of velocity distribution with distance 
along axis of submerged jer 
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fig./4. Comparison of profiles of free 
and submerged Jers 
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MAXIMUM AMPLITUDE OF OSCILLATION OF A 
CYLINDER DUE TO ITS VORTEX WAKE 


B. Guyton,17 A, M. ASCE, B. W. LeTourneau,18 and R. L. Player, Jr.19 
(L. Landweber,20 Adviser) 


INTRODUCTION 


When a circular cylinder is situated with its axis normal to a stream, ed- 
dies are shed periodically from alternate sides of the cylinder, generating a 
so-called Karman vortex trail.21 This gives rise to a fluctuating velocity 
field and hence also a fluctuating pressure field which exerts a periodic 
lateral force on the cylinder. 

The frequency of generation of the eddies by a non-vibrating cylinder has 
been determined experimentally.22 It was found, however, that, when the 
cylinder is allowed to vibrate, the eddy frequency and other properties of the 
wake are strongly affected by the oscillation;23 i.e., the oscillation modifies 
the force which causes it. 

An approximate analysis of the vortex-induced vibration of a cylinder, tak- 
ing into account the interaction between the vibration and the wake, was given 
by Steinman,24 who devised an ingenious graphical method for displaying the 
modification of the vibration frequency due to the vibration. The section of 
his paper entitled “Effect of Vibrations on Wake and on Eddy Frequency” con- 
cludes with the statement 


“In these vortex-induced oscillations, there is no ‘catastrophic 
range’ of increasing amplification with unlimited increase of stream 
velocity. The greatest amplitude attainable is little more than the 
diameter, and may occur only in a limited critical range just above 
the reasonance velocity.” 


17. Hydraulic Engineer, Waterways Experiment Station, Vicksburg, Missis- 
sippi. 

18. Atomic Power Division, Westinghouse Electric Corporation, Pittsburgh, 
Pennsylvania. 

19. Department of Civil Engineering, Duke University, Durham, North Caro- 
lina. 

20, Research Engineer, Iowa Institute of Hydraulic Research, State University 
of Iowa, Iowa City, Iowa. 

21. Karman, Th. von, and Rubach, H., “Ueber den Mechanismus des Flussig- 
keits- und Luftwiderstandes, Physikalische Zeitschrift, January, 1912. 

22. Relf, E. F., and Simmons, L.F.G., “The frequency of the eddies generated 
by the motion of circular cylinders through a fluid,” British Advisory 
Committee for Aeronautics, R. and M. 917, 1924. 

23. Landweber, L., “Flow about a pair of adjacent, parallel cylinders normal 
to a stream,” David Taylor Model Basin Report 485, July, 1942. 

24, Steinman, D. B., “Problems of Aerodynamic and Hydrodynamic Stability,” 
Proceedings of the Third Hydraulics Conference, State University of lowa 
Studies in Engineering, Bulletin 31, 1946. 
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Steinman defines amplitude as the maximum displacement of the center of 
the cylinder from its undisturbed position. However, he furnishes no evi- 
dence for his statement concerning the greatest attainable amplitude. 

The purpose of the present study is to enquire further, by both analytical 
and experimental means, into the question of the maximum amplitude of a 
vortex-induced vibration. 


Properties of the Karman Vortex Trail 


In the range of Reynolds numbers in which laminar separation occurs on 
the cylinder, it is found experimentally23 that the frequency n of vortex 
generation by a vibrating cylinder is given by 


nD _ 0,185 
(1) 


where D is the diameter of the cylinder, V is the free-stream velocity, 
and A is the ratio b/bg of the width of the vortex street with and without 
vibration. If a is the amplitude of the oscillation, we have the experimental- 
ly obtained relation 


A= 1+ 1.545 (2) 
The force per unit length F acting on the cylinder is an oscillating quan- 
tity whose amplitude is proportional to pV2h , where / is the density of 


the fluid; or, since b = Abg and since, for the case of laminar separation, it 
is found experimentally that bo = 1.3 D, it follows that 


F= (3) 


where § is a constant of the order of unity. 


Analysis of Cylinder Vibration 


The vibration of the cylinder will now be analyzed on the assumption that 
it is a linear harmonic oscillator satisfying the equation 


2 
+¢ ky = F cos 2nnt (4) 
dt 


where m is the apparent mass per unit length of the cylinder (actual mass 
plus virtual mass), c is a damping coefficient, and k a stiffness coefficient. 
When the vibration is steady, at frequencies near resonance, there results 
for the maximum amplitude 


(5) 


or, from Eqs. (1) and (3), 
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ay = 0,86 


But, from Eq. (2), 
&, = 0.65 D( Any 1) (7) 


Hence, upon equating the two expressions for ayq in Eqs. (6) and (7), one ob- 
tains a quadratic equation in A yw whose solution is 


The largest permissible value of 5.38 pVD/c is seen to be unity, correspond- 
ing to which A y = 2 and, from Eq. (7), ayg = 0.65 D. This indicates that the 
maximum attainable amplitude is about 2/3 of the diameter. 


Experiments 


An attempt was made to support experimentally either Steinman’s state- 
ment or this simplified theory. Since the theory as developed is based on the 
two-dimensional case of a rod vibrating without deformation, experimental 
confirmation should have been sought with a stiff cylinder supported by 
springs or elastic cantilever beams which would vibrate at right angles to 
the stream without deformation. However, time did not permit the design and 
construction of such an apparatus, and experiments were conducted instead 
on the vibration of a string. Since the string is flexible and deforms in a 
combination of characteristic harmonic shapes when vibrating transversely, 
and in addition is bowed in the direction of flow, the analysis presented above 
is no longer truly applicable, although it may serve as a first approximation. 
No analysis for the case of the string was attempted. 

For the assumed condition that the foregoing theory was applicable to the 
string, resonance would occur when the natural frequency of the string was 
equal to the frequency of generation of the vortices—i.e., 


0,185 KD Bs n 2, 3, 


where T is the tension in the string and s its length. 
Attempts were made to secure resonance using both 1/16-inch piano wire 
and 0,054-inch nylon string in a small air tunnel, but without success. This 
may have been due to the fact that the lowest fan speed available resulted in 
an air velocity in excess of 20 feet per second. Resonance was achieved, 
however, by stretching the same nylon string across the 14-inch diameter 
test section of the Institute water tunnel. The velocity was varied from 0 to 
15 feet per second, and the vibration of the string was clearly observed 
through a glass window. The fundamental mode and the first and second har- 
monic shapes, respectively, were clearly distinguishable as the velocity was 
increased, 
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The velocity was measured with a mercury manometer reading differen- 
tial pressure across the bell entrance to the test section. The amplitude of 
vibration, relative to the string diameter was determined by measuring with 
a surveyor’s transit the distance between the upper and lower extreme loci 
of the string as it vibrated. The transit was located on the same level in or- 
der to eliminate extraneous drag vibrations, and the tripod was rubber- 
insulated from the vibrations of the floor. The frequency was measured by 
observing the vibrations with a Strobotac and recording the highest synchro- 
nized frequency. The tension of the string was measured with the water at 
rest by pulling normally to the string center with a spring balance, applying 
a known force, measuring the resultant sag, and calculating the stress by the 
usual methods of elementary statics. It was found difficult to apply a large 
tension to the nylon, because the tied end connections had a tendency to slip. 


RESULTS AND CONCLUSIONS 


The experimental results are presented graphically in Fig. 15, in which 
the dimensionless amplitude a/D and the frequency parameter or the 
Strouhal number nD/V are plotted against the Reynolds number VD/V¥. The 
Strouhal numbers observed were of reasonable but somewhat high magnitude 
and decreased slightly with the Reynolds number, as would be predicted by 
the relation nD/V = 0.185 and the amplitude data. The amplitude data showed 
considerable scatter, some of which was due to the rudimentary instrumenta- 
tion and some to the inherent nature of the experiment. The amplitude rose 
rapidly from zero at very low speed to a range of a/D from 0.25 to 0.55 at 
the fundamental resonance Reynolds number of 3000 predicted by using a cal- 
culated fundamental natural frequency of the string of 212 cycles per second 
corresponding to the 8.4-pound tension. The maximum observed dimension- 
less amplitude was 0.58 at velocities near fundamental resonance, A mean 
curve is plotted through the data. The results tend to support the foregoing 
simplified theory, since the predicted and observed maximum relative ampli- 
tudes are quite close: 0.65 and 0.58, respectively. 

The visually observed modes of vibration are also indicated on the graph. 
It is to be noted that, once the maximum amplitude was attained, it remained 
relatively constant and independent of the mode of vibration as the velocity 
was increased, This indicates a sort of sustaining resonance, as was graphi- 
cally explained by Steinman.2 The transition between vibratory modes was 
not smooth, but there were ranges of velocity when the vibration was observed 
to “jump” from one mode to another, 
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DEFLECTION OF A JET BY A NORMAL WALL 


R. L Murray,25 L. E. Nelson,26 M, ASCE, and J. R. Steven,2? J. M. ASCE 
(P. G. Hubbard,28 J, M. ASC#, Adviser) 


Although the deflection of a submerged jet by a normal wall is closely re- 
lated to practical problems such as the scouring tendency of a plunging nappe, 
it is of interest primarily as an example of free turbulent diffusion and mo- 
mentum transfer. The investigation of this phenomenon had been preceded at 
the Iowa Institute by a rather complete study of the flow pattern which is 
created by discharging a jet of fluid from an opening in a plane boundary into 
a semi-infinite reservoir of the same fluid.29 In the present case, the object 
was to determine the streamline pattern which results when the flow is de- 
flected by a wall normal to the initial direction of the jet. 

The experimental equipment for the study consisted of two carefully 
rounded nozzles 0.025 and 0.05 foot in diameter which could be mounted, in- 
terchangeably, flush with the surface of an octagonal table 30 inches across. 
Air was supplied to the 6-inch pipe preceding the nozzle by a centrifugal fan, 
and the discharge was controlled by a variable-area intake to the blower. A 
masonite board mounted 2.5 feet above the table provided the normal wall for 
deflection of the jet. Velocities throughout the region of interest were mea- 
sured by means of a hot-wire anemometer mounted on a mechanism for mak- 
ing longitudinal and lateral traverses. The pressure distribution on the nor- 
mal wall was measured by an alcohol-kerosene differential manometer which 
could be connected to a series of piezometers arranged radially outward from 
the stagnation point. 

Because the flow system was axially symmetrical, it could be described 
completely in terms of the velocity v(x,r) or streamlines y(x,r) ina 
meridian plane, x being measured from the plane of the issuing jet and r 
radially outward from its axis. From the available results of the earlier jet 
studies already mentioned, it was known that the determinants of the flow pat- 
tern would be the nozzle diameter D , the momentum flux of the jet 

M = Qpvj , the fluid density p , and the distance a between the bounding 
planes. In order to reduce the scope of the study to correspond to the time 
available, a was made so large compared to the jet diameter that the jet 
could be considered as a point source of momentum. 

A dimensional analysis of the problem, beginning with the variables 


25. Assistant Professor, College of Engineering, University of Santa Clara, 
Santa Clara, California. 

26. Hydraulic Engineer, U. S. Engineer Office, Kansas City, Missouri. 

27. Assistant Professor, Department of Civil Engineering, College of the City 
of New York, New York. 

28. Research Engineer, Iowa Institute of Hydraulic Research, State University 
of Iowa, Iowa City. 

29. Albertson, M. L., Dai, Y. B., Jensen, R. A., and Rouse, Hunter, “Diffusion 
of Submerged Jets,” Trans. ASCE, Vol. 115, 1950. 


ER 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
] 
| 
| 
| 


1038-27 
already defined, can be outlined as follows: 


v= f(x,r, a, 9, M) 


or, in dimensionless terms, 


Wipe * 5) 


An equally useful and more easily visualized method of presenting the re- 
sults is to plot streamlines in an axial plane, using the definition 


y rdr=Q,/2n 


where vx is the axial component of velocity and Q, is the discharge 
through a circle of radius r. A pattern of streamlines can be obtained from 
this relationship by evaluating graphically the radii corresponding to equal 
increments of w and then plotting curves of constant y or of the dimen- 
sionless value 


Because the hot-wire anemometer used in measuring mean velocities was 
oriented so that its axis was normal to an axial plane, it indicated the result- 
ant velocity v rather than the axial component vx required in the computa- 
tions for the streamlines. This meant that two or more successive steps had 
to be taken in determining the location of the streamlines in the region where 
they were not parallel to the axis of the jet. The first approximation was ob- 
tained by assuming the measured velocity to be axial, integrating graphically, 
and plotting a streamline pattern on this basis. With 6' as the angle between 
the axis of the jet and the tangent to one of these streamlines, then a second 
graphical integration could be made using v' = v cos @' . This process con- 
verged rapidly to the correct pattern, and even the second approximation was 
unnecessary except where the streamline curvature was pronounced. 

Near the wall, of course, the streamlines curved rapidly and became prac- 
tically radial where r/(a-x) >3. For this region a modified method of locat- 
ing streamlines was used. If vr represents the radial component of velocity 
and Qa- x the rate of flow through a cylindrical surface of radius r and 
length a - x , then the stream function can also be defined as 


This takes the dimensionless form 
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a= 
a 


M/p a2 


Where correction of the values for inclination of the streamlines is neces- 
sary, the process is similar to the one already outlined, except that the 
quantity v' = v sin 9' is substituted in the numerator. In the region 

x/a > 0.8, r/a<0.3, both methods were used as independent checks upon the 
accuracy of integration. 

The final result for this study is shown in Fig. 16. Identical curves were 
obtained using each of the two available nozzles, so that the assumption of a 
point source of momentum was justified. In the latter circumstance the 
stream tube enclosing the discharge from the nozzle becomes so small that 
it cannot be shown on such a plot, the flow pattern hence seeming to consist 
almost entirely of streamlines for the fluid entrained by turbulent diffusion. 
The dashed line indicating the nominal boundary of the jet has the slope 1:5, 
which applies to a jet without the normal wall, and which also defines the 
limits in this case quite well. 

As is evident from the diagram, the effect of the wall is felt primarily in 
the region x/a > 0.8. In spite of the concentration of the flow about the axis, 
the measured form of the streamlines in this region is almost identical with 
that computed on the assumption of irrotational flow. For the latter, the 
streamline equation is? 


= mr*(a-x) 


wherein the coefficient m for the present experiments is approximately 


52 . 

Dimensionless values of the pressure measured on the normal wall are 
also plotted in Fig. 16, together with that obtained from the foregoing rela- 
tionship for irrotational flow with a stagnation value of 27 M/a2. Integra- 
tion of this pressure distribution yielded the unexpected result that the total 
force on the wall is approximately 4/3 as great as the initial momentum flux 
of the jet rather than exactly equal to it. This apparent dynamic unbalance 
was probably due to the pressure reduction on the lower boundary accompany- 
ing the radial acceleration of the air drawn in by the jet. Because this nega- 
tive pressure occurred over a relatively large area, its local intensity was 
sO small that it could not be measured reliably. 

On the basis of the experimental results of this study, two important con- 
clusions can be drawn concerning the deflection of a jet by a normal wall. 
First, the wall has a negligible effect on the flow in an initial zone having a 
radius of at least one-fifth and a length of roughly four-fifths the boundary 
spacing. Second, in the region of pronounced streamline curvature near the 
stagnation point on the wall, the flow pattern may be assumed to agree with 
that of irrotational motion over a radial distance of about a/8. 


30. Forchheimer, P., Hydraulik, Leipzig and Berlin, 3rd Ed., 1930, pp. 28 
and 29, 
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INFLUENCE OF SLOPE AND ROUGHNESS ON THE FREE OVERFALL 


J. W. Delleur,31 J. M. ASCE, J. C. I. Dooge,32 A, M. ASCE, and K. W. Gent33 
(E. M. Laursen,12 A. M. ASCE, Adviser) 


This study is concerned with the depth of flow which occurs at the end of a 
long channel discharging freely over a fall under various conditions of chan- 
nel slope and bed roughness. The case of a horizontal channel ending in a 
free overfall was investigated in detail by Rouse (Civil Engineering, April 
1936). The effect of channel slope was studied by Fathy and Amin (A.S.C.E. 
Proceedings, Separate No. 564, December 1954). A listing of a number of 
publications on various aspects of the free-overfall problem appeared in the 


February 1955 issue of Civil Engineering. 


Momentum Analysis 


By examining flow conditions, assumed to be two-dimensional, between the 
point of overfall and an upstream point x at which the pressure distribution 
may be taken as hydrostatic (see Fig. 17), one can express the momentum 
equation for a direction parallel to the bed in the following form: 


= (1) 


If 6 is taken as positive downward in the direction of flow, the above equation 
will apply to both downward and adverse slopes. Upon evaluating each term 
in turn, one obtains 


=/ pdy = K, 7-3 cos @ 
where p dy /(}y,*cos 6/2); 


Ms ef V? dy = Ke 


where -["v" du ys) s 
R=4 3 cos 6 


31. Assistant Professor of Civil Engineering, Purdue University, Lafayette, 
Indiana. 

32. Senior Design Engineer, Hydrometric Section, Civil Works Department, 
Dublin, Ireland. 

33. Hydraulic Engineer, H. G. Acres & Co., Ltd., Niagara Falls, Ontario, 
Canada, 
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since pressure is hydrostatic; 


Sx 
W, sind = 7 sin = Ky sind yyt 


where = [ dx/ (y,L) ; and, finally, 


where S¢ is the resistance slope at depth y , Sy is that at point x , and 
Ks = y Equation (1) now becomes 


2 
2 


Division of each term by ¥ and elimination of velocities by the continuity 
equation = Vx¥x = 49 cos@ then yields 


2 yx? | 
K, cos + cos = cos 6 +. cos @ + y,L(K, sinO-KsS,) 


This equation is quite general, being identical with Eq. (1). 

If the slope does not exceed 1:10, cos @ can be taken as unity and sin 0 
written as Sg. Furthermore, if the flow is turbulent, Kg will be slightly 
greater than 1.0 and Ko will be between 1.0 and Kg , because the flow is 
convergent. Since these terms are on opposite sides of the equation, one can 
take Ko = Kg = 1.0 without appreciable error, These two changes give 


Ys? Ye? yx” 


If Sg and S, are not too large, it is reasonable to assume (subject to ex- 
perimental verification) that the last term in Eq. (4) can be neglected to give 


(5) 
Ye 2 Ux 


This is applicable to turbulent flow on both positive and adverse slopes under 
various roughness conditions. For steep slopes of sufficient length the depth 
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of the approaching flow is equal to the normal depth, so that 


us? Ue? 
=> + = 
Yo 


which can be written as 


2K YAP 


For mild or adverse slopes yy, can be taken as equal to y¢ , so that Eq. (5) 
becomes 


at = 2ix, (4e)? 


Experimental Equipment and Procedure 


The experimental investigation was carried out in a glass-walled tilting 
flume 30 feet long and 2 feet wide. Discharges were measured by means of 
an orifice meter and were selected to give simple numerical values of the 
critical depth for parallel flow. In the preliminary tests seven values of the 
critical depth ranging from 0.1 to 0.4 foot were used, while the final tests 
were confined to discharges corresponding to 0.2 and 0.4 foot. Ten slopes in 
the range +0.03 to -0.02 were tried. Two roughness conditions were investi- 
gated: (a) the relatively smooth stainless-steel bottom of the flume, and 
(b) an artificial roughness consisting of 1/4x1/4-inch brass bars fixed trans- 
versely along the bottom of the flume. Previous tests by Walter Rand of the 
Institute staff had indicated that with this artificial roughness in position the 
virtual bottom of the channel was at the mid-height of the roughness bars. 
Accordingly, the crest of the overfall was arranged as in Fig. 18, so that the 
crest level would be in line with the virtual bottom of the approach channel. 

The standard experimental procedure was as follows. The discharge was 
set to give the required value of the critical depth for parallel flow and kept 
constant while the slope was varied throughout its range. For the horizontal 
case and for the maximum positive and negative slopes, enough point-gage 
readings were taken to determine the complete surface profile. For the in- 
termediate slopes, enough readings were taken to determine yr , L , and the 
slope of the water surface at the overfall crest and at the point x. 

The scatter in the results obtained seemed to be within the limits to be ex- 
pected for the experimental set-up involved. The fluctuations in the air- 
water manometer used in measuring the discharge could easily result in 
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an error of 2% or more in the discharge and hence a corresponding error in 
the value of the critical depth and critical slope. For flows near the criti- 
cal and for the steeper slopes the water surface was quite disturbed and the 
depth as determined by point-gage readings could be in error by 1 or 2%. 


DISCUSSION OF RESULTS 


The experimental results can be first used to justify the dropping of the 
last term in Eq. (4) - namely, L y,(K4Sg-KsS,). From the detailed profiles 
it was found that, as yg increased towards yx , K4 varied from 0.90 to 0.98 
and Ks varied from 1.25 to 1.05. It was noted that the following approxima- 
tions based on a parabolic profile gave good results: 


3 


For steep slopes the term under examination is L yx Sg (K4 - Ks), which 
varies from -0,015 Vx" to -0,005 yx". For mild slopes its value varies from 
-0.02 Vx" to -0.04 Vx? and in sixteen of the twenty cases measured lies be- 


tween -0.025 yx? and -0.035 yx?. The term is clearly of the same order (and 
opposite sign) as the approximation involved in taking K3 = 1 and may be 
neglected in a first approximation to the solution. 

The critical slope Se for two-dimensional flow is given by 


In evaluating the critical slope in the present tests, the value for the resist- 
ance coefficient f for the rough case was taken from an equation derived by 
Walter Rand for the same artificial roughness placed in the same flume, 
while for the smooth case C was taken from a formula proposed by R. W. 
Powell. The two formulas are 


Rough: = 2.04 log +135 


Smooth: C= 42 log B 


The values of yg measured in the final series of tests are plotted in 
Fig. 19 interms of yf/yc and Sg/Sc. It will be seen that the points repre- 
senting widely varying conditions of roughness and slope define a single line. 
This is to be expected in the case of steep slopes, as the plotting is in terms 
of the factors occurring in Eq. (6). The valu2s for mild and adverse slopes 
also plot fairly well about a line, but the small differences involved make it 
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difficult to determine whether the scatter which occurs is systematic or not. 
In the case of adverse slopes, less scatter is obtained if yg/ye is plotted 
against So +S¢ (which will be negative when Sg < - Sc), but further experi- 
mentation would be required to settle this question. Accordingly, the results 
for steep, mild, and adverse slopes are all plotted in Fig. 19 in terms of 
and So/S, . 

The corresponding relationship based on Eqs. (6) and (8) cannot be deter- 
mined unless values for the pressure variable Kj are known or assumed, 
The experiments indicated that Kj decreased in value as the slope passed 
from adverse to horizontal and from horizontal to mild and then to steep. A 
good correspondence with the experimental points (note the plotted curve) is 
obtained if the following values of Kj are assumed: 


< = 0.00 
<5,/8,<1 K, =0,30 + -/1-8,/8, /s 
K, = 0.80 


However, the validity of assuming these values must still be checked by 
direct pressure measurements. 
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(a) Positive slope (ob) Adverse slope 
Definition sketches 


Fig. 18 Details of bar roughness iy ] + 
and overtall crest 
Vb 
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Fig./9 Composite plot of test results 
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Paper 1039 


Journal of the 


HYDRAULICS DIVISION 


Proceedings of the American Society of Civil Engineers 


A NEW DEVELOPMENT IN FLOW MEASUREMENT: 
THE DALL FLOW TUBE 


Andre L. Jorissen,* M. ASCE 
(Proc. Paper 1039) 


Several new types of flow-metering devices of the pressure-difference or 
inferential type for pipe flow measurements have been developed in recent 
years. 

Of particular interest to civil engineers is the Dall Flow Tube which, as 
shown by Fig. 1, is a modified Venturi Tube. Originally built and studied in 
England! its characteristics were first discussed before an American audi- 
ence by I. O. Miner of Builders-Providence, Inc.2 

Essential differences between the Dall Flow Tube and the classical Venturi 
Tube are: 

1) The steepness of both the converging cone and of the diverging or 
pressure-recovery cone, the former having an included angle of 
approximately 40° and the latter an included angle of approximately 159°, 

2) The presence of a dam at (a) and the characteristic location of the up- 
stream pressure connection, similar to the type commonly used in 
Europe and known as “corner taps.” 

3) The sharp edges at (b), (c), (d) and (e). 

4) The slot between (c) and (d), through which the downstream or throat 
pressure connection is made. 

5) The sudden enlargement at (f). By truncating the recovery cone and 
because of the steepness of the converging and recovery cones, a very 
compact instrument is obtained, the length of the whole device being of 
the order of two pipe diameters only. 


Experimental data indicate that the differential pressure produced by a 
Dall Flow Tube is higher than that of a conventional nozzle or Venturi Tube. 
The ratio of Dall Flow Tube differential to Venturi Tube differential varies 
from 1.91 for a diameter ratio B of 0.53 to 2.43 for B = 0.8. 


Note: Discussion open until January 1, 1957. Paper 1039 is part of the copyrighted 
Journal of the Hydraulics Division of the American Society of Civil Engineers, Vol. 
82, No. HY 4, August, 1956. 

* Head, Dept. of Hydr. and Hydr. Eng., Cornell Univ., Ithaca, N. Y. 

1. By H. E. Dall, Research Engineer for George Kent, Ltd. 

2. I. O..Miner: “The Dall Flow Tube,” Trans. A.S.M.E., 78, no. 3, April 1956. 

3. Bis the ratio of throat diameter to pipe diameter: # = d/D. 
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An unexpected characteristic of the Dall Flow Tube is that the nonrecover- 
able pressure head, or head loss of the instrument is extremely low. Fig. 2 
gives a comparison, established by Miner, between the head loss, expressed 
in per cent of the differential, for the Dall Flow Tube and two types of Ven- 
turi Tubes. In the standard Herschel-type long form Venturi Tube, the in- 
cluded angle of the recovery cone is 5°5', In the short form Venturi Tube, 
this angle is appreciably larger. Furthermore, the recovery cone is trun- 
cated to obtain a more compact instrument. The solid lines on Fig. 2 indi- 
cate the head loss in per cent of differential as a function of the diameter 
ratio 8. The dotted line gives the loss through a Dall Flow Tube producing 
the same differential as that cof a Venturi Tube with the value of 8 given by 
the abcissa. 

The advisability of expressing the head loss in per cent of the differential 
may appear questionable because of the high value of the differential itself. 
Another comparison, based on the ratio of head loss in Dall Flow Tubes to 
head loss in Herschel type Venturi Tubes, long form for equal flow rates (at 
pipe Reynolds Numbers larger than 350,000) indicates that, for a given value 
of 8 , the head loss for the Dall Flow Tube is smaller when the diameter 
ratio is greater than 0,552. No satisfactory explanation has been given yet 
of the low pressure loss of the Dall Flow Tube. Several studies are in pro- 
gress, trying to obtain a better understanding of the fluid mechanics princi- 
ples involved in the flow through a Dall Flow Tube. 


Dall Flow Tube Discharge Coefficients 


The relationship between rate of discharge and differential head is ex- 
pressed in the most general case of compressible fluids by the classical 
equation for inferential flow meters: 


CYa Yy2g.h 


in which 
q = volumetric rate of 
a = throat diameter = a 


C = coefficient of discharge 

Y = expansion factor 

B = diameter ratio = d/D 

g = acceleration of gravity 

h = differential pressure head 


It is well known that the coefficient of discharge C is a function of the 
type and proportions of the instrument and of flow characteristics which com- 
bine in the expression of the pipe Reynolds number: 
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JORISSEN 
V = average velocity of flow in the pipe 


D = pipe diameter 
v = coefficient of kinematic viscosity of the fluid 


The expansion factor Y takes into account the compressibility of the 
fluid. It is a function of the diameter ratio 8, of the ratio of throat to inlet 
pressure and of the adiabatic exponent k, For incompressible fluids, Y has 
the value of unity. 

Values of C, as obtained by Miner from his experiments are given by the 
solid line on Fig. 3. Those values are for Reynolds numbers larger than 
350,000. 

The author has had the opportunity to calibrate several Dall Flow Tubes in 
the Hydraulic Laboratory at Cornell University. One of them was of large 
size: 24.02 in. x 17.544 in. (8 = 0.7304). This tube was calibrated by mea- 
suring flows between 5,62 and 20.50 cfs volumetrically in a standpipe 6 ft. in 
diameter and 55 ft. high providing a tank of 1500 ft3 useful capacity. A suit- 
able diverter equipped with an electrical timing device was used. Differen- 
tial pressure heads ranged from 0.282 to 3.80 ft. of water. The results of 
that calibration are given on Fig. 4. 

The curve shows the characteristic trend indicated by Miner, the coeffi- 
cient of discharge decreasing for increasing Reynolds number and reaching 
a constant value: Ce at high Reynolds numbers. 

The constant value of the coefficient of discharge is Cc. = 0.661, which 
agrees very well with Miner’s data from Fig. 3. 

Other calibrations were part of a research project on inferential meters 
conducted in the Hydraulic Laboratory at Cornell University with the spon- 
sorship of Builders-Providence, Inc. The tests were all conducted with 8 in, 
pipe size and the flows measured volumetrically. The results of those tests 
are summarized in Table I giving the dimensions of the Dall Flow Tubes 
tested, the value of Reynolds number at which the coefficient of discharge be- 
comes constant: (Rp)n and the constant value of the coefficient of discharge: 
Cc. The last three lines apply to the instrument known as the Dall Insert 
Nozzle which is a Dall Flow Tube inserted between flanges in the pipe line. 

It should be noted that the values of (Rp)p differ somewhat from the value 
of 350,000 indicated by Miner. The difference, however, is small except for 
two of the largest diameter ratios. 

A comparison of the experimental values of Cc; with the data from Fig. 3 
indicates a very satisfactory agreement except for the Dall Insert Nozzle with 
the largest diameter ratio. Part of this difference may be accounted for by | 
the fact that inferential flow meters of large diameter ratio are more sensi- 
tive to the conditions of installation than those of small diameter ratio. 


Effect of Two Short-Radius Elbows in Orthogonal Planes 


Miner has investigated the effect of various upstream disturbances on the 
behavior of the Dall Flow Tube, namely the effects of a single elbow, of a 
gate valve partly opened, and of an increaser or a decreaser. 

The author studied the effect of two short-radius elbows in orthogonal 
planes. 

The test run was an experimental pipe, 8 in. in diameter, under a static 
head of 80 ft. 

Tests were conducted both without and with straightening vanes at the 
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upstream end of the approach pipe, immediately after the two elbows. These 
straightening vanes consisted of four radial baffles approximately 16 in. (2D) 
long. 

The short-radius elbows themselves were preceded by sufficient lengths of 
straight pipe to insure that they alone would affect the behavior of the instru- 
ment tested. 

The results of the tests are summarized in Table IL. La is the length of 
approach pipe between the two elbows and the tube tested, D the pipe diame- 
ter, C the coefficient of discharge and a the difference (in per cent) between 
the coefficient of discharge for a particular condition and that obtained for the 
same tube with a length of approach pipe 27 ft. 7/8 in. (41D) long. 

It should be noted that any value of @ smaller than 0.2 per cent is insig- 
nificant, being smaller than the possible experimental error. 

The study indicates that the effect of two short radius elbows in orthogonal 
planes on the performance of Dall Flow Tubes is only slightly larger than 
their effect on standard Venturi tubes. The disturbance appears more signifi- 
cant for larger diameter ratios. The effect of straightening vanes is to raise 
the coefficient. 

Comparison of the values of the coefficient of discharge obtained without 
and with straightening vanes gives the following differences: 


0.15 per cent for the 8 in. x 4 in. Venturi Tube 
0.47 per cent for the 8 in. x 4.4 in. Dall Flow Tube 
0.31 per cent for the 8 in. x 5 in. Venturi Tube 
0.90 per cent for the 8 in. x 5.2 in. Dall Flow Tube 
0.73 per cent for the 8 in. x 6 in. Venturi Tube 


Tests are now under way to study the effect of two elbows on a Dall Flow 
Tube of very large diameter ratio. 


Expansion Factor Tests 


A determination of the value of the expansion factor Y is essential to the 
use of Dall Flow Tubes with compressible fluids. Since no theoretical analy- 
sis is available to this effect, the values of the expansion factor must be ob- 
tained experimentally. This was done by the author, using air flows. The 
procedure was as follows: 

Various Dall Flow Tubes were first calibrated in the Hydraulic Labora- 
tory, using water. The Dall Flow Tubes were installed in series with Venturi 
Tubes of known characteristics. Suitable lengths of approach pipes were pro- 
vided both for the Venturi Tubes and for the Dall Flow Tubes to eliminate the 
effect of approach conditions. After the water calibration, each installation 
was installed in an air compressor line. Since the expansion factors for 
Venturi Tubes are known‘ the Dall Flow Tube expansion factors could be de- 
termined by comparison, The results are given in Fig. 5. Values of the ex- 
pansion factor Y are plotted as a function of x/k. 


x = ratio of differential pressure to absolute inlet static pressure 
k = adiabatic exponent 
In each case, the Dall Flow Tube expansion factor may be compared with that 


4, A.S.M.E. Power Test Code, Part 5, Chapter 4. 
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of a concentric orifice of equal diameter ratio.5 The conclusion is that under 
the conditions investigated the values of the expansion factors for Dall Flow 
Tubes and concentric orifices are quite comparable. 


CONCLUSIONS 


The Dall Flow Tube, as now constructed, thus offers the two-fold advantage 
of a great compactness and a low head loss. It is only slightly more sensitive 
to the conditions of installation than standard Venturi Tubes. 

For Dall Flow Tubes with throat to inlet diameter ratios 6 not more than 
0.75 (usual limit for Venturi Tubes), the coefficient of discharge can be pre- 
dicted within 1%. For larger values of 8, direct calibration under upstream 
installation conditions is recommended, unless larger tolerances on coeffi- 


cient value are acceptable. 


5. A.S.M.E. Power Test Code, Part 5, Chapter 4. 
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DALL FLOW TUBE EXPANSION FACTOR 
FIG. 5 
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FLOOD PLAIN ASPECTS OF RIVER PLANNING 


Anthony M. Lunetta,* A.M., ASCE 
(Proc. Paper 1040) 


A proper evaluation of the role of the flood plain in river planning involves 
a consideration of all aspects of water resources. Such river problems as 
trar .portation, power, recreation, water supply, sewage disposal, pollution, 
flood control and sedimentation affect the welfare of the nation. A study of 
all the problems connected with and involving streams is beyond the scope of 
this paper. Instead, only one aspect of the flood control problem will be dis- 
cussed, namely, the use and control of the flood plains through the medium of 
zoning and the problems encountered. 

Regulating or prohibiting the use of the flood plain for dwellings or other 
property which would be seriously damaged by flooding has been proposed at 
various times since 1900. With the acceleration of house construction in the 
late 1940’s, the use of zoning to control the flood plain has increased in the 
United States. The concept of flood plain zoning has its greatest use as a 
preventive measure in areas which are undeveloped and sparsely developed. 
The presence of buildings and other structures of various types may make it 
more economically feasible to provide protection to the area rather than to 
evacuate the area. The nature of the zoning procedure takes many forms and 
is administered and enforced by numerous governmental units and agencies. 
The range of flood plain zoning varies from control of encroachments on the 
stream itself to the prohibition of the use of any portion of the flood plain for 
any activity or structure which will adversely affect flood flows. In some 
cases, where zoning powers were believed to ve insufficient to restrict the 
use of the flood plain, the areas were purchased as public lands and developed 
for parks and other recreational facilities. 

Encroachment on the flood plains of rivers and streams has occurred since 
the beginning of civilization. The earliest settlements, requiring water for a 
multiplicity of uses, such as transportation, potable water, protection and 
waste disposal, grew up along the banks of rivers, often on the flood plains. (1) 
These settlements established economic values and social roots which have 
withstood, in most instances, the destructive forces of floods. But floods are 
periodic and persistent, and necessitate action by man to minimize or prevent 


Note: Discussion open until January 1, 1957. Paper 1040 is part of the copyrighted 
Journal of the Hydraulics Division of the American Society of Civil Engineers, Vol. 
82, No. HY 4, August, 1956. 

* Hydraulic Engineer, Bergen County Engineering Department, Hackensack, 

New Jersey. 
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the damages caused by excessive storm runoff. Sometimes this action has 
been directed towards protecting these values and roots without determining 
if the protection, when all aspects are considered, is warranted. The net re- 
sult is a continuous contest of man opposing the forces of nature and in many 
instances the results are far from satisfactory. 

The expense involved, initially and subsequently, of providing the protec- 
tion against floods has reached the point where the wholesale evacuation of 
the flood plain is being considered as the most economical solution to the 
man-made problem of flood damages. Some areas, benefiting from the ex- 
periences of their neighbors, are attempting to forestall future flood prob- 
lems by maintaining as much of the flood plain as is economically feasible 
for the passage of flood waters. 


The Role of the Flood Plain 


The flood plain plays an important part in the development of a society. 
The environmental features of the flood plain lend themselves to manifold 
uses which have a direct bearing on the growth of cities and the development 
of the area. The proper uses of the flood plain varies from area to area and 
depends on the climate, geography and natural resources, as well as on the 
existing and potential development of the non-flooded area. The solution of 
stream problems is usually directed towards the maximum utilization of the 
flood plain. This may involve the inundation of agricultural lands by storage 
reservoirs on one part of the stream to protect highly developed industrial 
and residential areas situated on the flood plains of another portion of the 
stream. The social, psychological and economical problems arising from 
such a solution are exemplified in the controversy which resulted from 
the proposed Tuttle Creek dam on the Big Blue River in Kansas to protect 
the cities of Topeka and Kansas City. 

Planning, in its best sense, entails full knowledge of the problem, its 
causes, its symptoms and it cures. So with river planning, it is necessary to 
understand the behavior of streams and the characteristic features of their 
drainage basins in order to plan effectively for the maximum utilization and 
protection of the water resources available. In addition, procedures must be 
established to minimize the dangers from floods, an excess of water within a 
relatively short period of time. 

A stream is a constantly moving body of water providing drainage for all 
parts of the watershed. Storm runoff, from roofs, lawns, streets, farms and 
other surfaces, eventually finds its way to a natural drainage ditch or man- 
made storm sewer and then into a stream. Past experience has repeatedly 
shown that isolated sections of a stream cannot be altered without regard for 
adjacent reaches of the stream. A drainage basin is a closely correlated 
natural unit and should be studied and treated as a unit. The land along any 
stream may be divided into three zones:(2) 


1) The normal stream channel. 
2) The flood plain, which is usually dry but occasionally flooded. 
3) The high land, which is never flooded. 


In general, the stream channel is encroached upon less frequently than the 
flood plain. Many of the stream channel encroachments are bridges neces- 
sary to provide continuous and uninterrupted transportation. These bridges, 
if properly designed from a hydraulic viewpoint, nullify, to a considerable 
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extent, any adverse effect on flood flows. Most of the remaining stream 
channel encroachments occur along banks of streams. Structures, including 


physical obstacles involved in constructing structures adjacent to streams 
have deterred a greater number of encroachments on streams. The flood 
plain, on the other hand, does not have these obstacles unless it is swampy. 
Since the flood plain is usually dry, sometimes for years, there is consider- 
able encroachment on this area. The damage to these encroachments during 
times of flood account for a major portion of the total flood damages. There- 
fore, if future encroachment on the flood plain is not controlled, flood dam- 
ages will increase rapidly and the cost of providing protection will also in- 
crease. 

In many places, encroachment in the flood plain has been encouraged by a 
legal doctrine which permits man to invade not only the flood plain but even 
the stream channel itself, despite recurring and increasing flood losses. 
While some of this encroachment has been made in ignorance of the flood 
hazard, most of it has resulted from attempting to utilize the resources of 
the terrain, soil and water.(3) Farmers found easy and rich soil to till, 
manufacturers found cheap land, convenient transportation and power, rail- 
road companies found easy grades on which to lay their tracks. Following 
the usual tendency, businesses and people moved into the flood plain area to 
be near the center of activity. Largely through inertia, since some activities 
which are now in the flood plain are there uneconomically, many of these ac- 
tivities have remained in the area. 

Regulation of flood plain use has not been extensive because of: 


1) Lack of sufficient flood data. 

2) Lack of land use studies in flood plains. 

3) Objections of flood plain property owners and users to any measures 
limiting land use. 

4) Lack of legal precedent to prevent or restrict land use in the flood 

plain. 


Hydrologic and Hydraulic Aspects 


The engineer, administrator and planner, in considering flood plain zoning, 
is interested principally in the following hydrologic and hydraulic information 
and data: 


1) Flood seasons. 
2) Duration of flooding. 

3) Flood stages. 

4) Flood velocities. 

5) Extent of flooding. 

6) Minimum and/or economical channel widths. 
7) Frequency and magnitude of the greatest flood. 
8) Potential greatest flood. 


Knowledge of the flood seasons and probable duration of flooding will 
govern to a large extent what uses may be permitted on a flood plain. Flood 
stages, flood velocities, and the extent of flooding will assist in delineating 
the flood plain and in determining the ground floor elevations and the type of 


fills, which extend into the existing stream channel frequently restrict the E 
stream to the point where higher flood stages result. The economic and he 
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buildings to be permitted in the flood plain. The extent of flooding and flood 
stages are generally required for the zoning plan. The last three items are 
required to arrive at an economical appraisal of the value of zoning a flood 
plain. 

Nearly all floods are caused by excessive rainfall, which may be augment- 
ed by water from melting snow in some areas. Heavy rains result from a 
combination of meteorological factors which may occur at any time in an 
area. Meteorologic and climatic factors vary from area to area and deter- 
mine the type of storms to be expected, the time of occurrence of such 
storms, and the probable severity of each type of storm in each area. The 
two principal types of storms which produce enough rainfall to cause floods 
are the thunderstorm and the stagnant-cold-front storm. Combinations or 
variations of these two types of storms, in addition to simple lifting of warm, 
moist air masses over mountains, may result in flood-producing rains. 

The physical characteristics of drainage basins vary considerably from 
one basin to another, The upper portions of many drainage basins have rela- 
tively steep gradients. In these areas floods will have a rapidly rising and 
falling flood hydrograph. The flood plains are narrow with little storage and 
velocities are high. Large quantities of debris and sediment may be trans- 
ported for considerable distances. As the streams approach the sea, the 
gradients decrease and the flood hydrograph has a long period of rise, a 
broad peak, and a longer period of fall. Considerable flood plain storage is 
usually present; inundation of the flood plains may be protracted and sedi- 
ment and debris transported from the upper reaches may be deposited on the 
flood plains. The rivers tend to meander and frequently change their course. 

It is impossible to separate the flood plain from the river itself in order 
to consider its hydrologic and hydraulic aspects. All streams have flood 
plains along their entire length, although their width may vary from zero, in 
reaches which have perpendicular side slopes, to many miles in nearly flat 
plains, The behavior of the river in flood stages is determined not only by 
the physical properties of the normal channel, but also by the physical 
properties of the flood plain, 

Protection against the maximum flood which has occurred or is likely to 
occur is not considered feasible or desirable in some areas, particularly in 
small drainage basins.(4) Instead, protection is provided against a flood of a 
given frequency. The question of whether or not partial protection against 
floods is justifiable is beyond the scope of this paper. The fact remains, 
however, that some flood control projects are being designed and constructed 
for partial protection only. Therefore, a very important consideration in any 
analysis of a flood situation is the frequency with which floods of certain 
magnitudes are likely to occur. The frequency of floods is usually estimated 
by statistical analysis of the available data. The reliability of the analysis 
depends on the length of flood record, with respect to time, the accuracy of 
the data used, and the extent to which the flood frequency curve is extrapo- 
lated. The frequency and magnitude of floods, together with an analysis of 
the damages experienced during each flood and the estimated cost of protec- 
tion for each flood, determine to a large degree the flood control design. 


Economic Considerations 


Until the advent of the motor age, it was economically advantageous for 
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many industries to be located in the flood plain.(5) Convenient water trans- 
portation, usable processing water, a ready and inexpensive method of dis- 
posing of wastes, and the low cost of land more than compensated for the 7 
damages resulting from occasional floods. Some of these reasons no longer oe 
have an economic advantage. Today, many industries are seeking locations, a 
on or near major highways, which have a dependable water supply and attrac- 
tive surroundings. The problem of getting to work is no longer a major con- 
sideration since, with the widespread use of the automobile and good roads, 
employees are not only able and willing to commute to work but prefer the 
wider choice of home sites this mobility gives them. However, some indus- 
tries and utilities do have a valid reason for remaining in the flood plain. 

Flood plains have generally been considered as the ideal location for sew- 
age treatment plants and wells for water supply. These facilities are neces- 
sary for the continued existence of modern cities. Those industries which ps 
depend on water transportation for their supply of raw materials, and as a ee 
means to ship the finished product, regard the flood plain as the only feasible 
location for their operations. In these instances, structural adjustments and 
emergency measures would alleviate the extent of flood damages. On the 
other hand, telephone exchanges, hospitals, penitentiaries and essential in- 
dustries definitely should not be placed in the path of any conceivable flood 
waters. a 

The rich alluvial soil in some flood plains is considered by some persons x 
to mean that the flood plain was ordained to have an agricultural role. In fo 
other instances, the swampy conditions of the flood plain has led to the filling : 
in or draining of these swamps to reclaim the land for industrial and residen- ge j 
tial use, a role which is encouraged by the growth of cities and increasing a 
land values. 2 

In some areas, the extensive use of the flood plain makes it more economi- 2 { 
cal to provide protection to that area than to evacuate it. The protection, a 
however, is only against a flood of a certain magnitude. There is no assur- 
ance that a greater flood will not occur causing considerable damage because 
the feeling of security generated by the flood control works reduces the pre- 
cautionary measures which were taken before the flood control works were 
installed. This feeling of security is likely to occur regardless of whether or 
not partial or full protection against floods is provided. The existence of the - 
protective works, in addition, usually encourages further use of the flood - 
plain. 

For most towns and cities, which have local flood problems, the economic 
value of the flood plain does not warrant extensive and expensive control 
structures. Frequently, it is more economical to develop a stream valley in- 
to a park rather than to enclose the stream in a limited-capacity culvert. (6) 
The flood plains adjacent to the stream, utilized for parks, playground areas, : 
golf courses and other recreational activities either public or private enhance a 

4 


the beauty of the city and contribute materially to the welfare of its citizens. 
The damage experienced during flood periods is minor and usually can be 
taken care of with the regular maintenance force. However, park buildings 
and facilities subject to extensive damage should be placed in areas which 
will not be flooded. On most small streams, the flood plain is inundated less = 
than one per cent of the time per year. It is doubtful that this minor inunda- : 7 
tion would affect the lawn areas to any extent if no oil or harmful chemicals * 
were present in the water. 

In determining the feasibility of flood plain zoning from an economic point 
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of view, the following items should be considered: 


1) Existing and potential flood damages. 

2) Cost of protection. 

3) Benefits from protection. 

4) Cost of evacuation. 

5) Benefits from evacuation. 

6) Allowable land use. 

7) Community benefits. 

8) Structural adjustments of existing encroachments. 

9) Emergency measures such as flood warning services, rescue squads 
and elevating goods and machinery above possible flood heights. 


Social Implications 


In ancient times, the most favored locations for cities were those where 
conditions were suitable for agriculture and for transportation. The rich 
alluvial soil in the flood plains with the convenience of nearby water trans- 
portation satisfied these requirements. These cities have expanded, in most 
instances, to higher surrounding areas above possible flood heights. How- 
ever, the flood plain location of the original portions of many cities has been 
and still is reflected in many of the social problems facing cities today. 

Periodic flooding of residential areas in flood plains fosters the develop- 
ment of blight. The unrealistic economic values, the high cost of providing 
services, the deterioration of the areas due to unwholesome topographic and 
physical surroundings—all contribute to the conditions existing in many de- 
veloped flood plains. The existence of some of these factors is evident in the 
difficulty in evacuating the area or in providing adequate protection to the 
area. This protection usually requires extensive and costly flood control 
works, the proportioned share of which is usually beyond the financial ability 
of the people directly concerned. Therefore, the community, at large, is 
faced with three alternatives: 


1) Do nothing to the area and continue to pay for the high cost of 
services required and to absorb the cost of flood damages. 

2) Evacuate the area and retain it as public lands to be used for 
suitable functions. 

3) Provide adequate flood control works to protect the area. 


Protection of the Flood Plain 


The extent to which flood plains are zoned determines the degree of flood 
damage prevention which will result. The minimum amount of regulation in a 
flood plain prohibits land development which causes a construction of the 
natural stream channel. This is sometimes called the “channel-capacity”(3) 
regulation since the degree of restriction depends on the capacity of the chan- 
nel to discharge a flood of a certain frequency. The design frequency depends 
either on an economical analysis of the flood damages and benefits from pro- 
tection or on the degree of protection desired. Naturally, a flood greater than 
the design flood will inundate the adjacent flood plain. The maximum regula- 
tion of the flood plain limits the use of the flood plain to purposes which do 
not suffer from the action of a flood or to purposes which, if damaged, cause 
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no serious inconvenience to persons or interests located outside the flood 
plain. This regulation is often referred to as the “land-use”(3) regulation, 
which makes possible complete prevention of flood damage when this is con- 
sidered necessary or desirable. 

Most cf the existing flood plain ordinances in existence today are based on 
the “channel-capacity” regulation. The reasons for this are: (1) it is easier 
to determine, (2) its enforcement is simple, (3) legal precedents have been 
established, (4) the maximum amount of land is available to owners of flood 
plain property for use and (5) it directly affects only riparian owners. “Land- 
use” regulations are becoming more prevalent despite the restriction placed 
on the land. Often, however, complete restriction is considered unnecessary 
and some use of the flood plain is permitted provided certain safeguards are 
met. 

Inherent in the establishment of any area as unsafe because of possible 
flooding is the implication that the areas which are not thus delineated are 
secure from possible damage from floods. The degree to which the danger of 
this implication exists depends on the criteria used in establishing the limits 
of flooding. Flooded areas delineated by past floods may or may not include 
the area which could be potentially flooded by a future flood. It is very pos- 
sible that, as drainage basin characteristics and channel conditions are al- 
tered, a greater flood height than had been previously experienced may be 
expected. In fact, it is always possible to have a greater flood than any pre- 
vious flood regardless of whether or not there is any topographic change 
within the basin. 
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Discussion of 
“SYNTESSIS OF RAINFALL-INTENSITY-FREQUENCY REGIMES” 


by D. M. Hershfield, L. L. Weiss, and W. T. Wilson 
(Proc. Paper 744) 


D. M. HERSHFIELD,! L. L. WEISS,2 AND W. T. WILSON,? A.M. ASCE.— 
We appreciate Mr. Sammons’ complimentary reference to multi-graphical 
correlation and his suggested improvements. In the evolution of figure 1 the 
primary requirement was to estimate 2-year, 1-hour rainfall for localities 
lacking even daily rainfall data, The daily rainfall parameter was added to 
the rest of the diagram for application to localities having daily data. Before 
applying the diagram the writers tested it with independent data from many 
parts of the world, and noted that it had rather broad applicability with toler- 
able error. 

In arctic and subarctic regions this diagram was not as successful as an- 
other, shown here in figure 3. The large proportion of snow in the mean an- 
nual precipitation of arctic regions does not contribute much to a rainfall in- 
dex, and days-of-precipitation is a difficult observation to make in snow 
regions. The two parameters, mean annual precipitation and number of days 
of precipitation per year, combine to give a measure of rain per rainy day. 
In a large measure the daily precipitation parameter which the writers added 
to the diagram in figure 1 is redundant to these two parameters. 

It has been observed that over much of the world the average ratio of 2- 
year, 1-hour rainfall to 2-year, 24-hour rainfall is about 0.40. This ratio 
may be 0.20 or lower in maritime regions having a low incidence of thunder- 
storms, while in continental regions where there is a high incidence of con- 
vective activity the ratio may exceed 0.60. Figure 3 shows the relationship 
between annual extremes of hourly and daily rainfall as a function of thunder- 
storm incidence. The diagram is based on data from stations in Alaska, 
Canada and northern United States. 

To apply this diagram to localities having more than 25 thunderstorms per 
year, the diagram would have to be extended. The shape of the diagram may 
vary slightly from region to region. We believe that this type of diagram, 
where the required data are available, and derived for the pertinent climatic 
region will be at least as good, simpler, and possibly superior to the one the 
writers published in the paper. 

While the foregoing discussion is not specifically a reply to Sammons’ 
discussion, the writers’ paper was presented as a progress report, and it is 
felt that recent progress is pertinent. 

Coming, now, to some of Sammons’ points: The writers have found that 
1.13 is a better average factor for converting annual series to partial- 


1, Weather Bureau, U. S. Dept. of Commerce, Washington, D. C. 

2. Weather Bureau, U. S. Dept. of Commerce, Washington, D. C. 

3. Meteorologist, Weather Bureau, U. S. Dept. of Commerce, Washington, 
D. 
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duration series 2-year, 1-hour values; the 1.09 factor was based on a rela- 
tively small sample of data and, of course, even at the 2-year level there is 
some subjectivity in most methods of fitting curves to partial-duration data. 

In deriving the diagram that Sammons discussed, the writers were not 
concerned with the distribution of any of the parameters; they were dealing 
with a statistic, essentially the mean. The relationship of the 2-year value 
to the mean varies trivially among the several methods for analyzing annual 
series. In figure 1, no particular distributions were explicitly assumed. The 
scales used were largely a matter of convenience, with their inter-relation- 
ships fixed, of course, by the relationships of the data themselves. The em- 
pirical relationships of the figure are not dependent on subsequent discoveries 
as to the nature of the distributions of the variables used. 

It is not clear how Sammons obtained “Weather Bureau ratios” for the last 
table of his discussion. In the paper the writers explicitly suggested 5-year 
and 10-year ratios to the 2-year value of 1.4 and 1.7 as provisional average 
values in lieu of data for computing them. Weather Bureau Technical Paper 
No. 251 shows these ratios for about 200 United States stations. There is 
considerable variation among stations, but it is uncertain how much of this 
variation can be ascribed to physical and climatic causes, and how much is 
sampling error. 

Sammons raised a question about the writers’ stated record length averag- 
ing 40 years with respect to a base period. To avoid secular trends it isa 
common practice to restrict the data to a common period. This means elimi- 
nating data from outside this period, and stations whose records do not occupy 
all of the period. This period is sometimes made rather short to avoid ex- 
cluding too many stations. The writers felt that most of their records were 
long enough and for nearly enough the same period of time to be included. 
They felt that they would lose more by restricting the volume of data than by 
trying to get a homogeneous period. At present the writers are testing this 
issue, and at the same time are experimenting with the effect of duration on 
the frequency distribution. The writers might possibly find that the frequency 
distribution of daily amounts represents a different population than of hourly 
amounts, in some parts of the country. 

The writers might expect to find for some localities that hurricane rainfall 
defines the regime of daily extremes, while thunderstorm rainfall defines the 
hourly regime. In other localities they may find that the same storms fre- 
quently provide both the hourly and daily extremes. The writers want to find 
out under what conditions the relationship of daily to hourly extremes is a 
climatic index, and under what conditions it is a physical relationship. 


1. U. S. Weather Bureau, Department of Commerce, Weather Bureau Techni- 
cal Paper No. 25, Rainfall Intensity-Duration-Frequency Curves for Select- 
ed Stations in the United States, Alaska, Hawaiian Islands and Puerto Rico, 

Washington, December 1955, 
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Discussion of 
“GROUND WATER PHENOMENA RELATED TO BASIN RECHARGE” 


by Paul Baumann 
(Proc. Paper 806) 


PAUL BAUMANN,! M. ASCE.—Mr. Lee correctly evaluates the writer’s 
equation 1, although expressing the acceptance rate in c.f.s. per square foot 
instead of per wetted acre would have been preferable for the following rea- 
sons: 

In the study of two-dimensional ground water movements a strip one foot 
in width is normally considered. Hence, the spreading grounds and the 
stream channels, the width of which rarely is constant, may be divided up into 
as many such strips as are required to reasonably reflect the effect of the 
variable width on the total acceptance rate. 

In constructing the flow net Mr. Lee apparently was unduly influenced by 
the classical textbook example of percolation from a reservoir per unit length 
of a floating dam of infinite length on a saturated, porous foundation, with the 
less featured case of flow from a spreading ground under a river levee. The 
two cases are fundamentally different. The get-away of percolation from a 
saturated stream channel is approximately at right angle thereto, namely, in 
the direction of residual pore storage and of a control, as shown in the writ- 
er’s Figure 3-a. In the longitudinal direction of the saturated stream channel 
there is no residual pore storage and no control and therefore no get-away for 
additional sub-surface flow. Hence, in the absence of lateral escape, all per- 
colation under a floating dam from a reservoir must rise to the surface, 
whereas, percolation under a levee from an adjacent spreading ground must 
not so long as there is a get-away, by means of residual pore storage and a 
control. However, return flow to the river from the spreading ground may 
take place, depending upon the gradient between the water surfaces of spread- 
ing basin and stream channel as stated by the writer. 

The flow net shown by Mr. Lee does not consistently satisfy the basic re- 
quirements therefor. These are: (1) the vertical drop 4h between the 
equipotentials along any phreatic line must be constant; (2) the flow 4q ee 
through each of the flow paths must be constant; (3) the flow 4q must be . 
equal for all of the flow paths; (4) stream lines and equipotential lines must 
form quasi squares between themselves except near the boundaries where 
distortion is unavoidable. Not one of these conditions is reasonably satisfied 
by Mr. Lee’s flow net. It therefore is unbalanced and cannot be accepted. 

Fig. 11 shows a balanced flow net for the rather extreme condition similar 
to that chosen by Mr. Lee. The gradient between basin and stream channel 
water surfaces used in Fig. 11 is 2:1, as compared to the 5:1 gradient re- 
ferred to by the writer as the limiting value for insignificant return flow to a 
the stream channel, By assuming a number of phreatic lines through the levee it was 


1. Asst. Chi. Engr., Los Angeles County Flood Control Dist., Los Angeles, 
Calif. 
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possible to find the one for which the return flow is a maximum. For this 
condition the return flow was found to be 7.5% of the total flow toward the 
stream channel of percolation from the spreading basin. Hence, for a gradi- 
ent of 5:1 the return flow would indeed be insignificant as stated by the writer. 
In arriving at the return flow to the stream channel, Mr. Lee computes the 
flow through two vertical sections on the premise that the respective gradi- 
ents are equal to the difference in water surface elevation divided by the 
horizontal distance to the section. This is assuming horizontal flow through- 
out the entire reach in accordance with the Dupuit theorem. Under “Introduc- 
tion” the writer pointed out that... “the Dupuit theorem is applicable except 
in special cases” ... Subsequently, under the heading “Strip Basin” he stat- 
ed, “It is apparent that the Dupuit theorem is applicable with fair approxima- 
tion except near the spreading basin and near the control.” There the flow 
conditions cannot be ascertained without a flow net. The total flow computed 
from the flow net checks that in accordance with the writer’s eq. (1) within 
15%, which is a close check considering the extremity of the example. 

For the water year 1940-41 (October 1 to September 30) the rainfall for 
Los Angeles County was in excess of 200% of normal. Sustained runoff from 
the San Gabriel Mountains took place throughout the winter months and re- 
charge of ground water basins through percolation from stream channels and 
spreading grounds caused mounds to rise to the surface in the “forebay” of 
the Central Coastal Basin by the Spring of 1941. Fig. 12 shows the ground 
water contours and ground surface contours in the vicinity of the San Gabriel 
and the Rio Hondo Spreading Grounds of the Los Angeles County Flood Con- 
trol District as of May 6, 1941. These spreading grounds are located in the 
forebay. At stream channels the ground surface contours show average ele- 
vations between banks. Low water channel beds would be two to three feet 
lower. Of special interest is the merging of the mound beneath the San 
Gabriel Spreading Grounds with the channel mound, particularly in view of 
the fact that the get-away of percolation the east of the San Gabriel stream 
channel was more restricted than the get-away to the west of the Rio Hondo 
channel. Throughout this period there was no evidence of rising water in the 
saturated San Gabriel stream channel adjacent to the levee separating it from 
the spreading grounds, although the gradient from the spreading mound to- 
ward the stream channel is conspicuous in Fig. 12. Certainly there is no 
evidence of blocking by the channel mound of percolation from the spreading 
grounds. Of further interest is the curvature of the streamline trajectories 
normal to the ground water contours as referred to by the writer in connec- 
tion with the channel flow net. A typical trajectory is indicated by a dash line 
in Fig. 12. It is hoped that this will also serve to clarify the question in re- 
gard to curvature of flow nets raised by Prof. Brooks. 

Mr. Lee questions the similitude between the model test results and those 
of a prototype expressed by the writer’s eq. (1). The model test which was 
relatively free from capillary flow!3 showed a maximum discharge parame- 
ter q/K = 0.097 whereas eq. (1) for W = 2.4" (not 1.2""as used by Mr. Lee) and 
D = 12" gives q/K = 0.092, a difference of 5-1/2%. 

Prof. Brooks undertakes to correlate the evaluation of the writer’s eq. (1) 
and his eq. (2) for the basic parabola. This is in keeping with fundamentally 
sound reasoning. However, in developing his eq. (2) so as to appear in terms 
of D/W Prof. Brooks apparently overlooked the fact that the depth D 


13. bid, p. 1059. 
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in his Fig. A can never be equal to D in the writer’s Fig. 3a. As previous- 
ly14 shown the phreatic line between the spreading ground and the control 
can be approximated by a double parabola. D in Prof. Brooks’ Fig. A could 
therefore be considered as the depth below the point (P) of counter flexure, 
common to both, the lower and upper parabolas. Hence, for identification, D 
and L in Prof. Brooks’ Fig. A should be changed to d and 1. Prof. Brooks’ 
example for arbitrary values of D/W = 2 and L/D = 10 and the apparent dis- 
crepancy between the results based on the writer’s eq. (1) and on the basic 
parabola are therefore misleading because of the fact that for potential flow 
between the spreading ground and the control the ratio q/K must be common 
to both the spreading ground and the lower basic parabola. As explained by 
the writer, eq. (1) was developed from flow nets which necessarily terminat- 
ed in basic parabolas. Hence, all corresponding values are found as the in- 
tersection of the 1/d lines with the (1 - e~4/W) curve in the writer’s Fig. 13. : 
It is interesting to note that the 1/d = 3/4 line becomes tangent to the (1 - ha . 
e-d/W) curve as q/K approaches zero. In developing e-4/W into a series it e 
also follows that as 


d/jw — 0, g/K — d/W 


Actually, as may be seen from Fig. 13 the q/K ratio closely checks the 
d/W ratio if the latter is less than 0.1. 

In designing a spreading ground it is necessary to ascertain the economic 
limitations. Ordinarily, the depth D = 2d to the impervious stratum is known. 
Hence, by means of the writer’s eq. (1) and Prof. Brooks’ revised eq. (2) cor- 
responding values of q/K may be determined as shown in Fig. 14. 

For the case quoted by Mr. Lee, namely, D ~ 100 feet, it follows from 
Fig. 14 that for D/W = 0.10 and therefore W ~ 1000 feet, the infiltration 
parameter per square foot q/K = 0.05 or q = 3.35 c.f.s. per wetted acre for a 
Meinzer Unit = 1000 as used by Mr. Lee. 

Since Figs. 13 and 14 are correlated, the absolute value of 1 (which is half 
the distance between the edge of the spreading ground or stream channel and 
the control) may be ascertained from Fig. 13. 

As previously stated, extraction of ground water through pumping is to be 
expected where recharge through spreading is performed. Coordination of 
the two operations should include the locating of water wells as closely as 
possible to the theoretical control. However, the curvature in plan of the 
flow paths, previously discussed and indicated in Fig. 12, mitigates the need 
for precise location of the control, based on two-dimensional flow. For, 
prototype ground water movement, contrary to that in models, is seldom 
strictly two-dimensional but can, within limits, be approximated as such. 

Pumping in the spreading area proper, although effective in regard to in- 
filtration, may not be desirable because of a short path of percolation and ; 
possible contamination especially in shallow ground water basins. > 

The dash line near the division of flow or axis of symmetry in the writer’s 
Figs. 3a and 3b as questioned by Prof. Brooks is not intended to be a bound- 
ary. If it were, the equipotential lines would not extend across it. Hence, 
flow takes place along the vertical division line and the dash line merely 
meant to indicate the effective limit of the flow Aq in the flow path next to 
the boundary. Admittedly, the effective limit of flow could be modified to 
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14. Ibid, p. 1042-1045. 
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advantage as shown in Fig. 15, the symmetrical flow net for the free water 
surface of stream channel and spreading ground at the same potential. 

Prof. Knapp, in presenting his interesting and constructive discussion has, 
indeed, found a flaw in the crystal ball in that it prophesied ultimate condi- 
dions when, in fact, they were just “prior to ultimate.” Of particular signifi- 
cance is Prof. Knapp’s reference to interface slopes under the influence of 
pumping. As may be judged from the writer’s Fig. (4) the draw-down gradi- 
ents adjacent to the test site are of the order of 1/2 of 1%. Hence, in using 
the 40 to 1 ratio as correctly suggested by Prof. Knapp, the slope of the in- 
terface would be of the order of 20%. This phenomenon is somewhat startling 
at first glance. However, normally sea water spread along the bottom of an 
aquifer is “sucked” into wells at moderate rates as compared to fresh water 
because of the lowest well perforations being at such a distance above the 
bottom of the aquifer as to require the interface to assume its maximum 
slope in rising to the well. 

This clearly brings out the difference between the saline wave and sea 
water intrusion proper, questioned by Prof. Hotes and Dr. Harder, in that 
during the passing of the saline wave the lowest perforations of affected wells 
may be temporarily blocked, whereas no measure of protection will save a 
well from abandonment if sea water intrusion is permitted to continue un- 
checked. Also, those wells adversely affected by the saline wave might be 
temporarily used to extract sea water cut off from the ocean by the barrier. 
In more timely terms the dissipation of the cut off part of sea water intrusion 
may perhaps be compared with containing the cut off part of an invasion army 
as compared to containing the whole invasion army. 

It was the writer’s privilege to observe the ingeniously conceived model, 
shown in Prof. Hotes’ Fig. 1, at the Richmond Field Station of the University 
of California. As is clearly brought out in Prof. Hotes’ discussion, this 
model lent itself to simulate well-nigh all probable prototype conditions and 
to evaluate results quantitatively within reasonable limits of accuracy. As in 
all cases where ground water studies are concerned, the time scale would 
have to be used on the usual assumption of homogeneity of the aquifer as well 
as uniformity in magnitude. 5; 

Prof. Hotes’ eq. (1) which was developed by Dr. Harder and which is an 
adaptation of an example of flow through an earth prism analyzed by Dr. 
Muskat, is based on the assumption of a constant potential or equipotential 
along a vertical line through the injection well; a constant potential of equal 
magnitude along the horizontal sole and a variable potential along the verti- 
cal line through the seaward control. Hence, the fresh water flow passes 
through the control between the bottom of the clay cap and the sole of the 
aquifer with no flow along the latter. Exactly the same result as obtained by 
Dr. Harder can be obtained by applying the Dupuit theorem to flow between the 
clay cap and the interface. The shape of the interface is such as to have bal- 
anced pressure between fresh water and sea water at all points. Its shape is 
expressed by 


---+-+-- (2) in which, as 


shown in the writer’s Fig. 16 
z = Ordinate of interface above sole (x - Axis) 
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M = Magnitude of aquifer 

x = Abscissa relative to center of recharge well 
q = Rate of flow per unit width of aquifer 

K = Permeability coefficient 

s = Specific Gravity of sea water 

1 = Specific Gravity of fresh water. 


For any stable condition of the interface, the rate of flow q must be constant 
and a maximum. Hence, because the radical in the writer’s eq. (2) must be 
real, the two terms under the radical for qmax must be equal. Thus for 
x=L 


and q=1 (s-1) KM’ = 1 (s-1) MT --- 
2 L 


which checks Prof. Hotes’ eq. (1). 
It now follows that for the above condition to obtain, the product qL of the two 
variables must remain constant. Hence, as q increases, L decreases propor- 
tionately and vice versa as is evident from Prof. Hotes’ Fig. 3. From the 
writer’s eq. (2) it would follow that for any stable condition of the interface 
and therefore any corresponding product qL, the height z of the interface at 
the control is equal to M. As this would leave no aperture for the discharge 
of fresh water, it is evident that eq. (2), based on the Dupuit theorem, is an 
approximation of the true solution. To satisfy z = o for x = o the sign before 
the radical must be negative. 

From the writer’s eq. (2) it follows that the slope of the interface 


dx K(s-1) - 


K(s-1) 


and that for x = L the slope for any corresponding q and L values is 90°. 
While this answer again reflects the afore-mentioned approximation, it ap- 
pears to be reasonable if we remember Prof. Knapp’s 40 to 1 ratio between 
the fresh water gradient and the slope of the interface. In accordance with 
this ratio, any fresh water gradient in excess of 2-1/4 degrees would call for 
a slope of the interface at the control of 90° which, of course, cannot be ex- 
ceeded. 

If the interface were treated as an inverted, basic parabola and if the 
boundary condition z = M - q/K at x = L were imposed to bring the Dupuit 
solution in accord with the Kozeny!5 solution, then the writer’s eq. (2) would 
lead to 


----@) 


15. Ibid, p. 1040. 
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and q=- KL + ) + 


s-1) (s-1)2 


While no authenticity is claimed for eqs. (6) and (7), it is encouraging to 
find that the boundary conditions L=o °.q= TandL > ©.*. q > Oare 
satisfied. Also that q is of the same order as that from eq. (4). 


Putting eq. (6) inthe form q/K= ( G-1)*-1 ) 
2 


8-1 
L 
and expanding the radical into a binomial results in 


q/K = M*(s-1) - 1 M* (s-1)° 


from which it follows that for any finite discharge face the rate of flow q 
must be smaller than in accordance with Prof. Hotes’ eq. (1). Hence, the 
latter equation should be corrected to read 


q = K (s-1) M2 
2nL 


in which n > 1 and increases as the ratio M increases. 
L 

Obviously, if fresh water from the recharge well were discharged from a 
pressure level, h feet above M.S.L. at the recharge well into a fresh water 
lake with the water surface at M.S.L. and if the interface were fixed, yo at 
the control would be close to zero and the pressure gradient or piezometric 
line would be as indicated in Fig. 16. If now it were possible to suddenly con- 
vert the back pressure due to fresh water to that due to sea water, the dis- 
charge would stop just as suddenly. But the pressure head h would now exert 
itself and would cause fresh water pressure to exceed sea water pressure ex- 
cept at the interface where the two pressures would be balanced. Hence, for 
a liquid interface the blocking of flow at the base of the control is offset, at 
least in part, by the pressure head y, on top of the control. 

This moves into focus the danger of “blow-out” at the control which grows 
with gq and diminishes with the increase in the distance L between the toe of 
the salt water wedge and the control. Furthermore, this danger is greatest 
during periods of extreme low tide, a feature particularly significant at loca- 
tions subject to large tidal ranges, Fortunately, the extreme tidal range 
along the coast of California is moderate, 

It is possible that offset in the model, shown in Prof. Hotes’ Fig. 1, be- 
tween the back of the sea water weir box and the control face gave rise to the 
trapping of air split out of the fresh water and that the q measurements were 
thereby affected to some extent. 

Of considerable interest would have been the determination of the equipo- 
tential lines and therewith the flow distribution between the interface and the 
“clay cap” of the model. This is just one of many facets meriting additional 
research. 
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The discharge through the control (16) to the ocean of fresh water could be 
more readily comprehended if the interface, as a basic parabola, were ex- 
tended so as to terminate at the horizontal discharge face, q/2K in width and 
flush with the bottom of the clay cap. The discharge of fresh water would 
then be vertically upward and would thereby satisfy the physical premise of 
buoyancy of fresh water relative to sea water. 

Of special interest are the 10 stages shown in Prof. Hotes’ Fig. 4 and par- 
ticularly those pertaining to high injection rates. It must be assumed that the 
condition as shown at the bottom of the right-hand pictures would represent 
close to the ultimate condition of the sea water wedge cut off by the barrier, 

Dr. Harder computes the contamination of 240,000 acre feet of water to a 
concentration of 250 parts per million of chloride ions (Cl~) because of dif- 
fusion of that part of the sea water wedge which had intruded landward from 
the injection wells. If sea water in intruding a fresh water aquifer mixed as 
readily with the fresh water through diffusion as Dr. Harder seems to as- 
sume, the correction of sea water intrusion would indeed be a next to hope- 
less case unless it were always possible to locate the injection wells land- 
ward from the extreme tip of a saline wedge as suggested by both Prof. Hotes 
and Dr. Harder. At Manhattan Beach this procedure was carefully considered 
but had to be abandoned because of three principal features. First, the land- 
ward pressure gradient of the aquifer would have made it necessary to create 
a barrier above mean sea level of approximately twice the height of that at 
the location where the barrier was finally located; second, the right of way 
necessary for the barrier including pipe lines, pressure regulators, chlori- 
nators, etc., would have been far too expensive to be justified in the light of 
an experimental period of operation; and third, facilities for the disposal of 
water produced during development of the recharge wells were completely 
lacking. However, the writer still maintains that the effect of the saline 
wave, under similar conditions as those at Manhattan Beach, is negligible as 
compared to that of the sea water intrusion proper. It is only necessary to 
visualize the irregularities normally existing along the bottom of an aquifer 
in the form of mounds and depressions to realize the facilities for ready ab- 
sorption of the sea water cut off from the original wedge by the injection and 
forced landward. The average depth of sea water so absorbed may readily 
be such as to require relatively short landward distance of spread. 

In referring to the writer’s Fig. 8, Prof. Hotes concludes that the saline 
wave is worse than indicated by the writer. Apparently, Prof. Hotes over- 
looked the significance of the writer’s Fig. 7 which shows more tellingly than 
Fig. 8 what happened when the cut off sea water wedge was forced landward 
by the barrier wells. For example, it shows that the salinity in well G-8, 
1180 feet from the line of recharge wells dropped from about 5000 p.p.m. of 
Cl” concentration to approximately 1000 p.p.m. in 16 months of barrier 
operation. Since October of 1955 this particular well has been immersed in 
fresh water of less than 250 p.p.m. of Cl“. Prior to commencement of the 
recharge test, sea water displaced fresh water throughout the entire depth of 
the aquifer at a rate of some 500 feet a year. Hence, in the three years and 
two months since the test was started unadulterated sea water would have 
overtaken well G-8 and would have advanced some 400 feet beyond it. The 
writer’s “optimism,” therefore, does not seem to be entirely unjustified. 


16. Ibid, p. 1045. 
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Dr. Harder and Prof. Hotes’ points are however quite well taken in that 
sea water that had intruded into a shallow aquifer for a considerable distance 
landward and had completely displaced the fresh water for such a distance 
could not be checked by means of injection wells without seriously affecting 
close-by inland wells by the displacement of the landward part of the saline 
wedge. On the other hand, wells close to the barrier would, under conditions 
similar to those at Manhattan Beach, be close to the ocean as well and there- 
fore subject to permanent and irreparable instead of temporary damage. 

The model tests conducted by the University of California at Richmond 
have certainly paid dividends far beyond the return that could reasonably 
have been expected. 

The writer is deeply appreciative of the spirited discussions through 
which great value was added to his paper. He only regrets his inability to 
agree with all of the thoughts therein presented. 
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Discussion! of 
“EXTENDING STREAMFLOW DATA” 


by W. B. Langbein and C. H. Hardison 
(Proc. Paper 826) 


F. A. JOHNSON, 2 A.M. ASCE.—A graphical method is presented for esti- 
mating the relationship between monthly discharges of nearby streams having 
similar characteristics. The relation curve then is to be used to extend the 
shorter record on basis of the other, or long-term record; and statistical 
measures are suggested as an indication of the expected accuracy of the 
estimates. 

Several experiments made by the writer with pairs of long-term records 
led to fairly reliable extensions in some cases, but not in others. The statis- 
tics, such as standard error of the relation curve and the coefficient of cor- 
relation, did not invariably give a dependable indication of the accuracy. 

It was apparent from these tests that in some regions the areal pattern of 
runoff was relatively stable from year to year, and over long periods. In 
other regions, however, the distribution varied greatly from year to year, or 
between irregular periods of several years. These instabilities or changes 
in pattern of runoff were noticed in a number of regions, both of the Pacific 
and Atlantic Slope basins. They apparently were caused by variations in the 
areal distribution of precipitation. 

An average relationship between the monthly discharges of two nearby 
streams, say for 5 years of concurrent records, thus may differ considerably 
from that during another period, and lead tq systematic errors in the esti- 
mates. An obvious way to minimize such errors, as well as the random 
errors, is to establish a relation between the discharge in the subject or 
satellite basin and the average discharge fromtwo or more nearby basins on 
different sides of the one under consideration. For correlation of monthly 
discharge, however, this possibility often is limited by lack of sufficient long- 
term records having similar characteristics as to monthly distribution. 

One purpose of this discussion is to describe a method for estimating an- 
nual or long-term runoff of the satellite stream from two or more long-term 
base records, regardless of the monthly distribution of runoff. In many parts 
of the Atlantic and Pacific Slope basins this leads to fairly accurate estimates, 
which may be compared with extensions of monthly records on an annual ba- 
sis. Such comparisons are, in effect, independent tests of the accuracy of 
the extensions, if the annual estimates are based on different index records; 
and they sometimes may be used for purposes of adjustment. 

The need for such a test is demonstrated by an example of the authors’ 
correlative estimates which is summarized in Table 1 and 2, and Figure 4, 

In a number of similar experiments it was found that substantial errors of 


1. Publication authorized by the Director, U. S. Geological Survey. 
2. Regional Hydr. Engr., Geological Survey, U. S. Dept. of the Interior, 
Tacoma, Wash. 
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the kind illustrated here were not uncommon, although generally of somewhat 
lesser magnitude. 

In this example, the monthly runoff of the Wading River near Norton, Mass. 
was estimated from that of the Taunton River at State Farm, Mass. The 
drainage areas are 42.4 and 260 square miles, respectively. The center of 
the Wading River basin is about 15 miles east of that of the Taunton River 
basin. Both basins are near the Atlantic Coast in the southeastern part of 
Massachusetts. A relation curve was defined on logarithmic paper by 60 
pairs of observations recorded from October 1945 to September 1950. (Run- 
off figures in inches, as adjusted for diversions and storage, were readily 
available and were used without allowance for differences in length of months). 
The graphical method illustrated by the authors in Figure 1 was used for fit- 
ting the curve and estimating the statistical measures. The standard error 
of estimate was determined as 0.10 log unit, and the coefficient of correla- 
tion as 0.975. 

Estimates made for the period 1930 to 1945 from this curve are summa- 
rized by annual totals in Column 2 of Table 1 for comparison with the re- 
corded figures of Column 3. The runoff for two 4-year periods was under- 
estimated by about 20 percent, which is approximately the standard error as 
determined for monthly estimates. The relationship during the sample peri- 
od, 1946 to 1950, evidently was not representative of other periods between 
1930 and 1945. 

The authors pointed out that the statistical measures estimated from the 
sample are applicable in other periods only if the mean line is correctly lc- 
cated. However, it is difficult to see how this can be judged from internal 
evidence. It was pointed out also that these measures have the usual signifi- 
cance only when deviations from the mean line have a normal distribution. 
Although adjustments such as the logarithmic transformation often appear to 
bring the deviations into a normal distribution, the deviations may resemble 
random variations only if regarded without relation to time. The writer has 
noticed that large positive or negative deviations of monthly values some- 
times persist for periods of from 4 to 6 months. These consistent deviations 
sometimes result in an accumulation of error in the annual totals almost as 
large as the percentage errors of the monthly estimates. Because of such 
errors it is evident that the computed standard error of the relationship may 
be an underestimate even for the sample period. It may be grossly in error 
during other periods if there is a bias of the kind illustrated in Table 1 and 
Figure 4. 

It was noticed also that the seasonal distribution of discharge of some 
streams of the Pacific Slope is such that the monthly values tend to cluster 
in zones above and below the mean and median. Under these circumstances 
the coefficient of correlation does not have the same significance as for a 
normal distribution, and is apt to be fictitiously high. Coefficients of more 
than 0.99 in fact did not appear to be uncommon. 

The recorded annual runoff for the example of Table 1 is shown in Figure 
4. As would be expected, the annual values of runoff computed from monthly 
estimates fall approximately on an average curve through the recorded an- 
nual figures for the years used in definition of the monthly relationship. The 
degree of correlation between annual discharges of long-term records in a 
given region thus furnishes an indication of the reliability to be expected of 
monthly estimates for other streams in that region. Such comparisons 
would seem to be a useful preliminary step in extensions of streamflow data 
such as proposed by the authors. 
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It has been the writer’s experience that the possibility of large systema- 
tic and random errors can be minimized most effectively by establishing a 
relationship between the subject stream and the average of several nearby 
streams, rather than between two streams alone. Along the Pacific Coast 
correlations seem to be the most reliable where the index basins are on op- 
posite sides of the subject basin, and on a line roughly normal to the west- 
east, or average direction of the storm paths. A similar condition may exist 
in other regions. In connection with their discussion of snow-fed streams, 
the authors mentioned the possible need, at times, for use of more than one 
index of discharge to improve the correlation. It would seem that such pre- 
caution should be taken as a matter of course, since it can seldom be known 
that a group of observations in one period is closely representative of condi- 
tions in another period. 

If the monthly distribution of runoff in the several basins is similar, the 
average of several records may be used in the same way as was illustrated 
by the authors with a single index record. However, where there are sub- 
stantial differences between the altitude characteristics of the basins, and 
runoff results both from rainfall and snow melt, the monthly distribution may 
differ considerably. In that event, average indices of annual runoff may be 
used for defining relationships of annual runoff. Estimates of annual runoff 
then will serve to verify the accuracy of annual sums of monthly figures es- 
timated in the way proposed by the authors. 

Although the error due to variations in the areal distribution may be mini- 
mized by using averages, a relation line defined by only a few pairs of annual 
values obviously may not be dependable for purposes of extension. Fortunate- 
ly, in humid regions where the precipitation is substantially more than the 
normal evapotranspiration, fairly reliable estimates usually can be made by 
a computed relationship described ir the following paragraph. This has been 
used by Knox and Nordenson(1) for studies in the northeastern part of the 
country, but its wider applicability may not be generally recognized. 

Ignoring systematic and random variations in water loss, and changes in 
natural storage, the relationship between the annual precipitation on nearby 
areas can be expressed as follows: 

Ry + Ly = b(R, + 
where Ry and Rx are values of annual runoff, L, and Lx are estimated values 
of mean annual loss, and b is the ratio of precipitation on basin “y” to that on 
basin “x.” Both terms of the equation are subject to error equal to the devia- 
tion of the actual losses from the assumed average losses, plus the changes 
in natural storage. However these errors tend to be the same for nearby 
areas, and thus tend to cancel out. The ratio, b, often can be estimated 
closely from the average of only a few pairs of annual observations, even if 
their range of magnitude is small or negligible, so that a relation curve could 
not be defined graphically. The relation, Ry = b(Rx + Lx) - Ly then can be 
used to estimate values of Ry from Rx during periods of missing record. 
The results usually are affected only slightly by fairly large differences in 
the assumed values of average loss, Lx and Ly, from actual values, or by 
variations in the amount of annual loss with variations in precipitation. In- 
formation regarding natural water losses can be obtained from a number of 
sources. Some data were compiled by Rafter(2) as early as 1903, later by 
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Williams and others(3) in 1940, and by Lee(4) in 1942. 

This method of estimation is illustrated in Table 2, as a supplement to 
Table 1. (In this example the ratio, b, happens to be 1.00. The method 
seems to work just as well where there are considerable differences in the 
values of areal precipitation). The annual estimates of Table 1, as well as 
for a number of unlisted experiments, are in general substantially more ac- 
curate and consistent than totals derived from monthly correlations between 
two streams alone. The improvements were due to the use of more than one 
index record for the annual estimates. The same improvement in monthly 
estimates could be obtained by this means, where the monthly distribution is 
closely similar in all of the basins. 

In an actual extension of records, the reliability of the estimates can be 
judged only on basis of the errors for the years of definition, or by compari- 
son with estimates made from independent correlations with other groups of 
records. The standard error of annual estimates computed for only a few 
years of record obviously may be subject to considerable error, but it at 
least affords a crude measure of the probable magnitude of error if bias has 
been substantially minimized. Since relationships between annual figures 
can be established rapidly, this provides a convenient means for selection of 
groups of index records that probably will be the most reliable for extension 
of monthly records in the way proposed by the authors. 

This method of analysis was applied to the example given by the authors 
for estimation of the discharge of Murder Creek near Evergreen, Alabama 
from that of Sepulga River. As shown in Table 3, there was a better correla- 
tion between the annual runoff of Murder Creek and the average of two 
streams, Sepulga River and Escambia Creek, than between that of Murder 
Creek and Sepulga River alone. The Escambia Creek basin is on the other 
side of Murder Creek, about 30 miles to the southwest. (The indices for the 
Escambia Creek basin were multiplied by a constant factor, 0.88, to avoid 
giving greater weight to the variations in that basin than to the variations in 
the Sepulga River basin, because of the relative magnitude of the figures). 
Since the monthly distribution of the runoff in the three basins is similar, 
the comparison of the standard errors for the period of definition, 1941 to 
1945, suggests that a substantial improvement in the monthly correlation also 
could be obtained by use of average indices. The relationship between the 
monthly discharge of Murder Creek and the average of the other two streams 
(in amounts per square mile) in fact is nearly a linear one for the entire 
range, and the standard error for the sample of 60 values is only 0.05 log 
unit, or about 60 percent of the authors’ value. This improvement is roughly 
the same as between the two sets of annual estimates listed in Table 3, for 
the years 1941 to 1945. 

Statistical tests, whether of monthly or annual values, will not show 
whether the relationship for the period of definition is biased in relation to 
that for another period, as in the example of Figure 4. For that reason it 
would seem desirable to compare the results of monthly estimates with an- 
nual estimates based on the average of several index records whenever the 
monthly indices themselves cannot be averaged to minimize the possibility of 
bias. 

If it appears that estimates of annual runoff for a given stream are apt to 
be more reliable than estimates of monthly runoff, the latter could be multi- 
plied each year by the ratios of the annual estimates to the sum of the month- 
ly estimates. It would have been desirable to make such adjustments in the 
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estimates of Table 1, for example. This would have brought the estimates of 
mean runoff for the 4-year periods, 1935-1938 and 1942-1945, from within 
roughly 20 percent to within a few percent of the correct values. If the peri- 
ods happen to be critical ones from the standpoint of dependable water sup- 
plies or regulation by storage, such improvement would be of considerable 
importance. 

The graphical method proposed by the authors for estimating the relation- 
ship between monthly values of runoff is simple and can be applied rapidly. 
It is much more practical than numerical analyses, and considering the pos- 
sible magnitude of sampling error in short-term correlations, more elabo- 
rate methods of analysis usually would not seem appropriate. The writer 
agrees that extension of short-term records is of growing importance and 
that the suggested methods constitute a useful approach to a sound procedure 
for extending records of monthly discharge. It is believed, however, that the 
statistical measures should be viewed with caution as indices to the reliabili- 
ty of such estimates and that there is a need for independent estimates and 
controls of the kind proposed in this discussion. 
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Estimated runoff of Wading River near 
Norton, Mass., compared with recorded runoff. 


Runoff, 
Runoff, inches Error inches Error 
a in @ of b/ in gof 
Water Sum of monthly estimated Annual estimated 
Year estimates Recorded values estimates values 


+ 

~ 


13.4 12.6 13.5 
20.4 23.7 24.5 
15.5 15.1 17.8 
31.0 29.35 31.4 
21.2 23.2 23.1 
17.2 23.5 24.2 
25.3 24.4 
24.8 26.0 
33.8 33.7 
23.9 26.2 
22.0 22.4 
14.8 15.7 
16.1 15.3 
21.2 21.5 
11.8 13.1 
27.7 25.2 


21.8 22.4 


+ 
w 
w 


* 


N 
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a/ Monthly estimates were based on the runoff of Taunton River. 


b/ Annual estimates were based on arithmetical averages for 
Taunton River at State Farm, Mass., Charles River at Waltham, 
Mass., and Blackstone River at Woonsocket, R. I. 

The relationship of the indices to the runoff of Wading 
River was established for the period, 1946-50, by the method 
illustrated in Table 2. The Charles basin (227 sq. mi.) is 
about 20 miles north of Wading River, and the Blackstone 
basin (416 sq. mi.) is about 25 miles northeast. 


Figures of recorded runoff were taken from Water-Supply Paper 
1301 and rounded to the nearest tenth in this table. 


Table 1 
i 31 
32 
33 
34 
1935 
38 
39 
1940 
41 
42 
43 
44 
1945 4 
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Annual runoff of Wading River as 
estimated from relationship to 
indices of areal precipitation; 1946 to 1950 


All data in inches 


xX Ole Ol. 
Water R R R plus times minus R Error 
Year 25.0 1.00 25.0 ay 


26.5 e 28.6 53.6 53.6 28.6 29.6 | 
47 19.9 17.1 18.8 18.6 43.6 43.6 18.6 16.7 +1.9 . 
48 32.8. 25.7 26.9 27.1 82.2 29.6 -2.5 
49 22.5 15.4: 36:6 2642 48.2 4303 18.1 18.6 -0.5 
1950 12.4 138.6. -38.7 36,9 --36.9 13.9 11.6 +2.3 
Estimated loss: 125.0 125.0 
Estimated precip. 21.3 


Annual runoff (adjusted for diversions and storage ) 


Taunton R. at State Farm, Mass. 
Charles R. at Waltham, Mass. 


Blackstone R. nr. Woonsocket, R.I. 
: (R, +R, + 

R_: Wading R. nr. Norton, Mass. 

Ratio, b: 231.1/231.3 = 1.00 


Standard error of estimate: V Z ay” /(n - 1) = 2.02 inches 


4 Estimates made for the period 1950 to 1945, as in column 8, i 
| are listed in Table 1. The standard error during this 
] period was 1.40 inches. 
i Note: Columns 6, 7 and 8 are unnecessary in this example, - 
| but are shown to illustrate the usual procedure. 2 
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Table 3 


Annual runoff of Murder Creek near Evergreen, Ala., 
estimated from indices of areal precipitation. 


1 


2 


57. 

56.9 
68.4 
53.6 
66.2 
62.3 
60.5 
69.8 
49.3 


Average los 


49.9 50.4 
49.3 49.7 
65.5 62.9 
42.5 44.8 
68.0 63.2 
53.4 54.2 
50.8 52.1 
62.8 62.2 
37.4 40.4 


23.5 
22.9 
36.8 
17.6 
37.1 
27.6 
25.3 
36.1 
12 
2 


21.3 . 23.1 
22.6 22.5 
35.4 
19.9 
36.2 
28.4 
25.5 
33.2 
14,1 
20.4 


Sak dt 
*eeee 


ee 
Var 


8: 30 in./yr. (estimated) 


= 1.065, 


b= P = 1.062 
R+150)/ P,, 


where values of P and R are totals for period, 1941-45. 


1.065 


P - 30 


- 30 


= Annual runoff of Escambia Crk. + 30 inches 


= 0.88P 


1 


= Annual runoff of Sepulga River + 30 inches 


= Average of P_' and P 


1 2 


= Annual runoff of Murder Creek 


dy and dy' = Errors of computed values, Ry and Ry'. 


Standard error of estimate 


Estimated from deviations dy 


c es, 
1.69 inches, 1946-50 


Estimated from deviations, dy' 


nec es, 
4.35 inches, 1946-50 


| 
1041-28 HY 4 
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(All data in inches) 
oe ater P P P P d d 
1941 72.6 72.8 15.5 17.0 -1.5 15.5 
42 51.0 
Ee 1945 47.3 
46 58.4 
48 53.4 
49 61.5 
ean : 
| 
| 
a 
1 
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LEGEND 


Circled points: Runoff for 
indicated years ending 
September 30. 

Solid circles: For period 
used in definition of 
monthly relation. 

Dashed line: The average 
of annual totals (shown 
by crosses) as derived 
from monthly estimates. 
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RUNOF?, INCHES, TAUNTUN R, AT STATS FARM, MASS. 


FIGUKZ & CORKELATION OF ANNUAL RUNOFF 
October 1929 to September 1950 


ASCE 
i 
3 
Ng 
| 


. 
1 
id 
4 
4 
: : 
; 
1 
| 
+ 


ASCE 1041-31 


Discussion of 
“FLOOD PLAIN ASPECTS OF RIVER PLANNING” 


by Anthony M. Lunetta 
(Proc. Paper 1040) 


ROBERT L. SMITH,* A.M. ASCE.—The role of the flood plain in river 
planning is seldom limited to a specific type of water problem, nor is it al- 
ways limited to consideration of just the water resource. Under certain con- 
ditions problems of land, wildlife, or mineral resource management can be 
expected to influence materially the planning decisions governing water re- 
source development. By contrast the remarks in the foregoing paper have 
been limited to a specific water problem—floods. Likewise the remarks on 
the possible ways to minimize the damages caused by those floods were 
limited to the proposal that the use of land within the flood plain be subject 
to regulation. Further, those remarks were summarized in a very logical 
sequence and they set forth quite adequately the factors to be considered. 
Notwithstanding those various limitations, the most significant feature de- eae 
veloped in those remarks was the acknowledgment of the breadth and, in a 
essence, the complexity of the problem discussed. Of particular importance 
is the concept that this problem of wise utilization of the flood plain is inter- 
disciplinary in nature. In fact, engineering consideration of the problem can- 
not dictate the answer, but can only serve to crystallize the solution to 
fundamental economic, legal, and social questions. 

The subject presented is important. In every area of the country, in areas 
x both developed and underdeveloped, day by day expansion is continually in- ta 
a creasing the potential flood damage. The author suggests that this expansion ‘ 
is twofold. First, there is encroachment which is motivated by economic de- 
sire to utilize the resources of the flood plain. Secondly there is encroach- 
ment which did not seek the flood piain for a specific economic purpose, but 
which, because of ignorance of the risk involved, blissfully located in the 
green valley. It must be acknowledged, however, that all investments under 
our free enterprise system acknowledge and allow successfully for various 
risks, and it is submitted that in the final analysis either inability or refusal 
to evaluate the risk is responsible for all unwise encroachment. Thus, there 
is presented the need for educating John Q. Public on the virtues of determin- 
ing his calculated risk. ; 

Engineers and laymen confronted with flood plain losses can agree on a :. 
common objective, namely—“The reduction of flood damages by more intelli- _— 
gent use of the flood plains.” In the paper just presented it is implied that = 
one way of attaining this objective is through the utilization of land use regu- oe 
lations. It remains to be seen whether engineers, the courts, and the public 2 
will be willing to adopt a regulatory type approach. Four reasons were pre- : 
sented why widescale adoption of regulatory procedures has not been realized. 
Those reasons were briefly (1) lack of hydraulic data, (2) lack of land use 


* Director, lowa Natural Resources Council, Des Moines, Iowa. 
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data, (3) objections by property owners, and (4) lack of legal precedent. It is | 
submitted that these four points again illustrate the present inability of indi- ip 
viduals, courts, developers, and others to evaluate properly the risk involved ts 
in certain flood plain development. 

A secondary primary problem to be considered is tightly interwoven 
among the various disciplines. Simply stated the problem becomes one of 
defining what is the right of the public and what is the right of the individual. 

Primary arguments for public control hinge on protection of the general wel- 

fare from the direct adverse physical effect of unwise development, and the 
secondary effect of modern floods on the community and the regional economy me 
and public health. Countering these public needs is the desire for preserving - 
greatest possible individual freedoms. The only adverse effects sustained 
from many encroachments is the damage to the facility itself. If that same | 
facility can be shown to represent a long term economic gain can it be clas- 
sified as a non-conforming use? Further, floods come in many shapes and 
sizes—there is the extraordinary and the routine. Does the same public 
right hold for all possible areas of inundation? The answer to these and “4 
other questions go far beyond the scope of this discussion. They will not be 5. 
answered succinctly or quickly. Their answers will, at best, evolve with ; 
time for they are simply a portion of a more diverse question of local, state, = 
and national policy. A question which the Engineers Joint Council recognized 
when they stated “The application of desirable water policies can only be car- 
ried out . . . . when an overall policy of land utilization is available,” (1) 
Paradoxically, development of such a policy for flood plain areas must await 
a reasonable attainment in public understanding of the problem. 

The foregoing comments implying the need for an educational approach to 
the problem of flood plain use have not been asserted in opposition to the Za 
author’s implied plea for positive zoning procedures. However, the educa- 
tional approach has been grossly overlooked at the local, state, and national 
level. Presumably the responsibility for development of adequate education- 
al program could be a proper function of government. 

The writer once suggested the desirability of interpreting and presenting 
available hydraulic data in the form of river plan and profile sheets.(2) Ade- 
quate strip maps can be developed in many parts of the country from avail- . 
able aerial photos. Sufficient information such as river mileages, section ms 
centers, and major highways would have to be identified to allow correlation 
of plan and profile. The profile sheets might well include a low water profile, 
profiles for floods of various percent chance of occurrence, and a profile for 
a so-called provisional or project design flood. Seldom could a stream be 
found where the available data seemed sufficient in every detail. However, 
treated as a design problem the task would reduce to certain standardizing 
procedures. For example, recent utilization of regional analyses of stream 
flow data could be a useful tool for such work.(3) It seems apparent that such 
data would be in much demand once investment houses and others realized the 
significance of the message it contains. Particularly would this be true 
where those who seek flood plain development are required to pay their pro- 
portionate share, based on benefits received, of the cost of future protection 
works. 

Finally one cannot lose sight of the importance of the author’s last state- 
ment, Bigger floods are yet to come. Zoning or less formal approaches to q 
land use regulation can only be effective up to a certain limit. Thus land use 
regulations cannot be advocated as a panacea or complete cure, they simply 
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represent another tool for minimizing flood damage. If the thoughts ex- 
pressed in the paper are sufficient to reawaken in the minds of only a few the 
possibility of utilizing that tool, the paper will have served its purpose. 
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Paper 1052 
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HYDRAULICS DIVISION 


Proceedings of the American Society of Civil Engineers 


COAGULATION AND SEDIMENTATION@ 


J. M. DallaValle” 
(Proc. Paper 1052) 


SUMMARY 


By assuming initially a monodisperse system of particles or nuclei, the 
chance of collisions leading to 1, 2, 3, etc. aggregated particles is computed 
according to the procedure used by Smoluchowski. The effect of settling is 
ignored by assuming the process is very fast. A discussion is included of the 
general law of settling of dilute suspensions as well as highly concentrated 
suspensions. 


The breakup of matter into finer and finer particles can proceed only to a 
certain point before the phenomena of coagulation or coalesence begin to 
bring them together again. It may well be that the forces tending to disinte- 
grate matter are eventually counterbalanced by others which act to bring par- 
ticles together. In industry, fine particles are made to join together by a 
form of coagulation called polymerization and this is the basis of plastics 
manufacture, rubber synthesis and many other products. Often it is a transi- 
tory phenomena on the way to making a final product and thus receives little 
attention. Or, as in the case of water purification, it is a carefully controlled 
phenomena. Coagulation is not a phenomena characteristic of liquids alone. 
The manufacture of metallic oxides by burning pure metals exhibits this 
phenomena. In fact, with many metallic oxides, particle size control is 
achieved by letting the particles “grow”—in other words, to aggregate or 
coagulate while still in suspension. ** 

The tendency of particles to coalesce or coagulate is far more widespread 
than supposed. There can be no question but they act in all fine particle 


Note: Discussion open until January 1, 1957, Paper 1052 is part of the copyrighted 
Journal of the Hydraulics Division of the American Society of Civil Engineers, 
Vol, 82, No, HY4, January, 1956, 

a. Presented at the Dallas, Texas, Meeting of the ASCE, Feb. 13-17, 1956. 

b. Georgia Inst. of Technology, School of Chemical Eng., Atlanta, Ga. 

** We shall use the terms coagulation, coalesence, aggregation, cohesion and 
accretion synonymously and understand that they apply to the growth of 
particles by accumulation to their surfaces of other particles of the same 
species. 
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suspensions and are essential to the formation of stable sedimentary de- 
posits. The constant stirring up of tidal basins, channels and shorelines by 
wave motion, of seasonal upheavals of lake bottoms by thermal inversions 
and even atmospheric dusts as they are swept about by winds would make 
this a difficult world in which to live in the absence of coagulation. 

This paper proposes to outline the mechanism of coagulation and the bind- 
ing forces between particles once they stick together, and then to show how 
particles settle to form the tremendous sedimentary deposits which consti- 
tute so large a part of the earth’s crust. 


Smoluchowski’s Theory of Coagulation 


It is noteworthy that since the fundamental investigations of Hardy (1900), 
Wiegner (1911), Smoluchowski (1918) and more recently, Muller (1926, 1928), 
Tuorila (1927) and Fuchs (1934), the basic theory of coagulation has been only 
slightly modified. It is especially to Smoluchowski that we owe our present 
appraisal of the mechanism of coagulation. Muller and Fuchs have extended 
certain aspects of Smoluchowski’s theory, but for all its simplifications rec- 
tified in many respects by later investigators, Smoluchowski’s basic concepts 
remain unchanged. His theory is applicable to all fluids—liquids and gases— 
and his equations have more universal applications with regard to coagulation 
phenomena than any of those derived with greater attention todetail. Consider- 
ing the many extraneous complications affecting coagulations suchas (in liquids) 
the electrolyte atmosphere and accompanying double-layer effects, and (in 
gases) the existance of electrostatic charges, it is indeed remarkable that 
Smoluchowski was able to develop so successful a theory; for these factors 
affect the mathematical treatment so as to render a solution of the problem 
impossible. 

Smoluchowski used a very simple model for deriving his coagulation equa- 
tions. He assumed that the particles were initially uniform in size and ran- 
domly distributed. He regarded them as subject to Brownian motion during 
the whole period of coagulation. He assumed that there were collisions be- 
tween the particles as a result of their random motions and that on collision 
the particles formed a single particle. Rapid coagulation was regarded as a 
process in which particles on every contact adhered to each other on every 
contact; slow coagulation was attributed to the fact that only a fraction of the 
particles adhered on collision. We shall be concerned only with the first of 
these two types of coagulation in this paper. Why particles should adhere to 
each other on first contact in some instances and not in others will be dis- 
cussed briefly later. 

As indicated or implied above, coagulation is that process which leads to 
the growth of particles in a fluid atmosphere by accretion of other particles. 
If initially there are no particles per unit volume and no settling occurs, at 
atime t after coagulation begins, there will be no change in the mass of 
particles suspended, but there will be some kind of a distribution consisting 
of single particles, double particles, triple particles, etc. If the particles 
remain small even after the coagulation process has proceeded for some 
time, or if the process is very rapid, we may ignore the effect of gravity. 
However, with time, all particles come under the influence of the gravitation- 
al field and the suspension will no longer consist of randomly distributed 
particles. In general, if the particle diameter is a function of time, D(t), the 
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density of the suspension at any point z (from the bottom of the container) 
and time t is, after coagulation has ceased 


De 
Blz,t) = +6, = j F(Dg) (Dg) (1) 


° 


where Po and ps are the densities of the fluid and particle, respectively, 
Co is the concentration of the suspension at the moment coagulation ceases, 
De is the effective Stokes diameter and F(De) its size distribution. In using 
Eq. (1), we may take t = 0 at the instant the rapid coagulation process is 
completed and the particles begin to settle. Both ¢ and C can, of course, be 
readily expressed on a number basis when it is convenient to do so. 

We shall let n denote the number of particles in a unit volume at time t 
after coagulation has started. The reciprocal of this number, 1/n, can be re- 
garded as the average space occupied by 2 particle. If we have some method 
of measuring the change of n with time, we shall find that a plot of 1/n vs 
t yields a straight line with a positive slope. Experimental methods are 
available for determining n from instant to instant in a coagulating system 
and for both sols and aerosols, provided gravity has not affected the random- 
ness of the distribution, data have been obtained which yield the simple but 
important equation 


- (1/n,) = kt (2) 


in which K is called the “coagulation constant.” Since the n’s are usually 
quite large, it follows that the constant K for a given system is very small. 
The constant K will, of course, depend on temperature and whether the sys- 
tem is quiescent or turbulent. 

Actually, Smoluchowski derived Eq. (2) using the following considerations: 
First he regarded the system to be initially monodisperse and randomly dis- 
tributed and consisting of an average of no particles per unit volume. The 
volume of all the solids together was assumed negligible in comparison with 
the volume of the liquid. Next, with time increasing and with Brownian mo- 
tion, all or a stated fraction of particles on collision adhere to each other. In 
order to calculate the number of particles colliding with a particle (assumed 
fixed), we may consider the particle to be surrounded by an imaginary sphere 
of radius ro and apply the Fick diffusion equation giving the number of par- 
ticles n' diffusing through a surface. This is a procedure for a steady state, 
but it can be shown to be sufficiently accurate. The Fick equation is 


n' «= ) 


where D is the diffusion coefficient. We can integrate this equation subject 
to the conditions: 


for r = 00 


n=0 for 
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The last condition states that the particle concentration at ro is kept at 
zero, that is, the particle adheres. Integrating, subject to the conditions 
stated, the number of particles intercepted by the sphere is 


nt = (3) 


By way of simplification, we assume ro and D to be the same for all par- 
ticles. If each particle acts as a center at which coagulation takes place 
within a sphere of influence of radius ro, then 


(1/2) = 2nrgn*D 


since two particles are involved to make a single particle. It follows from 
Eq. (3) that 


2 
nt «= = - (4) 


an equation which is integrable if ro and D can be shown independent of 
time. At this point, Smoluchowski proceeds along lines well reasoned in the 
kinetic theory of gases. Since diffusivity is an additive property, and since, 
because of size differences as coagulation proceeds, we have in the case of 
two particles, diffusivities Dj and Dg, then we may write 


Also, on the basis of kinetic theory of gases, the new sphere of influence of 
two coagulated particles will be 


Po = (1/2) (Fo) + 


in which roj and rog are the radii of influence of the two particles before 
coagulation, We note that if ro, and rog are equal to twice the radii of the 
particles, the particles touch each other sufficiently close to pull together and 
coagulate, If the radii of the sphere of influence are greater than twice the 
radii of the particles’ field of influence, then obviously the particles will al- 
ways coagulate, With these assumptions 


dn n(D, + Dy ) (roy + (5) 


at 
Einstein has shown that for particles in Brownian motion, the diffusion 
coefficient is 
D = 
where R is the gas constant, T the absolute temperature, N is Avogadro’s 
number and t@ the mobility defined by Stokes law of resistance 
a = 1/3mD, 
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ut being the viscosity of the fluid. Thus, 


D = (RT/N) (1/3mD,) 


and from Eq. (5), if Dey and Deg are the actual effective diameters of the 
coagulating particles, we get 


- = (RT/3.") + * Toe yn? 


Now again, as a further simplification, we may say that the radii of influence 
are as their respective diameters—that is, roi = aDejy and ro2 = aDe2, a 
being a constant of proportionality. The last equation then becomes 


If the sol or suspension is initially regarded as monodisperse and if as time 
proceeds the differences in sizes is not very great, then we may put roj = 
To2 = etc., and we obtain an expression which is independent of particle size 


- (1fmg) © (20/3) (RT/WN) >t = Kt (7) 
K = (2/3) (RTa/uN) 


where the limits of integration are obviously n=ng when t=0 and n=n 
when t=t. Equation (7) is identical to Eq. (2). However, it includes the 
equation of state and an unknown constant a which must be evaluated by ex- 
periment. It is indeed curious that Eq. (7) with all the simplifications and 
assumptions made in deriving it should hold even for the early stages of 
coagulation. But there is much evidence that it applies very well to many 
coagulating systems until gravity begins to take over. 

Muller (1926, 1928) has re-examined the general theory and has shown 
how it can be applied to a system consisting of spherical particles of two dif- 
ferent sizes. He has also indicated the modifications needed for dealing with 
coagulation of needle-like particles. Muller’s treatment of these special 
cases are quite involved. 

Fuchs (1934) appears to have been the last investigator to tackle the the- 
oretical aspects of coagulation. However, he dealt with aerosols and by in- 
troducing concepts well-known in quantum theory advanced our knowledge of 
coagulation a great deal. His results, however, while more precise, remain 
to be tested. 

Equation (7) can also be applied to aerosols and this has been done with 
considerable thoroughness by Whytlaw-Gray (1932). For such applications, 
it is necessary to introduce a correction in mobility 1 as defined by Stokes 
equation, which takes account of particle slippage between gaseous molecules. 
In this case 


= (1 + 
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where A is the mean free path of the gas molecules, A is a parameter de- 
pending on the ratio 4 /De and other terms are as previously defined. The 
correction applied above is known as Cunningham’s correction and a discus- 
sion of its significance as well as values of A are contained in most texts 
on physical or colloid chemistry. 


Forces of Attraction 


Any discussion of the forces terding to hold particles together leads to the 
very nature of matter itself. The following discussion is intended to elucidate 
on mechanism, It is necessarily brief and intended only to give an idea of the 

magnitude of the binding forces. 

Let us assume that a particle in a coagulating system is momentarily fixed 
and is struck by another particle which has a kinetic energy of the order of 
10-8 ergs. This is a very small amount of energy easily acquired by a par- 
ticle 5 microns in size moving at a speed of about 10 cm per sec, But this 
amount of energy, small as it is, when given up in an inelastic collision with 
a cube of nickle 100 Angstroms on edge is sufficient to heat this cube to a 
temperature of about 100 deg C, (in air), Thus, if this energy were possessed 
by a particle of 100 Angstroms per edge and it struck a larger particie in- 
elastically, then it would appear that there is sufficient energy to cause the 
two particles to coalesce, Collisions in such instances “fuse” the particles 
together. 

But there are other forces of cohesion when particles come in contact that 
are referred to as van der Waal forces. These are forces exerted by unsatis- 
fied atoms (as might occur at surface discontinuities of a freshly made crys- 
tal) and which on contact with similar atoms in another particle tend to bind 
each other together. These binding forces are of great intensity and are of 
the order of several tons per square inch. Thus, the probability of very small 
particles joining together on contact is indeed very great. 


Particle Distribution Resulting from the Coagulation Process 


As coagulation proceeds, single particles will collide to form doublets. 
Single particles will collide with doublets to produce triplets, etc. Again 
doublets will be colliding with doublets to form quadruplets, or with triplets 
to form quintuplets, etc. Since Eq. (7) gives the law of diminution of particle 
numbers as coagulation proceeds, it is of interest to see whether it is possi- 
ble to determine the number of singletons, doublets, triplets, etc. present in 

a suspension at any time t after the process has begun. In this regard, we 
may again follow Smoluchowski. Let no as before be the number of particles 
initially present in a unit volume of the suspension, and n the number of par- 
ticles present at time t and composed of 


nj particles which are singles 
ng particles which are doublets 
ng particles which are triplets 
etc. 
nj particles which are i-lets 
i ny 
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In deriving Eq. (7), we have assumed that the coagulation constant K applies 
to all the particles. If each collision of two unit particles decreases the num- 
ber of units by 2, while collisions with other sizes diminishes the number by 
one, then 


2 
/at = + 2knjn, = 2Knjn 


For the decrease of doublets formed by coagulation of 2 unit particles 


| 2 
= “Kn,” + 2Knon, + 2kn, + -2Kn, + 2Knon 


and for triplets 
-dn;/dt = + 2Kn3n, + + + 2kn zn), tore 
-2Kn no 2Kn 3n 
Proceeding in this manner, we then obtain a series of expressions for i- 
particles. Since the rate of decrease of all the particles is controlled by 


eq. (7), we may use this equation to eliminate time 


- dn/ét « -2 dn, /at 


and thus obtain differential equations with nj, ng, ... , nj as variables. 
These equations are easily integrated and we obtain the geometrical sequence 


* (fn, - n)/n,? 


These equations, when added, form a geometric series with a common term- 
to-term ratio of (mo - n)/no. These are the most convenient form to work 
with, although (following Smoluchowski) the following arrangement for the 
general term is also useful: 


nm = n, (t/,)* ~ 4/2 + *2 


where we define 


t, * 


which is of course, the reciprocal of K with a = 2, a generally accepted 
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value of this constant. Figure 1 shows how the ratio n/ng changes with the 
ratio t/to with different values of nj as a parameter. We see, therefore, 
that if we have a means of determining the number concentration of particles 
at any time t and the initial concentration when t = 0, then the distribution 
of coagulated particles in terms of singletons, doublets, triplets, etc. is 
readily determined. 

We note in passing that the most frequent number of coagulated particles 
is given by the expression 


No = (2%, -n)f (10) 
or, in terms of the coagulation constant K and the time t 


N = 2Kn.t + 1 (11) 


Note also as shown in Figure 1 that the maximum of each grouping shifts to- 
ward the right as time increases. 


Sedimentation 


When the effect of gravity takes over, the coagulated particles will settle— 
provided there is no further accretions—in accordance with well defined laws. 
At the instant gravity may be said to take over, the density of the suspension 
at any point z measured vertically in the suspension and time t is con- 
trolled by Eq. (1).* 


Free Settling 


We shall now recognize that there are two kinds of settling: (1) that form 
in which the particles in the act of settling do not interfere with each others 
motion, This is usually the case when the volume concentration of the sus- 
pension is less than five per cent of the total volume. (2) Hindered settling 
which occurs with thick suspensions and is such that the particles settle as a 
mass with an interface between clear fluid and suspension. 

In free settling, the downward velocity of each particle is in accordance 
with Stokes law 


= ~ = sft (12) 


where z is the distance moved by the particle without acceleration in the 
time t and other terms have already been defined, If now in Eq. (1) z is 
fixed, we obtain by differentiation of and application of Eq. (12) that 


at 2t (13) 


in which @ is the ratio of concentration of particles at z at time t and the 


* Throughout this section time t= 0 is taken at the moment coagulation 
ceases and gravity takes over. 
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concentration when the particles were uniformly distributed in the suspen- 
sion (time, t = 0), that is, C(t)/Co=@. HE, on the other hand t is constant 
and we wish to obtain the distribution at various levels in the suspension 


(14) 
az 22 


and in this case, Z refers to the ratio of concentrations at various z-levels 
to the initial concentration, C(z)/Co. Equations (13) and (14) serve to give a 
complete history of the free-settling phenomena. 


Hindered Settling 


In the case of hindered settling, the interface between the liquid and the 
suspension moves at a constant velocity ug. This problem may be handled 
in much the same was as flow through a packed bed of particles with the par- 
ticles moving instead of the liquid. In fact, the Kozeny-Fair-Hatch equation 
in the form proposed by L. Hatch (1943) may be used. This equation applied 
to an expanded bed of fine particles with a porosity «. The form we use 
here applies to slow settling, but Hatch indicates well the general form of the 
equation for other types of flow. The equation is 


(ag? /eDg) ((1 -€)/E3} = 5.56 x 107% - 


where Dg is the specific surface diameter obtained by permeametric meth- 
ods (DallaValle, 1948). Other terms in the equation are defined elsewhere in 
this paper. Note that as Eq. (15) is written that the left hand portion is a _ 
modified Froude group and the right hand side a modified Reynolds group. 
Hence in general, we may expect the following equation to govern hindered 
settling: 


Fr, = (Re), (16) 


This equation holds very well up to porosities of 0.85 and for particles rang- 
ing from 5 to more than 1000 microns in water, oil or glycol suspending 
media. In fact, it is a convenient method for measuring the characteristic 
diameter Ds of particles. However, with coagulated or aggregated particles 
which in settling carry along with them an adhered layer of fluid, density cor- 
rections are necessary. In order to correct for this adhered liquid, it is 
necessary to follow along the lines indicated by Steinour (1944). 


SUMMARY 


By assuming initially a monodisperse system of particles or nuclei, the 
chance of collisions leading to 1, 2, 3, etc. aggregated particles is computed 
according to the procedure used by Smoluchowski. The effect of gravitation 
is ignored by assuming the process is very fast. A discussion is included of 
the general law of settling of dilute suspensions as well as highly concen- 
trated suspensions, 
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Figure 1. 
Relation between number of particles and time in rapid coagulation. 
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1956-18 HY 4 


DIVISION ACTIVITIES 
HYDRAULICS DIVISION 
Proceedings of the American Society of Civil Engineers 


NEWS BULLETIN-August 1956 


The April issue listed the committees for the Society business year 1955- 
1956 and referred to the Official Register for the names of members serving 
on Task Forces, The following additions and corrections should be made in 
your Register to make the record of those at work in the Division complete: 


Task Force On Energy Dissipators for Spillways and Outlet Works: 
Add Rex Elder 


Task Force on Hydraulic Aspects of the Design and Operation of Spillways 
and Crest Gates: 
Alvin J, Peterka is now chairman 
Marvin J. Webster is now a member 
Messrs. Nalder and Brown resigned. 


Task Force on Nomenclature for Hydraulics: (Note the change in the 


name of this task force) 
C. A. Betts is also on this task force 


Task Force on Water Measurements: 
Add A. L. Jorissen and C. A. Wright 


Task Force on Sediment Distribution in Reservoirs: 
Samuel Shulits, chairman 
Paul C. Benedict 
Lloyd C. Fowler 
Carl R. Miller 
John B, Stall 


The Hydrology Committee has established a Task Force on Hydrologic 
Data. The Task Force is staffed by the following members: 
D. W. Van Tuyl, chairman 


C. O. Clark W. W. Wheeler Victor A. Koelzer 
L. L. Harold Robert R. Balmer Edward J. Schaefer 
H, E. Hudson Verne Alexander Clayton H. Hardison 


Walter G. Schulz 
C. C. McDonald, Coordinator 


A Task Force on Spillway Design Floods has been authorized. 
Harvey O. Banks is chairman 


Note: No. 1956-18 is part of the copyrighted Journal of the Hydraulics Division of the 
American Society of Civil Engineers, Vol. 82, HY 4, August, 1956. a 
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The Committee on Research has two task forces at work. The Task Force 
on Aerated Flow in Open Channels consists of: F 
Warren W. De Lapp, chairman 
William J. Bauer 
Donald Colgate 
Joseph P. Lawrence 
(Frank B, Campbell is Contact Member) 


The Task Force on Cavitation in Hydraulic Structures consists of: 


F. R. Brown, chairman 

R. A. Elder 

V. E. Johnson, Jr. 
M. B. McPherson 
A. J. Peterka 
(J. M. Robertson is Contact Member) i 


Mr. Thornton J. Corwin, Jr., past chairman of the Executive Committee, 
recently resigned from the committee for reasons of health, It is regretted 
exceedingly that illness has made it necessary for him to give up his work in 
the Division. A recent note from Thor indicates that he is making a satisfac- 
tory recovery and is now able to be at his office in the Pacific Gas and Elec- 
tric headquarters. 

Dr. Arthur T. Ippen has been appointed by the ASCE Board of Direction as B 
a member of the Executive Committee of the Hydraulics Division. This ap- po 
pointment is effective immediately for a term ending in October of 1960. Art , 
brings with him a wealth of experience in hydraulics and committee work 
from which the Division will undoubtedly profit. 


Committee and Task Force Meetings 


The following meetings have been held so far this year since October 1955: 


Executive Committee Denver November 7-8, 1955 : 
Committee on Tidal Hydraulics Washington, January 17, 1956 i 
D.C. 
Committee on Hydrology Dallas, Tex. February 13-14, 1956 > 
Committee on Design Knoxville, May 3-4, 1956 x 
Tenn. 
Committee on Sedimentation . June 4-5, 1956 


Task Force on Sediment Distri- Urbana, Il. May 15-17, 1956 = 
bution in Reservoirs oa 


The Committee on Research will hold its next meeting at Chicago, 
August 25, 1956 


The Executive Committee will hold its next meeting at Madison, 
Wisconsin, August 21, 1956 
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HYDRAULICS DIVISION TO HAVE BIG PROGRAM 
AT PITTSBURGH CONVENTION 
October 15-19, 1956 


Place: William Penn Hotel 
Technical Sessions have been formulated by the following committees: 


Design 1 session 
Floods i session on flood insurance 
Research 1 session on turbulence in fluid flow 
1 session on cavitation in hydraulic 
structures 
1 general session 
Tidal Hydraulics 2 sessions sponsored jointly by the 
Committee on Tidal Hydraulics and 
the Committee on Coastal Engineer- 
ing of the Waterways Division 


Committeemen and Task Force chairmen are invited to use the Hydrau- 
lics Division News Bulletin to: 


a) Insert outline of proposed work and invite suggestions thereon 
b) Insert resume of work or report and invite suggestions thereon 
c) Insert questions on specific problems on which information is required 


Committee on Publications 


Chairman H, G. Dewey reports that since the publication of JOURNAL No, 
HY 2, the Division has passed on 14 papers to the Manager of Technical Pub- 
lications for publication in future issues of the JOURNAL; eight (8) papers 
are being revised by the authors; and three (3) papers are under active re- 
view. This indicates another very active year in publication of Hydraulic 


papers. 


HYDRAULICS DIVISION MEMBERSHIP 
Questionnaire 


It has been several years since the membership of the Hydraulics Division 
has been canvassed to determine the interest and availability of the individual 
members for service in Division activities. All who are interested in par- 
ticipation in Division activities are urged to fill out the following question- 
naire and return to the Secretary of the Division. 
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QUESTIONNAIRE OF PARTICIPATION OF HYDRAULICS DIVISION 


Name 
Address 
Position Division Committee Service | , 
Government Research q 
School Publications 
Stevens or 
Field of Interest Hilgard Prize 
Other Member 
Design 
Hydrology Review Papers 
Research 
Other No —— 
Floods 
Preparation of Technical Paper Sedimentation atch 
Gr’d Water 
By self Snow, Ice and 
Design 
Edit Papers Hydrology 
Yes Research General 
No prs Research Special- 
ized (Indicate 
specialty) 
Comments and suggestions 


Subjects recommended for discussion at meetings 


Suggested Division meeting place 


Please return to Harold M. Martin, Secretary, Hydraulics Division, ASCE, 
Bureau of Reclamation, Denver Federal Center, Building 56, Denver 2, 
Colorado. 
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